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The  concepts  of  topology,  which  are  inherent  in 
any  sort  of  thinking  about  interacting  systems,  are 
almost  totally  unfamiliar  to  the  ordinarily 
educated  person.  Though  extraordinarily  simple, 
the  ideas  of  topology  seem  not  to  enter  into  the 
experience  of  any  but  the  most  specialized  lives. 
As  a consequence  people  tend  to  regard  such  ideas 
as  strange  and  difficult.  Thus,  in  inviting  them  to 
adopt  the  new  viewpoint  of  microstructology,  one 
is  met  with  a double  resistance,  arising  from  a 
disinclination  to  abandon  a familiar  way  of 
thinking  and  from  difficulty  in  comprehending 
the  new  way.  This  is  a particularly  difficult 
obstacle  to  overcome  and  it  is  not  likely  to  go 
away  of  its  own  accord.  Some  positive  action 
will  be  required  to  stablish  microstructology  in 
the  thinking  of  metallurgists. 
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The  purpose  of  this  research  is  to  obtain  information  on  both 
structural  and  physical  properties  of  sol-gel  silica  monoliths  as  the 
material  goes  through  densification,  and  to  seek  possible  correlations 
between  the  properties  measured. 

The  structural  properties  are  divided  into  two  groups:  the  metric 
properties  and  the  topological  properties.  The  evaluation  of  metric 
properties  was  based  mainly  upon  nitrogen  gas  adsorption  experiments. 
In  order  to  obtain  topological  properties  some  geometric  models  of  the 
pore  structure  of  silica  gels  are  proposed,  which  are  used  in  conjunction 
with  the  gas  adsorption  data.  The  gas  adsorption  results  show  that  for 
three  types  of  silica  gel  monoliths  the  average  pore  radii  of  12  A,  32  A 
and  81  A remain  fairly  constant  during  densification.  The  application 
of  topological  modeling  to  silica  gel  monoliths  yields  the  number  of 


xxx 


nodes  (Nv),  number  of  branches  (Bv),  number  of  separate  parts  (Pv) 
and  the  genus  (Gv)  of  the  pore  network.  As  densification  proceeds  Nv 
and  Bv  go  to  zero  as  the  material  becomes  fully  dense.  Pv  is  initially 
zero  and  increases  as  pores  start  closing  in  the  last  stage  of  densification, 
and  the  genus  approaches  zero.  The  topological  models  proposed  yield 
consistent  results;  the  scale  of  the  structure  obtained  from  the  models  is 
the  same  as  is  observed  with  high  resolution  transmission  electron 
microscopy  for  dried  silica  gel  monoliths. 

The  mechanical  properties  explored  were  the  following: 
microhardness,  flexural  strength,  flexural  modulus,  uniaxial 
compression  strength,  compression  modulus,  tensile  strength  (from 
Brazilian  test),  tensile  storage  modulus  and  shear  storage  modulus  (from 
dynamic  mechanical  analysis).  In  general,  the  strength  and  modulus  of 
the  silica  gel  monoliths  increase  continuously  as  densification  proceeds. 
Measurements  of  dimensional  change  show  that  the  silica  gels  exhibit  a 
higher  thermal  stability  as  densification  proceeds.  Silica  gel  monoliths 
with  8 1 A average  pore  radius  show  greater  thermal  stability  than  silica 
gels  with  12  A average  pore  radius. 

The  strength  of  silica  gel  monoliths  with  12  A average  pore  radius 
increases  relatively  slowly  when  densified  from  180  °C  to  about  800- 
900  °C.  The  increase  in  strength  becomes  more  pronounced  when  the 
densification  temperature  is  further  increased,  in  the  range  of  900  °C  to 
1150  °C.  The  low  temperature  region  (180-900  °C)  is  associated  with  a 
high  pore  connectivity,  while  the  second  region  (900  °C  to  1150  °C)  is 
associated  with  pore  closure,  which  causes  stiffening  of  the  structure. 
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CHAPTER  1 

THE  SOL-GEL  METHOD 


1.1  Introduction 

Silica  (Si02)  in  its  various  forms  is  the  most  abundant  component  of 
the  crust  of  the  earth.  It  is,  therefore,  no  surprise  that  mankind  has 
been  using  this  component  ever  since  the  beginning  of  civilization.  The 
first  known  silica  based  glasses  are  dated  to  be  4,000  to  5,000  years  old 
[1,2].  Throughout  history  Si02  has  found  countless  applications  either  in 
its  crystalline  or  amorphous  forms.  SiC>2  has  been  widely  used  as  a 
major  component  of  most  glasses,  while  its  most  common  crystalline 
form,  naturally  occurring  quartz,  has  many  important  applications  such 
as  in  piezoelectric  transducers. 

In  glasses  SiC>2  has  an  important  characteristic  of  glass  forming 
because  of  the  tetrahedral  coordination  of  oxygen  atoms  around  the 
silicon  atom.  Extensive  studies  have  been  done  on  the  structure  of  Si02 
glasses  [3-23].  Figure  1.1  is  a schematic  of  a structural  segment  of  a 
pure  silica  glass.  Silica  glasses  can  be  divided  into  two  categories:  (1) 
multicomponent  silicate  glasses,  where  SiC>2  is  an  important  constituent 
but  there  are  more  elements  present;  and,  (2)  pure  SiC>2  glasses,  where 
Si02  is  the  only  component.  Applications  for  both  types  of  silica  glasses 
abound  [1,2,24-28]  and  use  of  pure  SiC>2  glasses  has  been  markedly 
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Figure  1.1  Tetrahedral  coordination  of  SiOi- 
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growing  recently  because  of  improved  processing  and  increased  ability 
to  control  the  final  purity  of  the  glasses. 

Traditional  processing  of  Si02  glasses  involve  the  melting  of 
naturally  occurring  quartz,  the  silica  glasses  resulting  from  this  process 
are  called  Types  I and  II,  according  to  Table  1.1  [29,30].  Type  I is 
obtained  by  electric  melting  of  natural  quartz  crystals  under  vacuum, 
while  Type  II  is  made  by  flame  fusion  of  natural  quartz.  Both  Types  I 
and  II  are  usually  called  fused  quartz.  According  to  Hench  et  al.  [30, 
P-77] 

deficiencies  of  fused  quartz  optics  can  include:  considerable  cation 
impurities  (Type  I),  hydroxyl  impurities  (Type  II),  inhomogeneities, 
seeds,  bubbles,  inclusions,  and  microcrystallites  [31].  The  relative 
extent  of  these  defects  depends  on  the  grade  of  fused  quartz  with  a 
corresponding  increase  in  cost  for  the  better  grades. 

Both  Types  III  and  IV  are  termed  synthetic  fused  silica.  The  cation 

impurities  in  Types  m and  IV  are  substantially  lower  than  in  Types  I 

and  II,  because  of  the  high  purity  of  SiCU  used  in  the  processing  of 

Types  III  and  IV  [30].  However,  water  contents  of  several  thousand 

ppm  are  present  in  Type  III  silica,  and  Cl  ion  in  a few  hundred  ppm  is 

retained  as  an  unreacted  residual  in  both  Types  III  and  IV  silica  [30]. 

1.2  Objectives 

Sol-gel  derived  colloidal  silica  systems  have  been  extensively  studied 
in  many  different  ways  [28,32-34].  Alkoxide  systems  by  other  hand  still 
lack  understanding  in  many  fundamental  aspects  of  processing. 
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Table  1.1  Silica-glass  processes  [29,30]. 

Fused  Quartz 

Type  I:  Electric  melting  of  natural  quartz  crystals 
heat 

Si02  (quartz)  > SiC>2  (glass) 

vacuum 

Type  II:  Flame  fusion  of  natural  quartz  crystals 
heat 

Si02  (quartz)  > Si02  (glass) 

oxy-hydrogen 

Fused  Silica 

Type  IQ:  Vapor  - phase  hydrolysis  of  pure  silicon  tetrachloride 
carried  out  in  a flame 

SiCl4  + 02  + 2H2  — > Si02  (glass)  + 4HC1 

Type  IV:  Oxidation  of  pure  silicon  tetrachloride  which  is  fused 
electrically  or  by  means  of  a plasma 

SiCl4  + 02  — > Si02  (glass)  + 2012 


Gel-Silica 

Type  V:  Gelation  of  alkali  silicate  colloidal  solutions  with  full 
densification  (1500  °C  to  1720  °C),  or 
Hydrolysis  and  condensation  of  alkoxide  precursor  with  full 
densification  (1 150  °C  to  1200  °C) 


Si(OR)4  + 4H20  — > Si(OH)2  + 4ROH 

Si(OR)4  — > Si02  (dense  gel-silica)  + 2H20 

Type  VI:  Hydrolysis  and  condensation  of  alkoxide  precursor  with  partial 
densification  (600  °C  to  950  °C) 

Si(OR)4  + 4H20  — > Si(OH)4  + 4ROH 

Si(OH)4  — > Si02  (porous  gel-silica)  + 2H2O 
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Colloidal  systems  generally  involve  larger,  sub-micrometer  scale 
structures,  while  alkoxide  systems  exhibit  smaller  nanometer  scale 
structures.  It  is  the  purpose  of  this  work  to  study  the  structural  and 
physical  evolution  of  nanometer  scale  alkoxide  derived  silica-gel 
monoliths  as  densification  proceeds. 

The  purpose  of  this  chapter  is  to  review  the  sol-gel  method, 
concerning  the  stages  of  mixing,  casting,  gelation,  agingand  drying. 
Chapter  2 discusses  densification  and  is  intended  to  describe  the 
experimental  techniques  used  in  this  study  of  sol-gel  silica  processing. 
Analyses  of  nitrogen  gas  adsorption  and  evolution  of  structural  metric 
properties  for  gel-silicas  with  different  pore  radii  are  reported  in 
Chapter  2. 

Chapter  3 is  concerned  with  the  concepts  of  structural 
characterization  and  topological  modeling  of  the  nanometer  scale 
structure  of  silica  gel  monoliths,  and  describes  the  topological  evolution 
of  the  structure  during  densification.  Data  on  the  evolution  of 
mechanical  properties  with  densification  are  reported  in  Chapter  4.  The 
mechanical  tests  under  consideration  involve  microhardness,  flexural 
testing,  uniaxial  compression  testing,  Brazilian  testing  and  dynamic 
mechanical  analysis.  Chapter  5 describes  experimental  work  on  thermal 
properties,  with  emphasis  on  dimensional  change  and  its  correlation 
with  the  sintering  process.  Chapter  6 ("structology")  analyses  and 
correlates  the  structural  data  (metric  and  topological)  with  the  physical 
properties,  with  the  goal  of  obtaining  additional  understanding  of  the 
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the  structure  of  sol-gel  silica  with  densification.  The  conclusions  of  this 
work  are  presented  in  Chapter  7. 

1.3  Sol-Gel  Technology 
1.3.1  Evolution  of  Sol-Gel  Technology 

Sol-gel  research  has  grown  markedly  due  to  potentially  unique 
applications  for  this  new  class  of  materials  [35-43].  The  beginning  of 
interest  in  sol-gel  derived  materials  dates  back  to  mid  1800s  with 
Ebelman  [44,45]  and  Graham’s  [46]  studies  on  silica  gels.  For  a period 
from  the  late  1800s  through  the  1920s  gels  received  special  attention  of 
chemists  because  of  the  phenomenon  of  Liesegang  Rings  [47-53]. 
Geffcken  and  Berger  [54]  mention  sol-gel  dip  coatings  in  the  literature 
in  1939.  In  the  1950s  and  1960s  Roy  and  co-workers  [55-58] 
recognized  the  potential  for  obtaining  high  levels  of  chemical 
homogeneity  in  colloidal  gels.  Roy  and  co-workers  [55-58]  synthesized 
a large  number  of  novel  ceramic  oxide  compositions  (involving  Al,  Si, 
Ti,  Zr,  etc.),  that  could  not  be  made  using  traditional  ceramic  powder 
methods  [28].  During  the  same  period,  Iler’s  [32]  pioneering  work  in 
silica  chemistry  led  to  the  commercial  development  of  colloidal  silica 
powders,  Dupont's  colloidal  Ludox®  spheres  [28].  Stober  and  co- 
workers [59]  extended  Iler's  findings  to  show  that  using  ammonia  as  a 
catalyst  for  the  TEOS  hydrolysis  reaction  one  could  control  both  the 
morphology  and  the  size  of  the  powders,  yielding  the  so  called  Stober 
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spherical  silica  powder  [28].  The  final  size  of  the  silica  powder  is  a 
function  of  the  initial  concentration  of  water  and  ammonia,  the  type  of 
silicon  alkoxide  (methyl,  ethyl,  butyl,  pentyl  esters,  etc.)  and  alcohol 
(methyl,  ethyl,  butyl,  pentyl)  mixture  used  [28,59]  and  reactant 
temperature  [60].  In  the  late  1960s  and  early  1970s  the  sol-gel  synthesis 
of  well  defined  multicomponent  oxides  was  reported  [57,61-65]. 

Nucleation  of  particles  in  a very  short  time  followed  by  growth 
without  supersaturation  yields  monodispersed  colloidal  oxide  particles 
[66-68].  Using  these  concepts  Matijevic  and  co-workers  [69-72]  have 
produced  an  enormous  range  of  colloidal  powders  with  controlled  size 
and  morphologies  including  [28]  oxides  (TiC>2,  aFe203,  Fe3C>4, 
BaTiC>3,  CeC>2);  hydroxides  (AlOOH,  FeOOH,  Cr(OH)3);  carbonates 
(Cd(0H)C03,  Ce20(C03)2,  Ce(III)/Y  HC03);  sulfides  (CdS,  ZnS); 
metals  (Fe  III,  Ni,  Co),  and  various  mixed  phases  or  composites  (Ni, 
Co,  Sr  ferrites),  sulfides  (Zn,  CdS),  (Pb,  Cds),  and  coated  particles 
(Fe304  with  Al(OH)3  or  Cr(OH)3). 

The  controlled  hydrolysis  of  alkoxides  has  also  been  used  to  produce 
submicrometer  TiC>2  [73],  doped  Ti02  [74],  Z1O2  [75],  doped  Zr02 
[75],  doped  SiC>2  [76],  SrTiCU  [77],  and  corderite  [77]  powders. 

Emulsions  have  been  employed  to  produce  spherical  powders  of 
mixed  cation  oxides,  such  as  yttrium  aluminum  garnets  (YAG)  and 
many  other  systems  [28,77]. 

Sol-gel  powder  processes  have  been  applied  to  fissile  elements  [78] 
where  spray  forms  sols  of  UO2,  and  rigid  gel  spheres  of  U0-Pu02  are 
formed  during  passage  through  a column  of  heated  liquid  [28]. 
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Both  crystalline  and  vitreous  ceramic  fibers  have  been  prepared 
using  the  sol-gel  method.  Compositions  include  Ti02-Si02  and  Z1O2- 
SiC>2  glass  fibers  [79],  high  purity  SiC>2  waveguide  fibers  [80-85], 
AI2O3,  Z1O2,  Th02,  MgO,  Ti02,  ZrSi04,  3A103-2Si02  fibers  [86-91]. 
Alumina  grains  sol-gel  derived  are  important  commercial  products  [91]. 

A variety  of  coatings  and  films  have  also  been  developed  using  sol- 
gel  methods  [28].  Of  particular  importance  are  the  antireflection 
coatings  of  indium-tin-oxide  (ITO),  and  related  compositions  applied  to 
glass  window  panes  to  improve  insulation  characteristics  [39,92,93]. 
Other  work  on  sol-gel  coatings  is  reviewed  by  Schroeder  [94], 
Mackenzie  [95,96]  and  Wenzel  [97]. 

1.3.2  Motivation  for  the  Sol-Gel  Method 

According  to  Hench  and  West  [28]  the  motivation  for  sol-gel 
processing  is  primarily  the  potentially  higher  purity  and  homogeneity 
and  the  lower  processing  temperatures  associated  with  sol-gels, 
compared  with  traditional  glass  melting  or  ceramic  powder  methods. 
Wenzel  [42]  lists  the  following  advantages  of  the  sol-gel  method: 

1)  High  purity  raw  materials  are  available.  This  is  specially 
useful  in  preparing  fiber  optic  preforms  with  very  low  (ppb 
range)  concentration  of  contaminants. 

2)  Sol-gel  methods  can  be  used  to  make  glasses  of  many 
compositions  and  glasses  with  cations  in  unusual  oxidation  states. 
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Oxidation  state  manipulation  opens  the  possibility  of  new  colored 
glasses. 

3)  Highly  homogeneous  multi-component  gels  and  glasses  may  be 
synthesized,  or  the  scale  of  this  homogeneity  may  be  coarsened 
voluntarily  to  yield  a block  copolymer  gel  [98,991.  Homogeneity 
is  required  for  the  production  of  optical  glass  and  fiber  optic 
preforms.  Control  of  structure  at  the  molecular  level  will  lead  to 
glasses  with  unusual  structures  and  to  new  glass-ceramics. 

4)  Dried  gels  can  be  made  in  a wide  range  of  densities,  surface  areas 
and  pore  sizes.  This  may  be  used  for  catalyst  design,  transparent 
insulation,  molecular  sieving  and  in  impregnation  applications. 

5)  The  sintering  of  sol-gel  derived  glasses  occurs  at  relatively  low 
temperatures.  This  opens  the  possibility  of  producing  refractory 
glasses  of  high  homogeneity  and  purity  [79]. 

6)  The  rheological  properties  of  sols  allow  the  formation  of  fibers, 
films,  and  composites,  by  using  techniques  such  as  spinning, 
dip-coating,  injection,  impregnation,  or  simple  mixing  and 
casting. 

Some  potential  applications  of  the  sol-gel  process  are  the  following: 

1)  Gel-silica  optics  and  porous  gel-silica  optics  [29,30,100,101]. 

2)  Catalyst  supports  [102-105]. 

3)  Fiber  optic  waveguides  [80,106]. 

4)  Thermal  insulation  [107]. 

5)  Fiber  reinforced  composites  [108-112]. 

6)  Optical  thin  films  [113-116]. 
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7)  Superionic  conductors  [117]. 

8)  Optical  fiber  coupling  [118]. 

9)  Thermal  application  fibers  [119]. 

10)  Structural  applications  [120]. 

11)  Low  thermal  expansion  glass  [29,121]. 

12)  Protective  thin  films  [105,113,122,123]. 

13)  Impregnated  glasses  [124-127]. 

14)  Novel  composition  materials  [79,128-130]. 

15)  Coatings  [113,114,131-134]. 

1.3.3  The  Sol-Gel  Process 

The  word  "sol"  describes  the  dispersion  of  colloids  in  liquids 
[96,135,136].  Colloids  are  described  as  solid  particles  with  diameters  in 
the  range  of  1 nm  to  100  nm  [136].  A gel  is  an  interconnected,  rigid 
network  of  sub-micrometer  dimensions.  The  term  "gel"  embraces  a 
diversity  of  combinations  of  substances  which  can  be  classified  in  four 
categories  according  to  Flory  [137]: 

1)  Well  ordered  lamellar  structures. 

2)  Covalent  polymeric  network;  completely  disordered. 

3)  Polymeric  networks  formed  through  physical  aggregation; 
predominantly  disordered. 

4)  Particulate  disordered  stmctures. 


According  to  Hench  and  West  [28]  three  approaches  are  used  to 
make  sol-gel  monoliths: 

Method  1 - Gelation  of  colloidal  powders. 

Method  2 - Hydrolysis  and  polycondensation  of  alkoxide 

precursors  followed  by  hypercritical  drying  of  gels. 

Method  3 - Hydrolysis  and  polycondensation  of  alkoxide 
precursors  followed  by  drying  under  ambient 
atmosphere. 

A silica  gel  may  be  formed  by  network  growth  from  an  array  of 
discrete  colloidal  particles  (Method  1),  or  by  formation  of  an 
interconnected  three-dimensional  (3-D)  network  by  the  simultaneous 
hydrolysis  and  polycondensation  of  an  organometallic  precursor 
(Methods  2 and  3)  [28].  When  the  pore  liquid  is  removed  as  a gas  phase 
from  the  interconnected  network  of  the  solid  gel  under  hypercritical 
conditions  (critical  point  of  drying,  Method  2),  the  network  does  not 
collapse  and  a low  density  material  called  "aerogel"  is  produced  [28]. 
Aerogels  can  have  void  volumes  as  large  as  98%,  and  densities  as  low  as 
0.08  g/cm3  [28,138,139]. 

When  the  pore  liquid  is  removed  at  or  near  ambient  pressure  by 
thermal  evaporation  (as  used  in  Methods  1 and  3),  the  material  is 
termed  a "xerogel"  [28].  If  the  pore  liquid  is  primarily  alcohol  based, 
the  material  is  called  an  "alcogel"  [28].  The  generic  term  "gel"  is 
usually  applied  to  either  xerogels  or  alcogels,  whereas  aerogels  are 
commonly  designated  as  such  [28]. 
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Figure  1 .2  shows  a schematic  of  a typical  sol-gel  (xerogel  or  alcogel) 
processing  sequence.  As  the  figure  shows,  the  processing  can  be  divided 
into  seven  steps,  as  follows: 

1)  Mixing. 

2)  Casting. 

3)  Gelation. 

4)  Aging. 

5)  Drying. 

6)  Stabilization. 

7)  Densification. 

In  the  next  sections,  steps  1 through  5 will  be  briefly  reviewed. 
Steps  6 and  7 will  be  reviewed  in  Chapter  2. 

1.3.4  Structure  of  Sols 

According  to  Brinker  and  Scherer  [140]  the  electrons  surrounding  a 
nucleus  do  not  constitute  a spatially  and  temporally  uniform  screen,  but 
every  atom  is  a fluctuating  dipole.  The  effect  of  the  dipoles  is  the 
creation  of  forces  of  attraction  between  atoms  [140-143]  which  are 
termed  van  der  Waals  (or  dispersion  forces).  The  van  der  Waals  forces 
result  from  three  types  of  interaction:  permanent  dipole-permanent 
dipole  (Keesom  forces),  permanent  dipole-induced  dipole  (Debye 
forces)  and  transitory  dipole-transitory  dipole  (London  forces)  [140]. 

The  attractive  dispersion  forces  are  approximately  additive,  so  that 
they  become  important  for  particles  of  colloidal  dimensions,  and  they 


Figure  1.2  Schematic  of  the  sol-gel  process. 
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tend  to  cause  flocculation  of  colloids  [140].  For  two  infinite  slabs 
separated  by  a distance  h,  the  attractive  potential  (Va)  is  given  by  [140] 


VA 


-A 

12  n h2  ’ 


(1.1) 


where  A is  the  Hamaker  constant. 

The  van  der  Waals  forces  have  been  investigated  for  silica  systems 
[144-147],  and  it  is  possible  that  they  play  little  or  no  role  in  the  silica 
system,  except  possibly  for  colloidal  particles  larger  than  100  nm  [33]. 

The  DLVO  theory  (named  after  Derjaguin,  Landau,  Verwey  and 
Overbeek)  describes  the  stabilization  of  colloids  by  electrostatic 
interactions  [33,140,148].  The  DLVO  theory  predicts  that  coagulation 
will  result  from  reduction  in  the  double  layer  repulsion,  which  varies 
depending  on  the  pH  and  type  and  concentration  of  electrolytes.  The 
effectiveness  of  a given  electrolyte  (counterions)  in  reducing  the 
electrostatic  repulsion  depends  strongly  on  the  valence  of  the  ion  [140]. 
The  empirical  Schulze-Hardy  rule  [140]  predicts  that  the  concentration 
of  electrolyte  required  to  cause  flocculation  (cf)  is  inversely 
proportional  to  the  sixth  power  of  the  valence  of  the  ion  [140]: 


C £3  (k  T)5 

cf  — \2  0 6 z6  ’ 


(1.2) 


where  C is  a constant; 

e is  the  dielectric  constant  of  the  solvent; 
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k is  the  Boltzmann's  constant; 

T is  the  absolute  temperature; 

A is  Hamaker  constant; 
e is  the  electron  charge; 
z is  the  valence  of  the  ion. 

According  to  Brinker  and  Scherer  [140],  silica  does  not  conform  to 
the  DLVO  theory,  because  Si02  is  apparently  stabilized  by  a layer  of 
adsorbed  water  that  prevents  coagulation. 

In  a sol  prepared  by  hydrolysis  of  a salt,  the  primary  particles  have 
diameters  of  3 nm  to  10  nm  and  are  nonporous  [140].  However 
polymers  that  grow  in  alcohol-water  solutions  do  not  form  dense 
particles,  but  are  fractal  objects  [140,149-151]. 

1.3.5  Mixing 

According  to  Hench  and  West  [28]  in  Method  1 a suspension  of 
colloidal  powders,  or  sol,  is  formed  by  mechanical  mixing  of  colloidal 
particles  in  water  at  a pH  that  prevents  precipitation,  as  discussed  by  Iler 
[32,33].  In  Methods  2 and  3 a liquid  alkoxide  precursor,  such  as 
Si(OR)4  (where  R is  CH3,  C2H5,  or  C3H7),  is  hydrolyzed  by  mixing 
with  water.  The  hydrolysis  reaction  of  TMOS  (tetramethylorthosilicate, 
Si(OCH3))  is  represented  as  follows: 
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CH3 

I 

0 

1 

H3C  - O - Si  - O - CH3  + 4H20 
I 

0 

1 

ch3 


OH 

I 

->  HO  - Si  - OH  + 4CH3OH 
I 

OH 

(1.3) 


The  hydrated  silica  tetrahedra  interact  in  a condensation  reaction 
forming  -Si-O-Si-  bonds,  as  shown  in  Eq.(1.4): 


OH  OH  OH  OH 

I I II 

HO  - Si  - OH  + HO  - Si  - OH  ->  HO  - Si  - O - Si  - OH  + H20 
I I II 

OH  OH  OH  OH  (1.4) 


Additional  linkage  of  -Si-OH  tetrahedra  is  referred  to  as  a 
polycondensation  reaction,  and  eventually  results  in  a Si02  network  as 
shown  in  Eq.(1.5)  [28]: 
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OH  OH 

I I 

HO  - Si  - O - Si  - OH  + 6Si(OH)4  — > 

I I 

OH  OH 

(1.5) 

OH  OH 

I H H I 

HO  - Si  - O O - Si  - OH 

I I 

OH  O O OH 

II  II 

HO  - Si  - O - Si  - O - Si  - O - Si  - OH  + 6H20 
II  II 

OH  O O OH 

I I 

HO  - Si  - O O - Si  - OH 

I H H I 

OH  OH 


The  -Si-O-Si-  bonds  are  termed  "siloxane"  bonds,  and  the  -Si-OH  bonds 
are  termed  "silanol"  bonds  [140]. 

The  hydrolysis  and  condensation  reactions  given  by  Eqs.(1.3)  and 
(1.4)  can  be  represented  in  a simpler  way  as  follows: 

1)  Hydrolysis: 


-Si-OR  + H20  — > -Si-OH  + ROH 


(1.6) 


2)  Condensation: 


2-Si-OH  — > -Si-O-Si-  + H20  (1.7) 

The  hydrolysis  reaction  rate  constant  is  kn,  and  the  condensation 
reaction  rate  constant  is  kc.  The  kinetics  of  hydrolysis  and 
condensation  is  more  complex  than  indicated  by  Eqs.(1.3)  and  (1.4) 
[28,121,140,152-168].  Many  species  are  present  in  the  solution,  and 
furthermore,  hydrolysis  and  condensation  occur  simultaneously  [28]. 

The  important  variables  that  control  the  kinetics  of  hydrolysis  and 
condensation  are  temperature;  nature  (acid  or  base)  of  electrolyte; 
concentration  of  electrolyte;  nature  of  solvent  and  type  of  precursor 
[28].  Pressure  also  influences  the  gelation  process  [154,157,159],  but  it 
is  usually  not  a processing  variable  [28]. 

Considering  the  hydrolysis  of  TEOS  (tetraethylorthosilicate, 
Si(OC2H5)4)  in  different  solvents,  Aelion  et  al.  [152]  showed  that  kn 
increases  linearly  with  the  concentration  of  H+  or  H30+  in  acidic  media 
and  with  the  concentration  of  OH-  in  basic  medium. 

According  to  Hench  and  West  [28],  the  nature  of  the  alkoxy  groups 
on  the  silicon  atom  also  influences  the  hydrolysis  kinetic  rates.  In 
general,  the  longer  and  the  bulkier  the  alkoxide  group,  the  slower  is  the 
rate  constant  [153,158,161]. 

There  are  few  data  on  the  kinetic  rates  for  condensation  reactions 
[28,157,165,169,170].  Among  the  problems  associated  with  the 
evaluation  of  kinetic  rates  one  finds  [28]:  (1)  the  order  of  the  reaction 


may  vary  with  time  [163];  (2)  there  are  several  hydrolysis  and 
condensation  reactions  possible,  each  having  its  own  rate  constant 
[156,170,171]. 

Orcel  et  al.  [160,165,166,172,173]  explored  the  effect  of  acid 
catalysis  and  formamide  (a  DCCA  or  drying  control  chemical  additive) 
on  the  kinetics  of  hydrolysis  and  polycondensation  of  a TMOS  silica 
system  using  29Si  NMR.  The  data  obtained  [160]  show  that  acid 
catalysis  increases  kn  by  a factor  of  2 and  that  formamide  decreases  kn 
and  slightly  increases  k(\ 

Si02  particles  depolymerize  in  an  acidic  medium  and  the  resulting 
monosilicic  acid  forms  a yellow  complex  with  Mo,  which  can  be 
measured  optically  [28].  The  depolymerization  rate  constant  (ko)  can 
be  related  to  the  particle  diameter  d,  through  the  following  empirical 
law  [28,174]: 

log  d = a + B log  ko,  (1.8) 

where  a and  b are  constants. 

Table  1.2  summarizes  Hench  and  Orcel's  [160]  data  of  the  effect  of 
DCCAs  and  the  presence  of  acidic  conditions  on  the  values  of  kn,  kc 
and  d of  different  silica-gel  solutions. 

Zerda  and  Hoang  [175]  determined  experimentally  the  presence  of  a 
pentavalent  transition  state  for  the  hydrolysis  of  Si(OR)4.  General 
expressions  for  the  hydrolysis  can  be  written  as  follows  [28,175]: 
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Table  1.2  Physico-chemical  characteristics  of  different 
gel  solutions  [160]. 


SAMPLE 

SW55 

SF25 

SF50 

SF23 

vol.  formamide 

0 

25 

50 

25 

vol.  formamide  + MeOH 

h2o 

DI 

DI 

DI 

pH=3 

103  kH  (1  moHh-1) 

12 

7 

2 

25 

kc  (1  mol-1  h_1) 

29 

31 

0 

6 

d (nm) 

2 

2.2 

2.5 

- 
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OH-  + -Si-OCH3  ._>  -Si-OCH-3  ->  -Si-OH  + OCH'3 

I 

OH 


(1.9) 


H+ 

I 

H3O+  + -Si-OCH3  — > -Si-OCH3  — > 

I 

OH  (1.10) 

-Si-OH  + CH3OH  + H+ 

Quantum  calculations  [176-189]  applied  to  exploration  of  hydrolysis 
reactions  provide  additional  understanding  of  the  process  and  predict 
the  existence  of  the  pentacoordinated  silica  species  in  the  reaction  path. 

1.3.6  Casting 

After  mixing  the  components  of  the  sol-gel  process,  a homogeneous 
solution  is  formed  upon  stirring.  Depending  on  the  types  of  reactions 
involved  and  the  kinetics  of  hydrolysis  and  condensation,  a specific 
mixing  time  must  be  selected  so  that  the  homogenization  is  optimum  and 
the  viscosity  of  the  solution  is  appropriate  for  casting. 

The  container  into  which  the  sol  is  poured  must  have  the  following 


characteristics: 


22 


1)  Appropriate  shape,  since  the  final  material  will  keep  the  shape  of 
the  container. 

2)  Appropriate  size,  considering  that  the  gel  will  shrink  throughout 
the  entire  process  and  only  stops  shrinking  when  it  becomes  fully 
dense. 

3)  Appropriate  composition,  so  that  there  will  be  no  reactions 
between  the  sol  and  the  internal  wall  of  the  container. 

4)  Appropriate  surface  quality,  which  assumes  special  importance  in 
the  production  of  optical  quality  materials. 

5)  Appropriate  cleanliness,  since  all  types  of  particles,  dust  and 
contaminants  are  undesirable  in  the  final  product  and  will  serve  as 
stress  risers  [100]. 

6)  Appropriate  thermal  properties,  compatible  with  the  required 
processing  temperatures  associate  with  chemical  reactions  and 
possible  heat  treatments  during  the  following  gelation  and  aging 
steps. 

7)  Interfacial  characteristics  that  prevent  formation  of  bubbles  on 
the  gel  surface  which  can  act  as  stress  risers  during  drying  [100]. 

1.3.7  Gelation 

The  hydrolysis  and  condensation  reactions  that  occur  in  the  sol  will 
eventually  lead  to  the  formation  of  a gel  caused  by  clustering  and 
linking  of  the  particles.  The  three-dimensional  network  thus  formed 
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will  contain  a continuous  liquid  phase  permeating  the  space  between  the 
solid  phase. 

The  gelation  point  (or  gelation  time,  tgei)  is  usually  defined  as  the 
point  at  which  the  clustering  and  linking  in  the  sol  forms  a three- 
dimensional  network  that  can  support  a stress  elastically  [28].  The  sol- 
to-gel  transformation  is  not  a sharp  and  well  defined  transformation, 
but  it  rather  occurs  continuously.  Therefore,  the  practical 
determination  of  the  gelation  time  is  subject  of  many  different 
approaches  [28, 33, 121,1 40, 1 65, 1 90-204] . 

Sacks  and  Sheu  [200]  use  the  viscoelastic  response  of  the  gel  as  a 
function  of  the  shear  rate  to  determine  tgei.  The  complex  shear  modulus 
(G)  at  a given  frequency  is  given  by  [28] 

G = G'(co)  + G"(co),  (1.11) 

where  G'  is  the  shear  storage  modulus; 

G"  is  the  shear  storage  modulus  loss; 
co  is  the  frequency. 

The  elastic  component  of  the  sol-gel  material  gives  rise  to  the  storage 
modulus  G',  while  the  loss  modulus  G"  comes  from  the  viscous 
component  [140].  The  relative  measure  of  the  viscous  energy  losses  to 
the  energy  stored  in  the  system  is  defined  as  the  loss  tangent  tan  5 [140]: 


(1.12) 


2 4 


Sacks  and  Sheu  [200,205]  pointed  out  that  the  gelation  time  (tgei) 
cannot  be  accurately  defined  by  a certain  value  of  viscosity  (rj),  since 
the  viscosity  depends  on  the  shear  rate.  They  suggest  that  tgei  would  be 
best  located  by  using  the  loss  tangent  tan  5 (say,  the  location  of  the 
maximum  or  a particular  value  of  tan  8). 

The  type  of  catalyst  used  in  a gelation  reaction  can  have  enormous 
effects  on  the  structure  of  the  gels  formed,  as  well  as  on  the  rapidity  of 
the  gelation  process  [198,206,207].  Table  1.3  [198]  shows  the  effect  of 
some  selected  catalysts  on  the  gelation  time  of  silica  gels  made  from 
TMOS  [140,208]. 

The  temperature  dependence  of  the  gelation  time  can  be  represented 
by  an  Arrhenius  type  equation  as  follows  [140]: 

In  (tgei)  = A + jpjT,  (1-13) 

where  A is  a constant; 

E is  the  activation  energy; 

R is  the  gas  constant. 

For  silica  gels,  E is  found  to  be  in  the  range  of  10  to  20  kcal/mole 
[190,202],  depending  on  the  catalyst  and  alkoxy  group  [140]. 

According  to  Brinker  and  Scherer  [140],  it  can  be  shown,  using 
nuclear  magnetic  resonance  (NMR),  that  the  silica  network  is  more 
highly  cross  linked  when  the  water/alkoxide  ratio  is  high  and  that  the 
structure  at  the  gel  point  (tgei)  is  highly  variable. 
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Table  1.3  Gelation  time  for  different  catalysts  [198]. 


Catalyst 

mole  catalyst 
mole  TEOS 

Initial  pH 

tgel  (h) 

HF 

0.05 

1.90 

12 

HOAc 

0.05 

3.70 

72 

HQ 

0.05 

<0.05 

92 

HNO3 

0.05 

<0.05 

100 

H2SO4 

0.05 

<0.05 

106 

NH4OH 

0.05 

9.95 

107 

HBr 

0.05 

0.20 

285 

HI 

0.05 

0.30 

400 

No  catalyst 

- 

5.00 

1000 
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Different  studies  [151,197,207,209-213]  show  that  acid  catalyzed  sols 
develop  a linear-like  structure  with  very  little  branching,  whereas  base 
catalyzed  systems  develop  highly  ramified  structures. 

According  to  Zallen  [214],  the  Flory-Stockmayer  theory  (also 
termed  classical  theory,  or  mean  field  theory)  is  the  standard  theory 
used  in  physical  chemistry  to  describe  the  sol-to-gel  phase  transition  and 
the  molecular  size  distributions  which  accompany  it.  The  basic 
assumption  of  the  Flory-Stockmayer  theory  [215-221]  is  that  gelation 
may  be  modeled  as  a branching  process.  The  Flory-Stockmayer  theory 
is  essentially  a theory  of  dendritic  polymerization  [214].  Figure  1.3 
shows  a schematic  of  a tri-functional  (z=3)  branched  molecule 
[214,215];  "z"  is  the  functionality  of  the  branch  point  (node)  of  the 
polymer. 

According  to  Zallen  [214],  the  Flory-Stockmayer  theory  is  successful 
in  predicting  a sharp  gel-point  and  the  critical  extent  of  the  reaction 
(cco)  at  which  the  gelation  occurs,  which  is  given  by 

etc  = -TTf.  (1.14) 


At  a given  time,  the  weight  fraction  of  gel  (wg)  is  given  by  [137,140] 


wg 


i - a -p)3 

- p3  • 


(1.15) 


where  p is  the  fraction  of  formed  bonds. 
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Figure  1.3  Schematic  of  a tri-functionally 
branched  molecule  [214,215]. 
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Thus  the  Flory-Stockmayer  theory  can  predict  the  distribution  of 
molecular  weights  at  the  gel  point.  However,  the  model  allows  for  no 
ring  formation,  leading  to  an  increasing  number  of  nodes  as  the  radius 
of  the  cluster  increases  [28]. 

According  to  Brinker  and  Scherer  [140],  the  percolation  theory  of 
gel  formation  does  not  exclude  the  formation  of  closed  loops,  while  it 
has  the  disadvantage  of  not  yielding  analytical  solutions  for  the 
percolation  threshold  and  the  size  distribution  of  polymers.  Reviews  of 
percolation  theory  have  been  written  by  Stauffer  et  al.  [222],  Zallen 
[214]  and  Brinker  and  Scherer  [140]. 

Using  the  number  of  bonds  surrounding  a site  as  equivalent  to  the 
functionality  of  the  monomer,  the  threshold  for  bond  percolation  (pc) 
can  be  represented  by  [140] 

Pc  - (7TT)-  a.16) 

According  to  Zarzycki  [223],  experimental  results  indicate  that  for 
silica  sol-gel  systems  pc  exists  in  the  following  range: 

0.6  < pc  < 0.84.  (1.17) 

The  fractal  model  of  structures  as  designated  by  Mandelbrot 
[224,225]  has  been  applied  to  diffusion  of  particles  [226-228]  and  sol-gel 
particle  growth  [151,223,229-231].  The  nature  of  fractals  requires  that 
they  be  invariant  with  scale,  meaning  that  a fractal  is  required  to  look 
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similar  no  matter  what  level  of  detail  is  chosen.  A second  requirement 
for  fractals  is  that  their  density  decreases  with  size,  as  shown  in  Figure 
1.4  [28]. 

The  mass  of  a fractal  can  be  related  to  the  fractal  dimension  (dF)  and 
its  size  (radius  R)  by  [28] 

M = K RdF,  (1.18) 


where  K is  a constant. 

1.3.8  Aging 

According  to  Brinker  and  Scherer  [140],  the  processes  of  change  in 
the  structure  after  gelation  can  be  categorized  as  follows: 

1)  Polymerization. 

2)  Syneresis. 

3)  Coarsening. 

4)  Phase  transformation. 

Polymerization  increases  the  connectivity  of  the  network  by  means 
of  poly  condensation  reactions: 

-Si-OH  + HO-Si > -Si-O-Si-  + H20  (1.19) 

Zerda  et  al.  [232],  Ketls  et  al.  [233],  Smith  et  al.  [234]  and  Orcel  et 
al.  [235]  show  that  condensation  reactions  in  silica  gels  continues  after 
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Figure  1.4  Schematic  of  the  density  variation  with  size 
for  typical  crystals  and  fractals  [28]. 


gelation.  The  condensation  continues  to  occur  because  of  the  presence 
of  hydroxyl  groups  after  gelation.  The  rate  of  reaction  depends  upon 
the  temperature,  and  the  concentration  and  the  pH  of  the  solution  [140]. 
Because  of  the  creation  of  new  bridging  bonds,  the  continued 
polymerization  reactions  strengthen  and  stiffen  the  network  [140]. 
Aging  can  give  rise  to  hydrolysis  reactions  [140]: 

-Si-OR  + H20  — > -Si-OH  + ROH  (1.20) 

The  reverse  reaction  is  the  re-esterification  [140]: 

-Si-OH  + ROH  — > -Si-OR  + H20  (1.21) 

The  re-esterification  reaction  can  be  suppressed  by  using  an  excess  of 
water  [140]. 

Syneresis  is  the  shrinkage  of  the  gel  network  resulting  in  expulsion 
of  liquid  from  the  pores  [28,33,140,236,237].  Syneresis  is  generally 
attributed  to  the  formation  of  new  bonds  resulting  from  condensation 
reactions  or  hydrogen  bonding  [140].  In  most  inorganic  gel  systems 
syneresis  is  an  irreversible  process  [140]. 

Since  in  particulate  gels  the  structure  is  controlled  by  the  balance 
between  electrostatic  repulsion  and  attractive  van  der  Waals  forces,  the 
extent  of  shrinkage  is  determined  by  the  type  and  concentration  of  the 
electrolyte  [140,238].  It  should  be  pointed  out  that  some  particulate 
systems,  such  as  alumina  gels  and  some  clays,  exhibit  swelling  rather 
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than  syneresis  [239],  because  of  invasion  of  water  in  the  layers  causing 
expansion  [140]. 

Vysotskii  et  al.  [240-242]  have  shown  that  the  shrinkage  rate  during 
syneresis  of  silica  gel  has  a minimum  at  the  isoelectric  point  (IEP).  For 
silica  systems  the  IEP  is  at  a pH  of  1.7  [33,140].  Vysotskii  et  al.  [241] 
and  Sheinfain  and  Neimark  [243]  concluded  that  condensation  reactions 
are  responsible  for  both  gelation  and  syneresis  in  silica  gel  systems. 

Acker  [244]  and  Jones  and  Fischbach  [245]  reported  that  there  is  a 
minimum  solids  content  (4g  to  5g  of  Si02  per  100  ml  of  sol)  below 
which  no  observable  shrinkage  occurs. 

Brinker  and  Scherer  [140]  and  Scherer  [237]  suggested  that  syneresis 
is  affected  not  only  by  condensation  reactions,  but  also  by  the  tendency 
to  reduce  the  large  solid/liquid  interfacial  area  of  the  gel. 

Assuming  that  the  activation  energy  for  formation  of  bonds  is  given 
by  AEb,  then  the  rate  of  formation  of  bonds  (dn/dt)  is  given  by  [140] 


where  k'  is  a constant.  Assuming  that  the  activation  energy  for  viscous 
flow  is  AEf,  then  the  viscosity  of  the  solid  phase  is  given  by  [140] 


dn  , f-AEb" 
, - k exp  R T 


(1.22) 


r|  = k"  exp  , 


(1.23) 
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where  k"  is  a constant.  Thus,  the  shrinkage  rate  dV/dt  can  be  written  as 
[140] 


dV  k dn  f-(AEb  + AEf)^ 

dt  “ dt  “ k exp  [ R T 


(1.24) 


where  k is  a constant  (k'k").  Equation  (1.24)  indicates  that  the 
activation  energy  for  shrinkage  is  the  sum  of  the  activation  energies  for 
bonding  and  flow  [140]. 

Brinker  and  Scherer  [140]  suggest  that  the  shrinkage  rate  decreases 
(as  the  process  continues)  because  of  stiffening  of  the  network,  rather 
than  depletion  of  reactive  groups.  The  absence  of  syneresis  in  gels  with 
very  small  solids  contents  [244,245]  may  indicate  that  the  low  modulus 
of  the  network  is  inadequate  to  squeeze  liquid  out  of  the  pores  at  an 
appreciable  rate  [140]. 

According  to  Brinker  and  Scherer  [140],  coarsening  (or  ripening)  is 
a process  of  dissolution  and  reprecipitation  driven  by  differences  in 
solubility  (s)  between  surfaces  with  different  radii  of  curvature  (r). 
The  solubility  is  given  by  [140] 


s = s0  exp 


(2  Ysl  V m 

R T r 


\ 


/ 


(1.25) 


where  sQ  is  the  solubility  of  a flat  plate; 

Ysl  is  the  solid/liquid  interfacial  energy; 
Vm  is  the  molar  volume  of  the  solid. 
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The  amount  of  shrinkage  that  occurs  during  drying  depends  on  the 
stiffness  of  the  network  [140,246].  If  the  gel  is  aged  under  conditions  of 
high  solubility,  the  network  may  resist  compression  by  capillary  forces 
during  drying  [140], 

Examples  of  phase  transformations  that  can  occur  during  aging 
include  microsyneresis  (in  which  the  solid  phase  separates  from  the 
liquid  on  a local  scale),  and  segregation  of  the  liquid  into  two  or  more 
phases  [140,209].  Aging  may  also  lead  to  crystallization  [247,248]. 

According  to  Brinker  and  Scherer  [140],  the  stiffer  and  stronger  the 
gel  network  becomes  the  better  it  can  withstand  the  capillary  pressure  so 
aged  gels  shrink  and  crack  less  during  drying.  The  structure  of  the  gels 
also  influences  the  sintering  behavior,  since  the  densification  process  is 
driven  by  the  interfacial  energy,  which  depends  on  the  structure  of  the 
gel. 

1.3.9  Drying 

The  drying  behavior  of  porous  solids  have  been  extensively  studied 
by  Sherwood  [249-251],  Moore  [252],  Keey  [253],  Whitaker  [254], 
Cooper  [255],  Mujumdar  [256,257],  Ford  [258]  and  Scherer  [259-266]. 
However,  most  of  the  data  reported  by  those  researchers  have  been  on 
powder  systems  with  relatively  large  pore  sizes.  The  studies  conducted 
by  Kawaguchi  [267]  and  Dwivedi,  [268]  have  also  been  on  large  pore 
size  gels. 
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Classical  theory  of  drying  divides  the  process  into  three  stages  [258]: 

Stage  1 : The  constant-rate  period,  during  which  the  evaporation  rate 
is  independent  of  the  moisture  content. 

Stage  2:  The  first  falling-rate  period,  during  which  the  evaporation 
rate  is  usually  a linear  function  of  the  moisture  content 
[269-271]. 

Stage  3:  The  second  falling-rate  period,  during  which  the 

evaporation  rate  is  a non-linear  function  of  the  moisture 
content. 

For  gel  systems,  during  the  Stage  1,  the  evaporation  of  liquid  is 
directly  responsible  for  the  decrease  in  volume  of  the  gel.  For  large 
pore  size  systems  this  first  stage  of  drying  is  called  the  constant  rate 
period  (CRP)  [140],  because  the  evaporation  rate  per  unit  area  of  the 
drying  surface  is  independent  of  time  [252,272].  The  rate  of 
evaporation  Re  is  given  by  [28,140,258] 

Re  = kE  (pv  - PA) , (1.26) 

where  kE  is  a constant; 

py  is  the  vapor  pressure  at  the  evaporating  surface; 

PA  is  the  ambient  vapor  pressure. 

Dwivedi  [268]  shows  that  the  rate  of  water  loss  was  constant  for  his 
gels  during  Stage  1,  as  seen  in  Figure  1.5.  Suzuki  and  Maeda  [273] 
show  that  the  evaporation  rate  is  constant  even  when  clay  patches  form 
on  the  surface  of  the  body. 


Rate  of  water  loss  (g  cm min 
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Figure  1 .5  Plot  of  the  rate  of  water  loss  from  alumina 
gel  versus  water  content  of  the  gel  [268]. 
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Brinker  and  Scherer  [140]  show  that  a tensile  stress  (ax)  is  created 
because  of  a shrinkage  rate  difference  between  the  inside  and  outside  of 
a drying  body.  The  tensile  stress  is  given  by  [140] 


__  L til  Re 
Gx  3D  ’ 


(1.27) 


where  L is  the  thickness  of  the  drying  body; 

T|l  is  the  viscosity  of  the  liquid; 

Re  is  the  evaporation  rate; 

D is  the  permeability  of  the  network. 

Quantitative  analysis  of  the  drying  kinetics  of  acid  catalyzed  alkoxide 
gels  by  Wilson  [274]  and  Wilson  and  Hench  [275,276]  show  that  for  gels 
with  average  pore  radius  < 20  nm  the  rate  of  evaporation  in  Stage  1 is 
not  constant,  but  decreases  substantially  with  time  as  shown  in  Figure 
1.6  [274].  As  shown  in  Figure  1.7  by  Wilson  and  Hench  [276],  there  is 
a pronounced  shrinkage  of  the  alkoxide  derived  silica  gel  during  Stage 
1.  During  Stage  2 the  dimensions  of  the  silica  gel  remain  practically 
unchanged.  The  relationship  between  relative  weight  (W/W0)  and 
relative  volume  (V/V 0)  of  the  silica  gel  during  drying  is  shown  in 
Figure  1.8  [274].  According  to  Hench  and  co-workers  [28,274-280]  the 
decreasing  rate  of  evaporation  throughout  Stage  1,  results  from  the 
effective  pore  radius  (the  pore  radius  minus  a layer  of  bond  water) 
decreasing  from  3.2  nm  to  1.3  nm. 

The  Stage  2 of  drying  starts  when  the  "critical  point"  is  reached. 
During  this  stage  shrinkage  stops  and  the  liquid  flows  from  the  pores  to 
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Figure  1.6  Time  dependence  of  (a)  temperature  and 
weight;  (b)  absolute  loss  rate;  and  (c)  loss  rate  per  unit 
area  for  silica  gel  [274]. 
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Figure  1 .7  Time  dependence  of  (a)  linear  shrinkage,  and 
(b)  volumetric  shrinkage  of  silica  gel  [2741. 
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Figure  1.8  Weight  reduction  (W/W0)  as  a function  of 
volumetric  shrinkage  (V/V 0)  for  silica  gel  [274], 
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the  surface  where  evaporation  occurs.  Whitaker  [281]  describes  the 
flow  of  liquid  as  being  driven  by  the  gradient  in  capillary  stress. 

According  to  Wilson  [274]  and  Wilson  and  Hench  [275,276],  shortly 
after  entering  Stage  2 the  gel  turns  opaque,  starting  at  the  edges  and 
progressing  linearly  toward  the  center.  Possible  causes  of  this 
phenomena  includes  phase  separation  of  the  pore  liquid  or  exsolution  of 
gas  from  the  liquid.  Shaw  [282,283]  suggests  that  the  opaque  stage  is 
caused  by  light  scattering  from  isolated  pores  (or  groups  of  pores)  in 
the  emptying  process. 

According  to  Hench  and  West  [28]  the  transition  to  Stage  3 is  the 
hardest  to  identify,  and  its  "start"  is  probably  best  defined  as  the  "end" 
of  the  opaque  stage.  Wilson  and  Hench  [276]  show  that  Stage  3 is  not 
affected  by  local  changes  in  temperature,  ambient  vapor  pressure,  flow 
rate,  etc. 

Hench  and  West  [28,  p.75]  point  out  that  "successful  drying  of  large 
gel  monoliths  require  control  of  the  drying  rate  through  the  opaque 
stage  and  elimination  of  processing  defects  during  mixing,  casting,  and 
gelation." 


CHAPTER  2 

DENSIFICATION  OF  GEL-SILICA 
2.1  The  Sintering  Process 


2.1.1  Introduction 

Based  on  the  pioneering  work  of  Rhines  [284]  densification  (used  as 
a synonym  of  sintering  in  this  work)  can  be  defined  as  a geometric 
change  in  which  a mass  of  separate  particles  becomes  a coherent  body, 
motivated  by  a decrease  in  surface  free  energy.  Even  though  this  seems 
a simple  and  clear  definition  (and  it  is),  Rhines  [285,  p.480]  points  out 
that  despite  of  the  fact  that  "a  number  of  theories  of  sintering  have  been 
offered,  no  single  theory  accounts  for  the  entire  range  of 
[experimental]  observations."  According  to  Johnson  [286,  p.137],  "the 
lack  of  understanding  of  sintering  is  manifested  in  the  diversity  of 
opinions  as  to  the  mechanisms  of  the  process,  and  the  inability  to  predict 
the  sintering  behavior  of  materials."  Exner  and  Petzow  [287,  p.  1 17]  go 
on  further  stating  that  "a  general  shrinkage  equation  presumably  will 
never  be  derived  due  to  the  high  complexity  and  the  large  variations  of 
the  sintering  systems." 


Despite  the  complexity  of  the  sintering  process,  the  driving  force  for 
the  study  of  sintering  is  enormous,  since  most  ceramic  materials  and  a 
large  number  of  metallic  systems  need  to  be  sintered  in  some  stage  of 
their  processing.  The  major  applications  of  the  study  of  sintering  are 

1)  Better  sintering  products. 

2)  Improving  sintering  processing. 

3)  Better  understanding  of  the  process. 

The  driving  force  for  densification  is  a decrease  in  surface  free 
energy.  Besides  densification  itself,  this  decrease  in  free  energy  also 
leads  to  other  effects  such  as  [284] 

1)  Rounding  and  smoothing  (Figure  2.1(a)). 

2)  Isolation  of  pores  (Figure  2.1(b)). 

3)  Segregation  of  pores  (Figure  2.1(c)). 

None  of  the  effects  above  can  by  themselves  increase  the  density  of 
the  material.  Densification  only  occurs  when  the  bulk  density  of  the 
material  decreases  as  a result  of  the  decrease  in  volume  fraction  of 
pores  (Vv)  [288]. 

The  pioneering  attempts  to  quantify  the  geometric  changes  during 
sintering  dates  back  to  the  1920s  with  Tamman  [289]  and  Hedvall  and 
Helin's  [290]  works.  They  were  followed  shortly  after  by  Balshin  [291], 
Huttig  [292-298],  Frenkel  [299],  Pines  [300],  Mackenzie  and 
Shuttleworth  [301]  and  Kuczynski  [302-304]. 
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Figure  2.1  (a)  Rounding  and  smoothing;  (b)  isolation  of  pores; 
and  (c)  segregation  of  pores. 
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2.1.2  Mechanisms  of  Sintering 

Some  of  the  primary  proposed  mechanisms  of  sintering  are 

1)  Vapor  transport  or  evaporation-condensation.  The  transport  of 
mass  occurs  due  to  a gradient  of  pressure,  which  tends  to  decrease 
the  total  surface  energy.  This  mechanism  normally  does  not  lead 
to  densification,  but  only  to  local  surface  rounding  [284]. 

2)  Liquid  transport  or  dissolution-precipitation.  This  mechanism  is 
very  similar  to  vapor  transport  and,  like  that,  liquid  transport 
leads  only  to  local  surface  rounding. 

3)  Surface  diffusion.  Like  the  previous  mechanism,  usually  this 
mechanism  leads  only  to  surface  rounding.  In  this  case  the 
transport  of  mass  is  done  via  the  motion  of  atoms  on  the  surface 
of  the  material  [284]. 

4)  Volume  diffusion.  In  this  mechanism  the  transport  of  mass  is 
done  via  motion  of  vacancies  in  the  volume  of  the  material.  The 
main  effect  of  this  mechanism  is  some  very  localized  change  in 
density  but  it  does  not  promote  overall  densification  [284]. 

5)  Bulk  movement.  This  is  the  main  mechanism  (in  fact  group  of 
mechanisms)  that  usually  leads  to  the  overall  densification.  This 
mechanism  requires  growing  of  the  neck  between  particles,  which 
is  caused  by  a decrease  in  the  surface  free  energy.  Depending 

on  the  type  of  material,  two  types  of  bulk  movement  can  occur: 
a)  Crystalline  materials:  the  transport  of  mass  is  done  by  means 
of  atomic  diffusion  in  the  bulk. 
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b)  Amorphous  materials:  the  mechanism  of  transport  of  mass  is 
by  viscous  flow  [3021. 

6)  Reaction  sintering.  It  is  a mechanism  in  which  the  densification 
occurs  by  means  of  chemical  reaction. 

2.1.3  The  Rhines-DeHoff  Model 

According  to  Rhines  [284,  p.35],  "despite  the  apparent  chaos,  a 
description  of  some  features  of  the  geometry  of  sintering  is  made  both 
possible  and  comparatively  simple  by  the  use  of  the  concepts  of 
topology." 

Chapter  3 presents  a review  on  topology  and  its  important  concepts. 
Here  it  is  enough  to  state  that  one  of  the  most  important  concepts  of 
topology  is  the  genus  (G),  or  connectivity,  which  can  be  described  as 
[284,305] 

G = B - N + 1,  (2.1) 

where  B is  the  number  of  interparticle  contacts  or  branches; 

N is  the  number  of  particles  or  nodes. 

Based  on  the  connectivity,  the  sintering  process  can  be  divided  into 
three  stages,  as  shown  in  Figure  2.2.  The  characteristics  of  the  three 
stages  of  sintering  are  [284] 

First  stage:  - the  genus  is  constant  at  (B  - N + 1); 

- weld  necks  grow; 
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Figure  2.2  The  three  stages  of  sintering. 
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- the  surface  area  is  maximal. 

Second  stage:  - the  genus  diminishes  from  (B  - N + 1)  towards  zero; 

- the  pores  become  isolated. 

Third  stage:  - the  genus  is  constant  at  zero; 

- the  pores  disappear. 

As  will  be  pointed  out  in  Chapter  3,  the  main  limitation  of  applying 
topological  concepts  to  sintering  studies  is  the  difficulty  in 
experimentally  obtaining  topological  parameters.  Because  of  this, 
relatively  little  work  has  been  done  is  applying  the  Rhines-DeHoff 
model  to  sintering  studies.  The  most  important  work  in  this  area  was 
carried  out  by  the  so-called  "Rhines  School"  [306],  including  the  work 
done  at  the  University  of  Florida  by  Buteau  [307],  Steele  [308], 
Aigeltinger  [309]  and  Pinheiro  [310]. 

2.1.4  Frenkel's  Model 

Frenkel  [299]  proposed  a sintering  model  based  on  the  occurrence  of 
viscous  flow.  This  model  is  based  on  the  observation  that  the  different 
mechanisms  of  sintering  can  be  distinguished  by  the  particular 
relationship  of  the  radius  of  the  neck  between  two  particles  (x)  and  the 
sintering  time  (t),  as  follows  [302]: 

Viscous  flow  x2  ~ t 

Vapor  transport  x3  ~ t 

Volume  diffusion  x5  ~ t 

Surface  diffusion  x7  ~ t 


The  basic  relationship  used  by  Frenkel  is  [299] 


6r  = :jf . (2.2) 

where  £r  is  the  strain  rate; 

f is  the  acting  stress  (caused  by  surface  tension); 

T|  is  the  viscosity  of  the  material. 

Frenkel  shows  that  for  a viscous  flow  regime  the  following  equation 
will  be  obeyed  [299,302]: 


'_x_N|  2 _ 3_  t y 

, a J “ 2 ari  ’ 


(2.3) 


where  a is  the  original  radius  of  the  bonded  particles; 
y is  the  surface  energy  of  the  material. 

Kuczynski  found  very  good  agreement  of  this  model  for  the 
sintering  of  spherical  glass  particles  to  glass  plates  [302]. 

The  major  limitations  of  Frenkel's  model  are  the  following 
assumptions  of  the  model: 

1)  The  radius  "a"  is  kept  constant  during  the  entire  process.  This 
implies  that  Eq.(2.3)  is  not  applicable  beyond  the  early  stages  of 
sintering.  It  is  usually  applicable  for  values  of  linear  shrinkage 
smaller  than  10%. 
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2)  The  viscosity  of  the  solid  (r|)  is  given  by 


k T 


^ D GJ1/3  ’ 


(2.4) 


where  k is  the  Boltzmann  constant; 

T is  the  absolute  temperature; 

D is  the  coefficient  of  self-diffusion; 

G5  is  the  atomic  volume. 

Nabarro  [311]  showed  that  Frenkel's  calculation  of  viscosity  is  not 
correct  for  crystalline  materials. 

3)  The  material  is  homogeneous  and  has  stable  values  of  "a"  and  7]. 

2.1.5  Mackenzie-Shuttleworth's  Model 

The  Mackenzie-Shuttleworth’s  model  (or  simply  the  MS  model) 
[301]  is  also  a viscous  flow  model.  It  is  based  on  a geometric  model 
(Figure  2.3)  proposed  by  Frohlich  and  Sack  [312].  The  basic  equation 
assumed  by  the  Mackenzie-Shuttleworth's  model  is  the  following  [301]: 


where  p is  the  density  of  the  compact  relative  to  the  true  density  of  the 
material; 

ri  and  T2  are  defined  in  Figure  2.3. 


(2.5) 


Figure  2.3  Geometric  model  used  in  the 
Mackenzie-Shuttle  worth's  sintering  model  [301]. 
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Depending  on  the  rheological  behavior  of  the  material  one  can  obtain 
different  equations  correlating  the  progress  of  the  sintering  with  time. 
Assuming  for  instance,  that  the  material  presents  Newtonian  viscosity, 
one  obtains  [301] 


(t  - to) 


2 ( 3 V/3  ( 1 , 

o In 


47C 

vv 


(l+x3) 

(1+x)3 


V 3 tan-1 


2x-r 

V5j 


, (2.6) 


where  n is  the  number  of  pores  per  unit  volume; 
x is  equal  to  (1  - p)/p; 
to  is  the  initial  time. 

Assuming  Bingham  viscosity,  the  equation  obtained  by  the  MS  model 
is  [301] 
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dt 
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where 


= V2 


( 3 M/3  Tc 


4 K 


2 yni/3  ’ 


(2.8) 


and  Tc  is  the  critical  shear  stress  for  Bingham  flow. 

Sacks  and  Tseng  [313]  found  good  experimental  agreement  of 
sintering  of  colloidal  silica  glass  spheres  (Stober  spheres)  with  the  MS 
model. 
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Important  limitations  of  the  MS  model  are  the  following 
assumptions: 

1)  The  pores  are  isolated.  This  means  that  this  model  applies  better 
to  the  last  stages  of  sintering,  when  the  body  contains  only 
isolated  closed  pores. 

2)  The  material  is  incompressible. 

3)  The  material  has  homogeneous  and  stable  mechanical  properties. 
2.1.6  Scherer's  Model 

Scherer's  model  [314-318]  is  also  a viscous  flow  model.  It  is  a 
geometric  model  which  assumes  a cubic  array  formed  by  intersecting 
cylinders  to  represent  the  structure  of  the  solid,  as  shown  in  Figure  2.4 
[318]. 

According  to  Scherer's  model  the  rate  of  densification  is  given  by 
[140] 

y n1/3  t 

JL-^~  = f*  (y)  - fs  (yo) , (2.9) 


where 


fs  (y)  = - 


hn 


{a2-ay+y2}  , 

n4r + V3tairl 


fly-a^ 

aV3 


(2.10) 


and 
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Figure  2.4  Geometric  model  used  in  the  Scherer's 
sintering  model  [318]. 
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(2.11) 


and 


a = (8  V2)1/3  . 


(2.12) 


Sacks  and  Tseng's  work  [313]  show  that  the  Scherer's  model  can  be 
used  to  describe  the  sintering  of  spherical  silica  particles.  The  model 
has  also  been  used  to  describe  the  sintering  of  a body  containing  two 
different  pore  sizes  [318]  and  a distribution  of  pore  sizes  [140]. 

An  important  limitation  of  Scherer's  model  is  that  the  model  assumes 
a very  simple  geometry  which  cannot  be  applied  to  the  later  stages  of 
densification  when  the  cylinders  start  touching  each  other. 

2.1.7  Densification  of  Gel-Silica 

The  first  stages  of  densification  of  gel-silica  are  termed 
"stabilization"  stages  [28]  which  can  be  divided  into  chemical  and 
thermal  stabilization.  Chemical  stabilization  involves  removing  surface 
silanols  from  the  surface  until  a critical  concentration  is  reached  so  that 
the  surface  will  not  rehydroxy  late  [28].  The  thermal  stabilization 
relates  to  the  material  being  able  to  be  used  at  a given  temperature 
without  reversible  structural  changes  [28]. 

Many  chlorine  compounds  can  react  with  surface  hydroxyl  groups  to 
form  hydrochloric  acid  [319-322],  which  then  desorbs  from  the  gel 
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body  at  a temperature  range  of  400  °C  (673  K)  to  800  °C  (1073  K)  [3]. 
Examples  of  chlorine  compounds  that  can  be  used  for  dehydroxylation 
of  the  silica  gel  surface  include  [28]  ClSi(CH3)3,  Cl2Si(CH3)2, 
Cl3Si(CH3),  silica  tetrachloride  (SiCU),  chlorine  (CI2),  and  carbon 
tetrachloride  (CCI4). 

Water  in  silica  gels  is  present  in  two  forms  [28]:  free  water  within 
the  porous  gel  structure  (i.e.,  physical  water)  and  hydroxyl  groups 
associated  with  the  gel  surface  (i.e.,  chemical  water).  The  work  done 
on  the  hydration/dehydration  characteristics  of  the  silica  gel/water 
system  by  Benesi  and  Jones  [323],  Kiselev  [324],  McDonald  [325], 
Young  [326],  Hockey  and  Pethica  [327],  and  Hair  [328]  can  be 
summarized  as  follows: 

1)  Physical  water  can  be  eliminated  and  surface  silanol  (Si-O-H) 
groups  condensed  at  temperatures  starting  at  about  170°C  (443 
K). 

2)  The  dehydration  process  is  completely  reversible  up  to  about  400 
°C  (673  K). 

3)  At  temperatures  above  400  °C  (673  K)  the  dehydration  is 
irreversible  as  a result  of  shrinkage. 

4)  Viscous  flow  occurs  above  850  °C  (1 123  K). 

With  a successful  stabilization  treatment,  it  is  possible  to  manufacture 
monolithic  dense  gel-derived  glasses  by  heating  at  temperatures  below 
the  melting  point  [30,100,107,206,329-334]. 

According  to  Hench  and  West  [28]  the  amount  of  water  in  the  gel  has 
a major  role  in  the  densification  behavior.  The  amount  of  water  affects 
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the  viscosity  [335],  which  influences  the  sintering  behavior.  As  shown 
by  Nogami  and  Moriya  [206],  a gel  prepared  in  acidic  conditions  has  a 
higher  surface  area  and  water  content  than  a gel  prepared  in  basic 
conditions,  and  starts  to  densify  about  200  °C  sooner  than  the  base 
catalyzed  gel. 

Young  [326]  found  that  the  decrease  in  surface  area  of  silica  gel  at 
high  temperatures  is  a function  of  the  time  and  temperature  of  the  heat 
treatment.  Many  other  authors  [164,336-340]  present  a similar  type  of 
result.  However,  Klein  et  al.  [331]  and  Brinker  et  al.  [341]  showed  that 
certain  systems  display  an  initial  increase  in  surface  area  until  about  400 
°C  (673  K)  and  then  the  surface  area  decreases  with  further  heating. 
The  initial  increase  in  surface  area  was  attributed  to  the  removal  of 
water  and  organics. 

According  to  Brinker  and  Scherer  [342],  there  are  at  least  four 
mechanisms  that  contribute  for  the  shrinkage  and  densification  of  gels. 
Those  four  mechanisms,  as  described  for  a borosilicate-gel  system,  are 

1)  Capillary  contraction,  due  to  increase  in  surface  energy  associated 
with  elimination  of  silanol  (Si-OH)  groups,  occurs  mostly  at 
temperatures  below  150  °C  (423  K)  and  the  shrinkage  associated 
with  this  region  is  about  3%  of  the  total  shrinkage  observed  in  the 
experiment. 

2)  Condensation-polymerization,  involves  formation  of  siloxane 
(Si-O-Si)  groups  and  related  loss  of  water  and  organic  residues. 

3)  Structural  relaxation,  occurs  via  diffusive  motion  of  atoms. 

This  mechanism  and  the  condensation-polymerization  lead  to 
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skeletal  (structural)  densification.  These  mechanisms  occur 
mostly  between  150  °C  (423  K)  and  525  °C  (798  K)  and  the 
shrinkage  associated  with  skeletal  densification  is  about  33%  of 
the  total  shrinkage. 

4)  Viscous  sintering,  associated  with  mass  flow  driven  by  decrease  in 
surface  energy,  occurs  above  525  °C  (798  K)  and  is  responsible 
for  most  of  the  shrinkage  (about  64%). 

For  systems  other  than  the  borosilicate-gel  studied  by  Brinker  et  al. 
[342],  the  temperatures  and  shrinkages  described  above  are  likely  to  be 
different. 


2.2  Sol-Gel  Processing 
2.2.1  Chemical  Processing  of  Pure  Silica 

A schematic  of  the  sol-gel  processing  used  in  this  work  is  shown  in 
Figure  2.5.  Based  on  the  theory  reviewed  in  Chapter  1,  the  silica-gel 
process  can  be  divided  into  seven  stages:  Mixing,  Casting,  Gelation, 
Aging,  Drying,  Stabilization  and  Densification. 

During  mixing  the  chemical  components  are  put  in  an  intimate 
molecular  contact  and  hydrolysis  takes  place: 


Si(OCH3)4  + 4H20  — > Si(OH)4  + 4CH3OH 


(2.13) 
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Figure  2.5  Schematic  of  the  sol-gel  processing  of 
gel-silica  monoliths. 
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The  gelation  stage  is  characterized  by  the  occurrence  of  condensation 
reactions  up  to  a point  when  a gel  is  formed.  The  gel  formation  is 
accompanied  by  a characteristic  increase  in  viscosity.  The  main 
condensation  reaction  can  be  expressed  as  follows: 

Si(OH)4  ->  Si02  (gel)  + 2H20  (2.14) 

Further  development  of  condensation  processes  gives  rise  to 
syneresis  or  network  shrinkage,  characterizing  the  aging  stage.  During 
the  drying  stage  water  is  eliminated  from  the  gels  driven  by  thermal 
energy  and  capillary  forces.  Shrinkage  during  drying  is  caused  mainly 
by  capillary  forces  but  chemical  reactions,  osmotic  forces  and  other 
factors  may  play  a role  as  well  [140]. 

During  densification,  OH  groups  are  eliminated  from  the  structure 
driven  by  increased  thermal  energy  and  by  chemical  reactions  such  as 
chlorination  of  the  surface.  The  shrinkage  and  increase  in  density  that 
occur  during  sintering  are  driven  by  the  decrease  in  surface  area  of 
pores. 

2.2.2  Gelation  and  Aging 

The  gelation  of  the  samples  prepared  in  this  work  was  done  by 
leaving  the  containers  as  cast  with  the  solution  at  room  temperature 
until  the  gel  formation,  which  typically  occurred  in  two  to  three  days 
after  casting. 
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After  gelation  the  samples  were  put  in  an  oven  in  the  same  closed 
polystyrene  (PS)  containers  used  for  casting.  One  of  the  aging 
temperature  profiles  used  is  shown  in  Figure  2.6. 

2.2.3  Drying 

After  aging  the  samples  were  taken  out  of  the  casting  PS  containers 
and  transferred  to  Teflon®  containers  which  were  filled  with  DI  (de- 
ionized) water.  The  Teflon®  containers  were  closed  and  a small 
opening  (capillary)  was  made  in  the  containers. 

Drying  temperature  profiles  like  the  one  shown  in  Figure  2.7  were 
used  for  drying  of  the  samples. 

Using  the  above  procedures,  four  sets  of  cylindrical  samples  with 
different  sizes  were  prepared.  These  samples  are  termed  "samples  type 
A",  and  the  dimensions  of  the  samples  are 

Al:  diameter  12  mm;  length  =12  mm. 

A2:  diameter  5 mm;  length  = 8 mm. 

A3:  diameter  8.5  mm;  length  =13  mm. 

A4:  diameter  5 mm;  length  = 28  mm. 
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Figure  2.6  Schematic  of  the  aging  temperature  profile. 
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Figure  2.7  Schematic  of  the  drying  temperature  profile. 
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2.3  Large  Pore  Size  Structures 

2.3.1  HNCh  Based  Gels 

In  this  experiment  HNO3  was  used  as  catalyzer,  according  to  the 
process  described  previously  (Figure  2.5).  Using  a special  drying 
chamber  like  the  one  shown  in  Figure  2.8  [274]  it  is  possible  to  obtain 
samples  of  larger  pore  size  and  larger  initial  (dried)  volume  fraction 
than  by  using  the  Teflon®  containers.  The  samples  obtained  using  the 
drying  chamber  shown  in  Figure  2.8  are  called  "samples  type  B".  The 
internal  volume  of  the  drying  container  and  the  total  external  surface 
area  of  the  gels  inside  the  drying  chamber  seem  to  affect  the  drying 
process.  It  is  possible  that  this  effect  is  associated  with  changes  in  the 
equilibrium  vapor  pressure  of  the  chemical  components  of  the  gel  inside 
the  drying  container. 

2.3.2  HF  Based  Gels 

Replacing  part  of  HNO3  as  catalyst  for  fluoride  acid  (HF),  it  is 
possible  to  change  the  average  pore  size  and  volume  fraction  of  pores  in 
the  initial  (dried)  samples. 

In  the  processing  used  in  this  work,  6 ml  of  HNO3  was  added  to  485 
ml  of  DI  water.  After  stirring,  the  solution  was  transferred  to  a plastic 
(Teflon®)  container  and  20  ml  of  HF  was  added.  After  mixing,  250  ml 


65 


K Teflon® 
septum 

Aluminum 
end  plate 


Figure  2.8  Drying  chamber  [274]. 
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of  TMOS  was  added  to  the  solution,  and  after  mixing  for  few  minutes 
the  solution  was  poured  in  polystyrene  vials. 

After  casting  the  solution  was  left  gelling  for  18  hours,  and  the  aging 
step  was  carried  out  in  a oven  at  60  °C  (333  K)  for  two  days.  After 
aging,  the  samples  were  dried  in  Teflon®  containers  filled  with  DI 
water  and  containing  a small  capillary  for  vapor  elimination.  The 
drying  schedule  used  was  similar  to  that  shown  in  Figure  2.7.  The  HF 
based  silica  gels  are  termed  "samples  type  C". 

2.4  Densification  Processing 
2.4.1  Partial  Densification  in  Dry-Air  Atmosphere 

The  densification  experiments  were  carried  out  in  tube  furnaces  with 
atmosphere  control.  Dried  samples  were  put  in  fused  quartz  crucibles 
and  introduced  in  the  furnaces. 

Samples  were  densified  using  the  temperature  profiles  shown  in 
Figure  2.9.  During  the  entire  process  a flow  of  1.0  liter/min.  of  dry-air 
(<  3 ppm  of  H2O)  was  used.  A small  pressure  of  12  mm  of  Hg  inside 
the  furnace  was  kept  by  using  wash  bottles  filled  with  Di-water  at  the 
end  of  the  furnace.  Figure  2.10  shows  a set  of  silica-gel  samples 
densified  in  a dry-air  atmosphere. 
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Figure  2.9  Schematic  of  the  densification  temperature  profiles. 


Figure  2.10  Examples  of  pure  silica  gel  monoliths  densified  in  the 
range  of  180  °C  to  1000  °C  in  a dry-air  atmosphere. 
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2.4.2  Stabilization  and  Densification  Processing  Using  Chlorination 

Chlorine  on  the  surface  of  silica-gels  reacts  with  OH  groups 
promoting  dehydroxylation  before  pore  closing.  A typical  chlorination 
temperature  profile  is  shown  in  Figure  2.11.  The  step  of  heating  to  an 
intermediate  temperature  followed  by  a cooling  to  about  400  °C  (673  K) 
is  intended  to  promote  some  dehydroxylation  and  strength  at  the  same 
time  that  pore  closing  is  avoided.  After  the  pre-stabilization  in  a dry- 
air  atmosphere  the  samples  were  held  in  a chlorinated  atmosphere  for  a 
period  of  48  hours.  After  the  chlorination  treatment  the  silica-gels 
were  taken  to  1150  °C  (1423  K)  under  a helium  atmosphere.  Figure 
2.12  shows  samples  prepared  in  dry-air  atmosphere  (500  °C  (773  K)  to 
1000  °C  (1273  K))  and  a full  dense  sample  (1150  °C,  1423  K)  prepared 
in  a chlorine  atmosphere. 

2.5  Gas  Adsorption  Analysis 


2.5.1  BET  Theory 

The  adsorption  of  gases  on  a solid  surface  can  be  divided  into  two 
broad  categories  of  chemical  or  physical  adsorption,  depending  on  the 
strength  of  the  interaction  between  the  adsorbate  (gas  or  vapor 
molecule)  and  adsorbent  (surface  of  the  solid).  The  chemical  adsorption 
(chemisorption)  involves  large  interaction  potentials  associated  with 
chemical  bonding  which  makes  the  process  irreversible.  Physical 
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Figure  2.11  Temperature  profile  for  a chlorination  experiment. 
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Figure  2.12  Silica  gel  monoliths  prepared  in  a dry-air 
atmosphere  (500  °C  to  1000  °C)  and  in  a chlorine 
atmosphere  (1150  °C). 
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adsorption  (or  reversible  adsorption)  has  the  following  characteristics 
[343]: 

1)  Low  heats  of  adsorption  and  no  disruptive  structural  changes  on 
the  solid  surface  during  the  adsorption  process. 

2)  It  may  lead  to  surface  coverage  by  more  than  one  layer  of 
adsorbate,  enabling  the  measurement  of  pore  volume. 

3)  Physical  adsorption  does  not  occur  at  elevated  temperatures. 

4)  Since  no  (or  extremely  small)  activation  energy  is  involved  the 
adsorption  process  can  be  achieved  rapidly. 

5)  The  adsorption  process  is  reversible,  enabling  the  study  of  the 
desorption  process  as  well. 

6)  Physically  adsorbed  molecules  are  free  to  cover  the  entire  surface 
and  can  be  used  to  calculate  surface  areas  (rather  than  number  of 
sites). 

The  types  of  interactions  that  lead  to  physical  adsorption  are  the 
following  [343]: 

1)  Dispersion  (or  London  [344])  forces; 

2)  Ion-dipole  interactions; 

3)  Ion-induced  dipole  interactions; 

4)  Dipole-dipole  interactions; 

5)  Quadrapole  interactions. 

When  analyzing  the  pore  structure  through  gas  adsorption  it  is 
convenient  to  classify  the  pores  according  to  their  size,  as  in  the 
Dubinin's  classes  of  pores  [343]: 

1)  Micropores,  for  pore  radius  < 15  A; 
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2)  Mesopores,  pore  radius  in  the  range  15  A to  500  A; 

3)  Macropores,  pore  radius  in  the  range  500  A to  1000  A; 

Nitrogen  gas  adsorption  analysis  is  usually  not  employed  for 

structures  with  pore  radii  larger  than  1000  A (100  nm). 

According  to  Brunauer  et  al.  [345],  the  adsorption  of  gas  on  solid 
surfaces  at  a fixed  temperature  follows  one  of  five  possible  types  of 
"adsorption  isotherms"  showed  in  Figure  2.13.  A brief  description  of 
the  five  types  of  isotherms  is  as  follows  [343]: 

Type  I:  This  type  of  isotherm  occurs  when  adsorption  is  limited 

to,  at  most,  only  a few  molecular  layers.  In  the  case  of 
chemical  driven  adsorption  (chemisorption),  this  type  of 
isotherm  indicates  that  all  the  surface  sites  are  occupied 
for  values  of  P/P0  tending  to  1.  In  the  case  of  physical 
adsorption,  type  I isotherms  are  encountered  for 
micropore  powder  systems  where  the  pore  size  does  not 
exceed  a few  adsorbate  molecular  diameters. 

Type  II:  This  type  of  isotherm  represents  the  most  common  type 

of  adsorption  in  nonporous  powders,  or  powders  with 
pore  diameters  larger  than  micropores. 

Type  III:  This  type  of  isotherm  occurs  when  the  heat  of  adsorption 
is  less  than  the  adsorbate  heat  of  liquefaction. 

Type  IV:  This  type  of  isotherm  occurs  on  porous  adsorbents 

possessing  pore  radii  in  the  range  of  15  A to  1000  A. 
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Figure  2.13  The  five  types  of  adsorption  isotherms  as 
classified  by  Brunauer  et  al.  [345].  W is  the  weight  of 
gas  adsorbed  and  P/PQ  is  the  relative  pressure. 
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Type  V:  This  type  of  isotherm  occurs  when  the  interactions 

between  adsorbate-adsorbent  are  small  or  when  the  pore 
radii  are  in  the  range  15  A to  1000  A. 

Langmuir  [346]  was  able  to  describe  the  type  I isotherm  using  a 
kinetic  approach  and  the  assumption  that  adsorption  was  limited  to  a 
monolayer.  Brunauer,  Emmett  and  Teller  [347]  developed  later  a more 
general  theory  (the  "BET  theory")  that  can  be  used  to  describe  all  the 
five  types  of  isotherms.  The  BET  theory  assumes  that  the  uppermost 
molecules  in  adsorbed  stacks  are  in  dynamic  equilibrium  with  the  gas  or 
vapor.  The  BET  theory  assumes  that  the  equilibrium  between  the  vapor 
and  the  adsorbate  in  the  first  layer  are  described  by  the  Langmuir 
theory.  The  final  form  of  the  BET  equation  is  given  by  [343] 

1 1 C - 1 

W[(P0/P)-1]  " WmC  + WmC 


Po 


(2.15) 


where  W is  the  weight  of  vapor  adsorbed; 

Wm  is  the  weight  of  vapor  in  a completed  monolayer; 

P is  the  adsorbate  equilibrium  pressure; 

P0  is  the  adsorbate  saturated  equilibrium  vapor  pressure; 
C is  the  BET  constant,  given  by  [343] 
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where  Ai  is  the  condensation  coefficient,  or  the  probability  of  a 
molecule  being  adsorbed  upon  collision  with  the 
surface; 

A2  is  the  condensation  coefficient  for  the  second 
monolayer; 

vi  is  the  vibrational  frequency  of  the  adsorbate  normal 
to  the  surface  when  adsorbed; 

V2  is  the  vibrational  frequency  for  the  second  monolayer; 
Ei  is  the  energy  of  adsorption; 

L is  the  heat  of  liquefaction. 

When  the  weight  adsorbed  (W)  is  plotted  versus  the  relative  pressure 
(P/P0)  for  both  adsorption  and  desorption  processes  a hysteresis  loop  is 
commonly  formed,  de  Boer  [348]  has  identified  five  types  of  hysteresis 
loops  as  shown  in  Figure  2.14.  de  Boer  [348]  correlated  each  type  of 
hysteresis  with  pore  shapes,  as  follows: 

Type  A:  due  principally  to  cylindrical  pores  open  at  both  ends. 

Type  B:  associated  with  slit-shaped  pores  or  the  space  between 
parallel  plates. 

Type  C:  produced  by  a mixture  of  tapered  or  wedge-shaped  pores 
with  open  ends. 

Type  D:  produced  by  tapered  or  wedge-shaped  pores  but  with 
narrow  necks  at  one  of  both  open  ends. 

Type  E:  results  from  McBrain's  [349]  "bottle  neck"  pores. 
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Type  A 


Figure  2.14  de  Boer's  five  types  of  hysteresis  loops  [343]. 
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Despite  the  qualitative  correlation  between  type  of  hysteresis  and 
shape  of  pores,  Broekhoff  and  Linsen  [350]  point  out  that  there  is  no 
satisfactory  quantitative  theory  of  hysteresis. 

Gas  (especially  nitrogen)  adsorption  techniques  and  the  BET  theory 
have  been  extensively  applied  to  the  study  of  the  structure  of  silica  gels 
[34,351-356].  The  applicability  of  BET  theory  to  the  study  of  the 
structure  of  silica-gels  has  been  well  established. 

2.5.2  Experimental  Results 

The  gas  adsorption  analysis  was  done  using  an  automatic  gas- 
absorption  instrument,  the  Autosorb-6  (Quantachrome,  Inc.).  Based  on 
the  average  pore  radius  of  the  dried  silica  gel  monoliths  obtained  from 
the  gas  adsorption  analysis,  samples  type  A are  also  referred  to  as  "12  A 
average  pore  radius  samples";  samples  type  B (dried  in  the  large  drying 
chamber)  are  also  termed  "32  A average  pore  radius  samples";  and 
samples  type  C (HF  based)  are  also  termed  "81  A average  pore  radius 
samples". 

Figure  2.15  shows  the  BET  plot  for  a dried  gel-silica  monolith 
sample  of  12  A average  pore  radius  (sample  type  A).  The  data  shown 
in  Figure  2.15  is  used  in  Eq.(2.15)  in  order  to  calculate  the  constant 
"C".  Figure  2.16  represents  the  volumes  of  gas  adsorbed  and  desorbed 
for  samples  of  12  A average  pore  radius  densified  in  the  range  of  180 
°C  (453  K)  to  1150  °C  (1423  K).  As  the  densification  temperature 
increases  the  volume  of  adsorbed  gas  decreases,  reflecting  the  smaller 
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Figure  2.15  BET 
12 


plot  for  a dried  silica  gel  sample  with 
A average  pore  radius. 
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Figure  2.16  Volume  of  adsorbed  gas  versus  partial  pressure 
for  a set  of  silica  gel  samples  with  12  A average  pore 
radius  densified  in  the  range  of  180  °C  to  1150  °C. 
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porosity  present  as  the  sintering  proceeds.  Figure  2.17  shows  the 
variation  of  adsorbed  gas  versus  partial  pressure  for  dried  samples  of 
12  A (type  A),  32  A (type  B),  and  81  A (type  C).  The  hysteresis  of  the 
adsorption-desorption  curves  for  the  32  A and  81  A samples  suggest 
that  the  pores  in  samples  type  B and  C have  ink-bottle  shapes  (see 
Figure  2.14). 

Typical  pore  size  distributions  of  dried  samples  of  12  A,  32  A and 
81  A average  pore  radii  are  shown  in  Figure  2.18.  The  lower  limit  cut- 
off at  about  3.5  A represents  the  lowest  possible  resolution  attainable  in 
the  Nitrogen  gas  adsorption  analysis  (a  Nitrogen  atom  cannot  penetrate 
a pore  smaller  than  its  own  size). 

As  the  densification  temperature  increases  the  volume  of  pores 
decreases,  as  shown  in  Figure  2.19  for  a set  of  gel-silica  monoliths  with 
12  A average  pore  radius.  As  the  total  volume  of  pores  decreases  the 
average  pore  radius  remains  practically  unchanged.  The  pore  size 
distribution  curves  for  sets  of  larger  average  pore  radius  are  shown  in 
Figure  2.20  and  Figure  2.21.  Again,  it  is  shown  that  while  the  total 
volume  of  pores  decreases  with  the  evolution  of  densification,  the  pore 
size  distribution  keeps  about  the  same  average  pore  radius. 

Langmuir  plots  of  dried  samples  of  12  A,  32  A and  81  A average 
pore  radii  are  shown  in  Figure  2.22. 

Figure  2.23  shows  the  average  pore  radius  as  a function  of  the 
densification  temperature  for  four  sets  of  samples  (Al,  A2,  A3  and  A4) 
of  12  A average  pore  radius.  If  can  be  seen  from  Figure  2.23  that  the 
average  pore  radius  remains  fairly  constant  as  densification  proceeds,  as 
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Figure  2.17  Volume  of  adsorbed  gas  versus  partial  pressure  for  dried 
silica  gel  samples  with  12  A,  32  A and  81  A average  pore  radii. 
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Figure  2.18  Pore  size  distribution  for  dried  silica  gel  samples 
with  12  A,  32  A and  81  A average  pore  radii. 
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Figure  2.19  Pore  size  distribution  for  silica  gel  samples  with  12  A 
average  pore  radius  densified  in  the  range  of  180  °C  to  1150  °C. 
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Figure  2.20  Pore  size  distribution  for  silica  gel  samples  with  32  A 
average  pore  radius  densified  in  the  range  of  180  °C  to  1150  °C. 
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Figure  2.21  Pore  size  distribution  for  silica  gel  samples  with  81  A 
average  pore  radius  densified  in  the  range  of  180  °C  to  1150  °C. 
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Figure  2.22  Langmuir  plots  of  dried  silica  gels  with  12  A, 
32  A and  81  A average  pore  radii. 
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Figure  2.23  Average  pore  radius  versus  densification 
temperature  for  samples  type  A. 
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suggested  by  the  pore  size  distribution  curves  shown  in  Figure  2.19. 
Figure  2.24  shows  the  variation  of  the  average  pore  radius  with 
densification  temperature  for  samples  with  initial  pore  radii  of  12  A,  32 
A and  81  A. 

2.6  Evolution  of  Metric  Properties  During  Sintering 
2.6.1  Structural  Density 

True  (or  structural)  density  of  three  series  of  samples  of  12  A,  32  A 
and  81  A average  pore  radii  were  evaluated  using  a Helium 
micropycnometer  (Quantachrome,  Inc.). 

The  variation  of  true  density  with  densification  temperature  for  sets 
of  12  A,  32  A and  81  A average  pore  radii  are  shown  respectively  in 
the  Figures  2.25,  2.26  and  2.27.  Each  point  in  those  curves  represents 
the  average  of  eight  experimental  measurements.  Figure  2.28 
summarizes  the  results  shown  in  the  previous  figures.  The  data-points 
plotted  in  Figure  2.28  represent  points  on  the  polynomial  curves 
obtained  in  Figures  2.25,  2.26  and  2.27.  From  those  figures  if  can  be 
seen  that  the  structural  density  of  silica-gel  monoliths  change  as  the 
densification  proceeds.  For  all  pore  sizes  tested  the  initial  structural 
density  of  the  dried  samples  is  around  2.1  g/cm3  to  2.15  g/cm3.  As  the 
densification  temperature  increases  the  structural  density  increases  to  a 
maximum  at  an  intermediate  temperature  in  the  range  800  °C  (1073  K) 
to  1000  °C  (1273  K).  The  true  density  decreases  again  when  the 
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Figure  2.24  Variation  of  the  average  pore  radius  as  a function 
of  densification  temperature  for  silica  gel  samples  with 
12  A,  32  A and  81  A average  pore  radii. 
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Figure  2.25  Variation  of  true  density  as  a function  of  densification 
temperature  for  samples  type  A. 
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Figure  2.26  Variation  of  true  density  as  a function  of  densification 
temperature  for  samples  type  B. 
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Figure  2.27  Variation  of  true  density  as  a function  of  densification 
temperature  for  samples  type  C. 
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Figure  2.28  Variation  of  true  density  (from  polynomial  equations)  as  a 
function  of  densification  temperature  for  silica  gel  samples  with  12  A, 
32  A and  81  A average  pore  radii. 
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processing  temperature  is  raised  at  higher  values.  The  reasons  for  the 
structural  density  variations  are  not  yet  fully  known  but  it  is  possible 
that  they  come  from  changes  in  the  defect  atomic  coordination  of  Si  in 
the  structure. 

2.6.2  Metric  Properties 

The  variation  of  volume  fraction  of  pores  (Vy)  with  densification 
temperature  for  the  four  sets  of  samples  with  12  A average  pore  radius 
is  shown  in  Figure  2.29.  It  is  seen  that  as  densification  proceeds  there  is 
a continuous  decrease  in  the  volume  fraction  of  pores.  This  decrease  is 
more  pronounced  for  temperatures  just  above  the  range  800  °C  (1273 
K)  to  900  °C  (1273  K). 

Figure  2.30  shows  the  variation  of  volume  fraction  of  pores  with 
densification  temperature  for  sets  of  samples  with  average  pore  radii  of 
12  A,  32  A and  81  A.  The  samples  with  a larger  average  pore  radii 
used  in  this  work  also  present  larger  initial  volume  fraction  of  pores. 
Comparing  the  sets  of  samples  with  12  A and  32  A pore  radii,  it  is  seen 
that  the  larger  pore  size  samples  maintain  a higher  volume  fraction  of 
pores  during  the  entire  process  of  sintering  and  the  sharp  decrease  in 
volume  fraction  occurs  at  higher  temperatures.  Similar  results  occur 
for  the  set  of  samples  of  81  A pore  radius  except  that  the  sharp  decrease 
in  volume  fraction  may  be  occurring  at  lower  temperatures  than  for  the 
set  of  samples  of  32  A pore  radius. 
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Figure  2.29  Variation  of  the  volume  fraction  of  pores  (Vy)  as  a 
function  of  densification  temperature  for  samples  type  A. 
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Figure  2.30  Variation  of  the  volume  fraction  of  pores  (Vv)  as  a 
function  of  densification  temperature  for  silica  gel  samples 
with  average  pore  radii  of  12  A,  32  A and  81  A. 
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The  variation  of  surface  area  of  pores  per  unit  volume  with 
densification  temperature  for  four  sets  of  samples  of  12  A average  pore 
radius  is  shown  in  Figure  2.31.  Similar  to  the  Vy  evolution,  it  can  be 
seen  that  Sv  decreases  with  increasing  processing  temperature  and  that 
decrease  is  sharper  beyond  the  range  800  °C  (1073  K)  to  900  °C  (1173 
K).  Figure  2.32  shows  the  variation  in  Sy  pores  with  processing 
temperature  for  samples  of  12  A,  32  A and  81  A average  pore  radii. 

When  surface  area  per  unit  volume  (Sv)  is  plotted  as  a function  of 
volume  fraction  (Vy)  of  pores  for  the  samples  Al,  A2,  A3  and  A4, 
Figures  2.33,  2.34,  2.35  and  2.36  are  obtained.  Those  figures  show  a 
linear  relationship  between  Sy  and  Vy,  representing  a similar  path  of 
densification  for  those  samples.  The  linearity  between  Sv  and  Vv  also 
indicates  a constant  pore  radius  during  densification,  which  was  shown 
before  in  Figure  2.23. 

Figure  2.37  shows  the  variation  of  surface  area  of  pores  per  unit 
volume  with  volume  fraction  of  pores  for  a set  of  samples  of  32  A 
initial  average  pore  radius.  This  figure  does  not  show  the  linear 
correspondence  between  Sy  and  Vv  as  seen  for  the  12  A pore  radius. 

Using  the  data  on  pore  volume  per  unit  mass  obtained  from  the 
nitrogen  adsorption  experiments  and  the  data  on  structural  density 
obtained  from  Helium  pycnometry  one  can  calculate  the  bulk  density  of 
the  samples  by  using  the  expression: 


Pb  - Vp  + 


1.0  gram 
Pt 


5 


(2.17) 
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Figure  2.31  Variation  of  the  surface  area  of  pores  per  unit  volume  (Sv) 
as  a function  of  densification  temperature  for  samples  type  A. 
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Figure  2.32  Variation  of  the  surface  area  of  pores  per  unit  volume  (Sy) 
as  a function  of  densification  temperature  for  silica  gel  samples 
of  12  A,  32  A and  81  A average  pore  radii. 
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Vv  pores 


Figure  2.33  Variation  of  the  surface  area  of  pores  per  unit  volume  (Sy) 
as  a function  of  volume  fraction  of  pores  (Vy)  for  samples  type  Al. 
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Vv  pores 


Figure  2.34  Variation  of  the  surface  area  of  pores  per  unit  volume  (Sy) 
as  a function  of  volume  fraction  of  pores  (Vy)  for  samples  type  A2. 
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Vv  pores 


Figure  2.35  Variation  of  the  surface  area  of  pores  per  unit  volume  (Sy) 
as  a function  of  volume  fraction  of  pores  (Vv)  for  samples  type  A3. 
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Vv  pores 


Figure  2.36  Variation  of  the  surface  area  of  pores  per  unit  volume  (Sy) 
as  a function  of  volume  fraction  of  pores  (Vy)  for  samples  type  A4. 
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Vv  pores 


Figure  2.37  Variation  of  the  surface  area  of  pores  per  unit  volume  (Sy) 
as  a function  of  volume  fraction  of  pores  (Vy)  for  samples  type  B. 
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where  pb  is  the  bulk  density; 

Vp  is  the  volume  of  pores  per  gram; 

pt  is  the  structural  (true)  density. 

Applying  Eq.(2.17)  to  the  data  obtained  for  the  four  sets  of  samples 
of  12  A pore  radius,  one  obtains  the  variation  of  bulk  density  with 
densification  temperature  as  shown  in  Figure  2.38.  The  variation  of 
bulk  density  as  a function  of  densification  temperature  for  sets  of  silica 
gel  monolithic  samples  with  12  A,  32  A and  81  A average  pore  radii  is 
shown  in  Figure  2.39.  The  variation  of  volume  fraction  of  pores  with 
bulk  density  for  a set  of  12  A average  pore  radius  samples  is  shown  in 
Figure  2.40.  There  is  an  almost  linear  relationship  between  Vv  and 
bulk  density.  The  deviation  from  linearity  is  attributed  to  structural 
density  variations  as  shown  in  Figure  2.25.  In  order  to  check  for  the 
appropriateness  of  applying  Eq.(2.17)  in  obtaining  bulk  density  data,  a 
set  of  direct  measurements  of  bulk  density  through  mercury 
pycnometry  was  performed.  In  these  experiments  the  bulk  density  is 
evaluated  by  measuring  the  weight  of  the  sample  and  measuring  the 
volume  of  the  sample  by  weighing  the  volume  of  mercury  displaced 
when  the  sample  is  introduced  in  a container  filled  with  mercury.  The 
results  obtained  are  shown  in  Figure  2.41  for  a set  of  12  A average  pore 
radius.  It  can  be  seen  that  both  ways  of  evaluating  bulk  density  match 
relatively  well.  Since  in  the  gas  adsorption  experiments  the  nitrogen 
molecules  have  to  pass  through  the  pores  in  order  to  probe  the  surface 
and  volume  of  pores,  that  means  that  the  pores  have  to  be  open  for  such 
a penetration  of  gas  occur.  The  good  agreement  between  the  bulk 
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Figure  2.38  Bulk  density  versus  densification  temperature 

for  samples  type  A. 
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Figure  2.39  Bulk  density  versus  densification  temperature  for  silica  gel 
samples  with  12  A,  32  A and  81  A average  pore  radii. 
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Figure  2.40  Variation  of  volume  fraction  of  pores  (Vv)  as  a function  of 
bulk  density  for  silica  gel  samples  with  12  A average  pore  radius. 
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Figure  2.41  Variation  of  bulk  density  as  a function  of 
densification  temperature  for  silica  gel  samples  of  12  A 
average  pore  radius,  as  measured  by  nitrogen 
gas  adsorption  (N)  and  mercury  pycnometry  (Hg). 


/ 


1 1 1 

density  values  obtained  from  gas  adsorption  and  from  the  direct 
measurements  of  bulk  density  by  mercury  pycnometry,  indicates  that 
the  pore  network  in  the  silica  gel  monoliths  tested  is  interconnected. 
For  a structure  of  a large  volume  fraction  of  closed  pores  (non- 
interconnected  pores),  the  bulk  density  calculated  from  gas  adsorption 
should  be  greater  than  the  actual  value  of  bulk  density. 

The  variation  of  bulk  density  with  densification  temperature  for  8 1 
A average  pore  radius  samples,  as  measured  by  nitrogen  gas  adsorption 
and  mercury  pycnometry  is  shown  in  Figure  2.42.  For  temperatures  up 
to  900  °C  (1173  K)  there  is  a very  good  correspondence  between  those 
two  techniques.  At  higher  temperatures  the  results  obtained  from  the 
techniques  differ.  That  difference  is  likely  to  be  associated  with  pore 
closure  at  temperatures  higher  than  900  °C  (1173  K)  for  these  larger 
pore  radius  samples. 

Figure  2.43  shows  the  variation  of  volume  fraction  of  pores  with 
densification  temperature  for  samples  kept  at  the  densification 
temperature  for  periods  of  1 h and  10  h.  There  seems  to  be  only  small 
variations  in  Vv  between  those  two  different  processing  conditions. 
Figure  2.44  represents  the  variation  of  surface  area  of  pores  per  unit 
volume  with  densification  temperature  for  the  two  different  holding 
times  of  1 h and  10  h.  Figure  2.43  suggests  a small  decrease  in  Sv  when 
the  time  at  a given  sintering  temperature  increases.  Further 
information  on  the  effect  of  time  upon  sintering  can  be  seen  in  Figure 
2.45.  This  figure  shows  the  variation  of  volume  fraction  of  pores  with 
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Figure  2.42  Variation  of  bulk  density  as  a function  of 
densification  temperature  for  silica  gel  samples  of  8 1 A 
average  pore  radius,  as  measured  by  nitrogen  gas 
adsorption  (N)  and  mercury  pycnometry  (Hg). 
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Figure  2.43  Variation  of  volume  fraction  of  pores  (Vy)  as  a 
function  of  densification  temperature  for  samples  type  A 
kept  at  the  densification  temperature  for  1 h and  10  h. 
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Figure  2.44  Variation  of  the  surface  area  of  pores  per  unit  volume 
(Sy)  as  a function  of  densification  temperature  for  samples 
type  A kept  at  the  densification  temperature  for  1 h and  10  h. 
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Figure  2.45  Variation  of  the  volume  fraction  of  pores  (Vy)  as  a 
function  of  time  for  samples  type  A2  and  A3. 
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holding  time  at  800  °C  (1073  K),  for  two  sets  of  12  A average  pore 
radius  samples. 


2.7  Humidity  Effects  on  Densification 

In  order  to  evaluate  the  effect  of  atmosphere  in  the  sintering 
behavior  of  silica  gel  monoliths,  samples  were  densified  in  a humid, 
ambient  atmosphere  furnaces,  as  opposed  to  using  a dry-air  atmosphere 
inside  the  furnace.  The  variation  of  volume  fraction  of  pores  with 
densification  temperature  for  12  A average  pore  radius  samples  is 
shown  in  Figure  2.46  for  both  types  of  atmospheres.  The  figure 
suggests  that  the  humid  ambient  causes  a slightly  faster  decrease  in  Vv 
of  pores.  Similar  suggestion  is  obtained  from  Figure  2.47  which  shows 
the  variation  of  surface  area  of  pores  per  unit  volume  with  densification 
temperature  for  both  atmospheres. 

The  variation  of  bulk  density  with  densification  temperature  is 
shown  in  Figure  2.48  for  samples  type  A1  sintered  in  a dry-air 
atmosphere  and  in  a humid  atmosphere.  Figure  2.49  shows  similar 
curve  for  samples  type  A2.  Both  Figures  2.48  and  2.49  suggests  a 
slightly  faster  increase  in  bulk  density  with  temperature  for  samples 
processed  in  a humid  atmosphere. 

The  relatively  faster  densification  in  a humid  atmosphere  can  be 
associated  with  the  experimentally  observed  fact  that  samples  densified 
in  humid  air  will  bloat  or  foam  easier  (meaning  at  lower  temperatures) 
than  samples  densified  in  dry-air.  This  fact  is  associated  with  trapping 
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Figure  2.46  Variation  of  volume  fraction  of  pores  (Vy)  as  a 
function  of  densification  temperature  for  samples  type  A 
densified  in  a dry-air  (DA)  and  in  a humid  (OF)  atmospheres. 
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Figure  2.47  Variation  of  surface  area  of  pores  per  unit  volume  (Sy) 
as  a function  of  densification  temperature  for  samples  type  A 
densified  in  a dry-air  (DA)  and  in  a humid  (OF)  atmospheres. 
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Figure  2.48  Variation  of  bulk  density  as  a function  of  densification 
temperature  for  samples  type  A2  densified  in  a dry-air  (DA) 
and  in  a humid  (OF)  atmospheres. 
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Figure  2.49  Variation  of  bulk  density  as  a function  of  densification 
temperature  for  samples  type  A1  densified  in  a dry-air  (DA) 
and  in  a humid  (OF)  atmospheres. 
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of  water  inside  closing  pores.  When  the  temperature  is  raised  further 
the  increasing  pressure  inside  the  closed  pores  will  cause  the  material  to 
crack  and/or  foam  (or  bloat).  For  that  reason,  most  of  the  samples 
prepared  for  physical  testing  in  this  work  were  densified  in  a dry-air 
atmosphere.  This  resulted  in  dense  monolithic  samples  with  no  defects 
such  as  cracks  or  foaming. 


CHAPTER  3 

TOPOLOGICAL  MODELING  OF  DENSIFICATION 


3.1  Description  of  the  Structure  of  Materials 
3.1.1  The  Metric  Properties 

According  to  DeHoff  [357,  p.4633]  "in  biology,  function  follows 
form;  in  materials  science,  structure  determines  properties." 
Therefore,  the  study  of  the  structure  of  materials  is  the  central  part  of 
materials  science.  In  order  to  describe  completely  the  structure  of  a 
material,  two  basic  characteristics  have  to  be  defined: 

1)  The  constitution  of  the  structure,  in  terms  of  its  chemical 
composition  and  components  (e.g.,  its  thermodynamic  phases). 

2)  The  geometry  of  the  structure,  which  is  described  by  two  types  of 
parameters:  metric  and  topological  parameters. 

For  the  case  of  silica  gels,  for  instance,  the  composition  is  pure  SiC>2 
(for  the  dense  material  only,  because  in  the  early  stages  of  densification 
one  has  to  take  into  account  the  presence  of  water,  acidic  solutions, 
organics,  etc.)  and  it  is  monophasic  (SiC>2  glassy  phase  for  the  fully 
dense  gel-silica).  For  intermediate  stages  one  has  to  recognize  the 
presence  of  two  phases:  the  SiC>2  solid  network  and  the  pores. 
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Since  the  structure  of  a material  exists  in  a three-dimensional  space 
(Euclidean  space  E3),  the  structural  parameters  are  usually  reported  per 
unit  volume  of  the  material.  A three-dimensional  structure  can  be 
decomposed  into  elements  of  dimensionalities  varying  from  three  to 
zero  (tri-dimensional,  bi-dimensional,  uni-dimensional,  and  zero- 
dimensional elements). 

In  the  case  of  a polycrystalline  material,  the  presence  of  grain 
boundaries  makes  it  straightforward  to  associate  structural  features  to 
each  dimensionality,  as  shown  in  Table  3.1  [357].  In  a glassy  material 
there  will  be  no  grains,  no  grain  boundaries,  no  edges,  and  no  comers. 

According  to  DeHoff  [357]  the  structural  state  of  a material  can  be 
described  by  three  different  states: 

1)  The  qualitative  state,  which  is  a list  of  the  classes  of  each  kind  of 
features  (cells,  surfaces,  lines,  and  points)  that  exist  in  the 
structure. 

2)  The  quantitative  state,  which  associates  numbers  and  quantities  to 
the  structural  features  present. 

3)  The  topographic  state,  which  describes  non-uniformities  in  the 
structure  such  as  (a)  gradients  (variation  of  property  with 
position);  (b)  anisotropies  (variation  of  property  with 
orientation),  and  (c)  variations  in  spatial  distributions  (e.g., 
clustering,  ordering,  or  associations  or  disassociations  of  one 
class  of  features  with  others). 
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Table  3.1  Features  encountered  in  a polycrystalline  material  in 
relation  to  their  dimensionalities  [357]. 


Dimension 

Feature 

3 

Cell 

Grain 

2 

Surface 

Boundary 

1 

Line 

Edge 

0 


Point 

Comer 


125 


Also  according  to  DeHoff  [305,  p.292], 

those  properties  of  a structure  which  may  be  classed  as  topological 
include  any  property  which  is  unchanged  when  the  structure  is 
stretched,  distorted,  or  twisted  without  limit,  so  long  as  it  is  neither 
tom  nor  joined  to  itself  during  the  deformation.  Topological 
properties  thus  have  to  do  with  the  connectivity  of  a structure,  and 
are  independent  of  details  of  shape  and  size.  It  is  precisely  because 
topology  abstracts  geometric  properties  without  concern  about 
metric  detail  that  these  concepts  are  so  valuable  in  characterizing 
complex  structures. 

On  other  hand,  the  metric  properties  of  a structure  change  when  the 
structure  is  stretched,  distorted  or  twisted.  Examples  of  metric 
properties  include  [305,357-365]: 

1)  Volume  fraction  (Vy),  which  represents  the  volume  of  a given 
phase  or  constituent  divided  by  the  total  volume  of  the  material; 
it  can  be  evaluated  either  using  stereology  (via  area  fraction, 
lineal  fraction  or  point  counting)  or  through  direct  measurements 
of  volume  on  the  structure  (e.g.,  using  gas  adsorption  techniques). 

2)  Surface  area  per  unit  volume  (Sy),  which  represents  the  surface 
area  of  a given  phase  or  constituent  divided  by  the  total  volume  of 
the  material;  it  can  be  evaluated  either  using  stereology  (via 
number  of  intercepts  per  unit  of  length  or  number  of  features  per 
unit  of  length)  or  through  direct  measurements  of  surface  on  the 
structure  (e.g.,  using  gas  adsorption  techniques). 

3)  Total  curvature  of  a surface  per  unit  volume  (Mv),  where  the 
total  curvature  is  defined  as  the  sum  of  the  surface  curvature  plus 
the  edge  curvature;  it  can  be  evaluated  stereologically  by 
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counting  either  the  number  of  tangents  per  unit  area  or  the 
number  of  features  per  unit  area. 

4)  Line  length  per  unit  volume  (Lv),  which  can  be  evaluated 
stereologically  by  counting  the  number  of  quadruple  points  per 
unit  area. 

5)  Integral  curvature  of  a surface  per  unit  volume  (Iv),  which 
corresponds  to  the  net  spherical  image  of  the  surface  in  a unit 
volume. 

6)  Integral  curvature  of  a line  per  unit  volume  (Cv). 

7)  Integral  torsion  of  a line  per  unit  volume  (Fv). 

Besides  those  seven  main  metric  parameters,  there  are  average 
properties  (or  combinations  of  metric  properties)  which  are  very  useful 
for  certain  applications.  Some  important  average  parameters  are  the 
following  [357,365-370]: 

1)  Mean  lineal  intercept  (A.3),  which  represents  the  average  surface- 
to-surface  distance  through  a cell,  and  for  each  type  of  surface  it 
is  given  by 


2)  Mean  free  lineal  intercept  (A,f),  which  represents  the  average 
surface-to-surface  distance  through  the  other  (or  others) 
constituent,  and  for  a given  type  of  surface  it  is  given  by 
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(3.2) 


3)  Average  mean  curvature  (H),  given  by 


(3.3) 


4)  DeHoff  s shape  parameter  (De),  defined  by 


De  = H\3  . 


(3.4) 


3.1.2  The  Topological  Properties 


According  to  Stillwell  [371,  p.2] 

topology  is  the  branch  of  geometry  which  studies  the  properties  of 
figures  under  arbitrary  continuous  transformations.  Just  as  ordinary 
geometry  considers  two  figures  to  be  the  same  if  each  can  be  carried 
into  the  other  by  a rigid  motion,  topology  considers  two  figures  to 
be  the  same  if  each  can  be  mapped  onto  the  other  by  a one-to-one 
continuous  function.  Such  figures  are  called  topologically 
equivalent,  or  homeomorphic. 

Those  properties  of  topological  spaces  that  are  invariant  under 
homeomorphisms  are  termed  topological  properties  [372]. 

Topology  is  divided  into  two  main  branches  [373]:  (1)  general 
topology  (or  point-set  topology)  and  (2)  algebraic  topology  (or 
combinatorial  topology).  While  general  topology  deals  with  set- 
theoretic  tools,  algebraic  topology  uses  algebra  (mainly  group  theory). 
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The  first  known  work  on  topology  was  carried  out  by  Euler  in  1752 
[374],  Listing  [375]  seems  to  have  been  the  first  to  use  the  word 
"topology"  in  1847.  Besides  Listing,  further  pioneering  work  on 
topology  was  done  by  Cantor  [376],  Riemann  [377]  and  Poincare  [378] 
by  the  end  of  the  19th  Century,  and  by  Riesz  [379],  Frechet  [380]  and 
Moore  [381]  in  the  early  1900s.  The  pioneering  work  of  introducing 
the  topological  concepts  in  materials  science  was  carried  out  by  Rhines 
[284]  in  1958.  Details  on  set-point  topology  and  algebraic  topology  are 
found  in  an  extensive  list  of  publications  [e.g.,  382-390]. 

Algebraic  topology  study  surfaces  by  "triangulating"  thep,  that  is,  by 
dividing  them  up  into  triangular  regions  in  such  a way,  that  two  regions 
which  intersect  do  so  by  either  contacting  along  a single  edge  or  at  a 
single  comer  (or  vertex)  of  a triangle.  These  rectilinear  triangulations 
are  termed  simplicial  complexes.  Giblin  [391,  p.7]  points  out  that  the 
advantage  of  reducing  a surface  to  simplicial  complexes  "lies  in  the  fact 
that  any  triangulated  surface  is  homeomorphic  to  one  in  which  the 
triangles  are  flat  and  the  edges  straight,"  which  allows  for  a rigorous 
algebraic  study  of  surfaces. 

Each  triangle  has  three  straight  edges  and  three  vertices.  The  so- 
called  intersection  conditions  that  are  imposed  on  the  geometric 
arrangement  of  the  triangles  are  [391] 

Two  triangles  either 

1)  Are  disjoint  (Figure  3.1(a)); 

2)  Have  one  vertex  in  common  (Figure  3.1(b)); 

3)  Have  two  vertices  in  common  (Figure  3.1(c)). 
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(a) 


Figure  3.1  Intersection  conditions. 
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Based  on  the  above,  a closed  surface  M is  defined  such  that  [391] 

1)  M satisfies  the  intersection  condition; 

2)  M is  connected; 

3)  For  every  vertex  v of  a triangle  of  M;  the  link  of  v is  a simple 
closed  polygon. 

For  a closed  surface  M,  the  Euler  characteristic  (x)  of  M is  defined 
as  [389] 


X(M)  = V - E + T , (3.5) 

where  V is  the  number  of  vertices  of  M; 

E is  the  number  of  edges  of  M; 

T is  the  number  of  triangles. 

The  Euler  characteristics  of  various  types  of  surface  are  shown  in 
Table  3.2  [389].  Figure  3.2  [390]  shows  schematics  of  the  geometric 
shapes  mentioned  in  Table  3.2  (sphere,  toms  and  Klein  bottle). 

Assuming  W as  a homology  group  of  a finite  complex  associated 
with  a topological  space,  a crude  measure  of  the  size  of  W is  given  by 
the  rank  (hk)  of  W,  or  the  number  of  elements  in  W.  The  numbers  hk 
are  termed  the  Betti  numbers  of  the  complex.  The  Betti  numbers  are 
described  by  the  symbol  Bi,  where  i is  an  integer  varying  from  0 to  n, 
and  n is  the  dimensionality  of  the  complex.  Table  3.3  shows  the 
possible  Betti  numbers  for  complexes  with  dimensionalities  varying 
from  0 to  3. 
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Table  3.2  Euler  characteristic  (%)  of  the  surface  M [389]. 

Surface  %(M) 

Sphere  2 

Connected  sum  of  n tori  2 - 2n 

Connected  sum  of  n projective  planes  2 - n 

Connected  sum  of  projective  plane  and  n tori  1 - 2n 

- 2n 


Connected  sum  of  Klein  bottle  and  n tori 
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(c) 


Figure  3.2  (a)  Sphere;  (b)  torus;  (c)  Klein  bottle  [390]. 
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Table  3.3  Betti  numbers  of  complexes  with  dimensionalities 

varying  from  0 to  3. 


Dimensionality  Betti  numbers 


0 


6 


o 


1 Bo,  Bi 

2 B0,  61,62 


3 


B0,  Bi,  62,  63 
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The  zeroth  Betti  number  (B0)  is  the  number  of  components  or 
disjoined  parts  of  a complex  [392].  The  first  Betti  number  (Bi)  is  called 
the  connectivity  number  of  the  surface  [392].  It  represents  the  largest 
number  of  closed  curves  that  can  be  drawn  on  the  surface  of  a complex 
without  dividing  the  surface  into  two  or  more  pieces.  The  second  Betti 
number  (B2)  is  one  for  all  compact,  connected  surfaces  [392].  The  third 
Betti  number  (B3)  is  always  zero  for  finite  complexes  imbedded  in  the 
Euclidean  space  E3  [392]. 

According  to  Rhines  et  al.  [363,  p.8], 

the  connectivity  C of  an  isolated  volume  is  defined  as  the  number  of 
independent  closed  curves  that  may  be  constructed  within  the  volume 
which  cannot  be  smoothly  deformed  into  each  other,  nor  shrunk  to  a 
point,  without  leaving  the  volume.  Thus  a sphere  has  a connectivity 
of  zero  since  any  closed  curve  constructed  within  its  volume  may  be 
shrunk  to  a point;  the  sphere  is  a member  of  a class  of  volumes  or 
bodies  called  simply  connected  bodies,  all  of  which  may  be  stretched 
or  deformed  into  a sphere. 

Figure  3.3  shows  a few  geometric  shapes  that  are  topologically 
equivalent  to  a sphere.  A torus  has  a connectivity  of  one,  since  one  can 
draw  a curve  (around  the  hole)  which  cannot  be  shrunk  to  a point 
without  leaving  the  volume  of  the  torus.  Figure  3.4  [390]  shows  a few 
three-dimensional  geometric  shapes  which  have  connectivity  of  one,  and 
therefore,  are  topologically  equivalent  to  a torus. 

For  sintered  structures,  the  connectivity  assumes  very  high  values, 
specially  in  the  initial  stages  of  the  sintering  process.  As  the  pores 
disappear  during  sintering,  the  connectivity  decreases.  Both  the  solid 
phase  and  the  gas  phase  in  a sintered  structure  have  the  same 
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Figure  3.3  Schematic  of  three-dimensional  geometric  shapes 
with  zero  connectivity  (C). 
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Figure  3.4  Schematic  of  three-dimensional  geometric  shapes  with 
connectivity  (C)  equal  to  one  [390]. 
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connectivity.  They  are  termed  [363]  conjugate  to  each  other,  because 
together  they  fill  space. 

While  the  connectivity  (C)  is  a property  of  the  volume,  the  genus 
(G)  is  a property  of  a closed  surface.  DeHoff  [305,  p.296]  defines  genus 
as  "the  maximum  number  of  non-self  reentrant  closed  curves  which 
may  be  constructed  on  the  surface  without  dividing  it  into  two  separate 
parts."  Another  equivalent  way  of  defining  genus  is  given  by  Massey 
[389]:  "a  surface  that  is  the  connected  sum  of  n tori  or  n projective 
planes  is  said  to  be  of  genus  n." 

According  to  Rhines  et  al.  [363,  p.8], 

it  can  be  shown  that  the  genus  of  a closed  surface  is  equal  to  the 
connectivity  of  either  of  the  two  volumes  (solid  and  gas)  of  space 
which  it  separates.  Thus,  for  a sintered  structure,  the  genus  of  pore- 
solid  interface,  the  connectivity  of  the  pore  phase,  the  connectivity  of 
the  solid  region,  . . . and  the  first  Betti  number,  are  all  equal. 

Therefore,  the  connectivity  of  the  pore  network  in  a sintered  gel-silica 

can  be  expressed  by  either  CyP  (connectivity  of  the  pores),  Cys 

(connectivity  of  the  solid  phase),  Gy  (genus  of  the  solid-pore  interface) 

or  by  6 1 (the  first  Betti  number  of  the  pore-solid  interface),  and  all 

those  parameters  possess  the  same  numerical  value. 

Since  a structure  is  usually  defined  in  the  Euclidean  space  E3 

(because  this  is  a three-dimensional  world),  it  is  useful  to  list  the 

topological  parameters  pertaining  to  each  dimensionality  varying  from 

zero  to  three  as  shown  in  Table  3.4.  Zero-dimensional  structures  (or 

zero-dimensional  components  of  a structure)  can  be  characterized  by  the 

number  of  zero-dimensional  entities  (N°),  and  there  is  no  connectivity 
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Table  3.4  Independent  topological  parameters  that  can  be 
associated  to  three-dimensional  structures. 

Dimensionality  Parameter 


0 N° 

1 N1 , |i 

2 N2,  G 


3 


N3,  C 
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associated  with  those  zero-dimensional  features.  One-dimensional 
structures  (or  one-dimensional  components  of  a structure)  can  be 
described  by  the  number  of  one-dimensional  features  (N1),  and  by  the 
so-called  "cyclomatic  number"  (ji),  which  is  given  by 

Ji  = b - n + p , (3.6) 

where  b is  the  number  of  branches; 
n is  the  number  of  nodes; 

p is  the  number  of  components  or  disjoined  parts  of  the  uni- 
dimensional network. 

Two-dimensional  structures  (or  two-dimensional  components  of  a 
structure)  can  be  described  by  the  number  of  two-dimensional  features 
(N2),  and  by  the  genus  (G)  of  the  two-dimensional  constituents.  Three- 
dimensional  structures  (or  three-dimensional  components  of  a structure) 
are  characterized  by  the  number  of  three-dimensional  features  (N3),  and 
by  the  connectivity  (C)  of  the  three-dimensional  constituents.  For  all 
dimensionalities  (0,  1,  2 and  3)  of  a structure,  the  number  of  features  is 
given  by  the  zeroth  Betti  number  (B0).  The  first  Betti  number  (Bi) 
corresponds  to  p,  (1-D),  G (2-D),  or  C (3-D).  A general  expression  for 
the  first  Betti  number  (Bi),  which  includes  Eqs.(3.5)  and  (3.6)  can  be 
written  as 


Bj  — B - N + Bq  , 


(3.7) 
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where  B is  the  number  of  branches; 
N is  the  number  of  nodes. 


3.1.3  Evaluation  of  Topological  Properties 


Figure  3.5  shows  a schematic  of  a structure  formed  by  an  array  of 
points  (nodes)  and  lines  (branches).  For  this  structure  one  can  count: 

B = 7, 

N = 7, 

Bo  = 2, 

therefore,  the  first  Betti  number  (Bi)  of  that  structure  as  given  by 
Eq.(3.7)  is 

Bi  = p.  = 7-  7 + 2 = 2,  (3.8) 

meaning  that  it  is  possible  to  cut  two  branches  of  that  structure  without 
increasing  the  number  of  separate  parts. 

Considering  the  two-dimensional  structure  shown  in  Figure  3.6(a), 
one  can  count: 

B = 6, 

N = 4, 

Bo=  1, 
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Figure  3.5  Schematic  of  a structure  formed  by  two  linear  features. 
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Figure  3.6  Two-dimensional  structure  containing  three  holes. 
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therefore,  the  first  Betti  number  (Bi)  of  that  structure  as  given  by 
Eq.(3.7)  is 


B1=G  = 6-4+l  = 3,  (3.9) 

meaning  that  it  is  possible  to  make  three  cuts  through  the  structure 
without  dividing  it  into  two  (or  more)  parts.  This  can  be  directly  seen 
by  inspecting  Figure  3.6(b),  which  shows  three  possible  cuts  in  the 
structure  without  separating  it.  Any  additional  cut  will  divide  the 
structure  into  at  least  two  parts.  Another  way  of  obtaining  the 
connectivity  of  three  is  by  directiy  counting  the  number  of  holes  present 
in  the  structure  (3  holes). 

For  three-dimensional  constituents  of  a structure,  one  can  describe 
the  connectivity  (C)  by  evaluating  the  genus  (G),  because  of  the 
numerical  equality  between  connectivity  anj  genus.  Considering  Figure 
3.7,  one  obtains 
B = 11, 

N = 7, 

Bo=  2; 

therefore,  Eq.(3.7)  yields 

B1=G  = C=  11 - 7 + 2 = 6,  (3.10) 

meaning  that  it  is  possible  to  make  six  cuts  in  that  structure  without 
increasing  the  number  of  separate  parts. 
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Figure  3.7  Schematic  of  a three-dimensional  structure 
formed  by  two  separate  parts. 
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In  the  examples  shown  above  it  is  very  easy  to  count  directly  the 
parameters  B,  N and  60-  hi  a real  structure  this  process  is  extremely 
difficult,  because  of  the  large  numbers  involved,  the  small  scale  of  the 
structure  of  materials,  and  because  of  the  inherent  problem  of 
identifying  a branch,  a node,  and  whether  two  parts  of  a structure  are 
indeed  separated  from  each  other. 

In  materials  science  and  biology,  one  needs  to  use  the  technique  of 
serial  sectioning  in  order  to  be  able  to  directly  count  N,  B and  B0-  This 
technique  is  a series  of  direct  observations  of  a polished  section  of  a 
structure  in  such  a way  that  after  preparing,  observing,  and  recording 
one  section,  one  polishes  out  a small  layer  of  the  material  and  observe 
the  structure  just  underneath  the  previously  observed  one.  This  process 
is  repeated  several  times  until  it  allows  for  a three-dimensional 
reconstruction  of  the  network  of  branches  and  nodes  based  upon  the 
combined  analysis  of  the  various  sections.  In  order  to  calculate  the 
topological  parameters  in  a unit  volume,  accurate  measurements  of  the 
thicknesses  of  the  removed  layers  have  to  be  carried  out.  Examples  of 
applications  of  the  serial  sectioning  technique  are  given  by  Aigeltinger 
[309]  and  DeHoffet  al.  [393]. 

The  application  of  the  technique  of  serial  sectioning  is  a time- 
consuming  and  tedious  process.  The  fundamental  limitations  of  the 
technique  is  that  one  has  to  be  able  to  not  only  see  the  structure,  but  to 
carefully  remove  a precise  layer  of  the  material  and  look  again  exactly 
underneath  the  previously  observed  one.  As  DeHoff  et  al.  [393]  and 
Aigeltinger  [309]  proved,  it  is  feasible  to  carry  out  that  task  for 
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micrometer  scale  structures.  However,  in  the  case  of  the  evaluation  of 
topological  properties  of  alkoxide  derived  silica  gels  such  a task  is  not 
yet  feasible.  The  difficulties  in  applying  serial  sectioning  to  alkoxide 
silica  gels  arise  because  of  the  small  scale  of  the  pore  structures 
(average  pore  radius  varying  from  12  A to  81  A,  as  shown  in  Chapter 
2).  There  are  high-resolution  transmission  electron  microscopes  with 
small  enough  resolution  to  identify  pores  in  the  nanometer  scale. 
However,  as  pointed  out,  the  application  of  the  technique  requires  the 
precise  elimination  of  layers  of  well  defined  thicknesses  (which  should 
be  of  few  nanometers  for  the  case  of  silica  gels)  and  the  observation  of 
the  structure  right  underneath  the  previously  recorded  observation.  All 
these  requirements  does  not  make  it  feasible  to  apply  the  serial 
sectioning  technique  to  alkoxide  derived  silica  gel  monoliths. 

Estimates  of  topological  parameters  can  be  carried  out  through 
modeling  the  structure  and  assuming  fixed  shapes  for  those  structures 
[394-398]. 

In  the  following  sections  the  modeling  of  the  pore  structure  of  silica 
gel  monoliths  will  be  described.  The  goal  of  such  a modeling  is  the 
obtaining  of  topological  parameters  that  characterize  the  pore  structure, 
so  that  a complete  description  of  the  structure  involving  both  types  of 
geometric  properties  (metric-Chapter  2,  and  topological)  becomes 
possible. 
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3.2  Topological  Modeling 


3.2.1  Introduction 

When  considering  the  possibilities  of  feasible  shapes  of  pores  to  be 
used  as  models  for  evaluation  of  topological  parameters,  the  limiting 
cases  are  the  triangular  shape  and  the  circular  shapes  as  seen  in  Figure 
3.8.  The  triangular  shape  has  a minimal  number  of  sides  (3)  while  the 
circular  shape  could  represent  any  polygon  with  an  infinite  number  of 
sides.  Any  other  geometric  shape  has  a number  of  sides  that  fall  in 
between  those  two. 

Assuming  the  triangular  shape  to  represent  the  cross  section  of  the 
pores  of  the  structure,  a three-dimensional  view  of  the  pore  structure 
would  look  like  the  schematic  shown  in  Figure  3.9,  when  the 
coordination  number  of  pores  is  equal  to  four.  Notice  that  the  shape 
formed  when  four  branches  of  the  triangular  pores  meet  is  a 
tetrahedron.  Therefore,  this  model  can  be  referred  to  as  tetrahedral 
geometric  model,  or  triangular  model. 

When  the  circle  is  assumed  to  represent  the  cross  section  of  the 
pores,  a cylindrical  model  is  created,  as  shown  schematically  in  Figure 
3.10. 

In  the  next  sections  tetrahedral,  cylindrical  and  spherical  models  are 
described.  The  models  are  numbered  as  follows:  Model  1 assumes 
spherical  geometry;  Models  2,  3 and  4 assume  Pv  = 1 (Pv  is  the  number 
of  separate  parts  per  unit  volume)  and  the  tetrahedral  geometry;  Models 
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Figure  3.8  Fundamental  shapes  used  for  geometric  modeling. 
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Figure  3.9  Tetrahedral  geometric  model  assumed 
for  topological  modeling. 
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Figure  3.10  Cylindrical  geometric  model  assumed 
for  topological  modeling. 
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5 and  6 assume  Pv  = 1 and  cylindrical  geometry;  Model  7 assumes  a 
variable  Pv,  a constant  Nv  (Nv  is  the  number  of  nodes  per  unit  volume), 
and  cylindrical  geometry;  Models  8,  9 and  10  assume  a variable  Pv  and 
tetrahedral  geometry;  and  Models  11  and  12  assume  a variable  Pv  and 
cylindrical  geometry. 


3.2.2  Tetrahedral  Modeling 


As  shown  in  Figure  3.11  the  tetrahedral  geometric  model  can  be 
decomposed  into  a set  of  branches  and  nodes.  One  can  associate  with 
each  branch  and  each  node  a volume  and  an  external  surface  which 
directly  depend  on  the  size  of  the  branch  (1)  and  the  height  of  the 
triangle  (h). 

The  total  volume  of  a unit  formed  by  one  node  and  one  branch  (Vu) 
is  given  by 


Vu 


V3  h2  (h  + 6 1) 
8 


(3.11) 


The  total  surface  area  of  a unit  formed  by  one  node  and  one  branch 
(Su)  is  given  by 


S 


u 


9 h 1 
" <2  • 


(3.12) 
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Figure  3.11  Structural  components  of  the  tetrahedral 
geometric  model. 
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Assuming  that  one  can  obtain  experimental  evaluations  of  the  total 
volume  of  pores  and  of  the  total  surface  area  of  pores  (using  BET  data 
for  instance),  it  is  possible  to  associate  the  total  volume  and  surface  of 
structural  units  (based  on  the  tetrahedral  model)  and  then  estimate  what 
the  total  number  of  units  (N)  should  be. 

The  ratio  between  the  volume  and  the  surface  of  a structural  unit  is 
given  by 


Vu  h2  V6  h 
Su  " 72  1 + 12  * 


(3.13) 


The  total  volume  of  pores  (Vt)  and  surface  area  of  pores  (St)  can  be 
written  as 


VT  = N Vu  , 


(3.14) 


and 

ST  = N Su  . (3.15) 

Using  Eq.(3.14)  and  Eq.(3.15),  Eq.(3.13)  can  be  rewritten  as 

Vt  _ V6  h2  V6  h 
ST  ' 72  1 + 


12  * 


(3.16) 
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Equation  (3.16)  can  be  expressed  in  terms  of  volume  fraction  of 
pores  (Vv)  and  surface  area  of  pores  per  unit  volume  (Sv).  Replacing  h 
by  D (average  diameter  of  pores)  and  1 by  L (average  branch  size), 
Eq.(3.16)  can  finally  be  written  as 

Vv  V6  D2  V6  D 

Sv  “ 72  L + 12  ' (3-17) 

Using  the  same  notation  in  Eq.(3.11)  and  Eq.(3.13),  and  Eq.(3.14) 
and  Eq.(3.15)  one  obtains 

V3 

Vv  = ^-NvD2(D  + 6L)  (3.18) 

and 

9 

Sv  = Nv  D L (3.19) 

where  Nv  is  the  number  of  nodes  per  unit  volume. 

Since  Eqs.(3.18)  and  (3.19)  have  three  unknown  parameters  (Nv,  D 
and  L),  it  is  necessary  to  assume  values  for  one  of  them  (D  or  L)  in 
order  to  solve  the  equations  and  obtain  Nv. 

Some  possibilities  for  solving  Eqs.(3.18)  and  (3.19)  are  shown  in 
Table  3.5.  Each  type  of  solution  gives  rise  to  a different  model,  which 
is  numbered  together  with  the  corresponding  assumption.  The  value  of 
Nv  can,  therefore,  be  obtained  by  using  either  Eq.(3.18)  or 
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Table  3.5  Set  of  assumptions  made  to  solve  the  equations 
corresponding  to  the  tetrahedral  geometric  model. 


Model  2 

Model  8 

L = kD 

L = kD 

Py  = 1 

Pv  = variable 

Model  3 

Model  9 

4(1-  Vy  ) 
L " Sy 

T 4(1-  Vy  ) 

L ~ Sy 

Py  = 1 

Py  = variable 

Model  4 

Model  10 

4 Vy 

D=  SV 

4 Vy 
D=  SV 

Py  - 1 

Pv  = variable 
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Eq.(3.19),  as  follows: 


Nv  V3  D2  (D  + 6 L)  ’ 


8 Vv 


(3.20) 


or 


(3.21) 


According  to  Table  3.5,  Models  2,  3 and  4 assume  that  the  number 
of  separate  parts  (Py)  is  equal  to  1 during  the  entire  process.  The  main 
implication  of  this  assumption  is  that  during  the  last  stage  of 
densification,  when  the  genus  is  equal  to  one,  the  number  of  nodes  has 
to  be  equal  to  the  number  of  branches  because  of  the  general  expression 


When  By  is  equal  to  Ny  the  coordination  number  of  pores  (CN)  is 
equal  to  two,  as  shown  in  Figure  3.12. 

Figure  3.13  shows  a schematic  of  the  pore  coordination  number 
varying  from  4 to  0. 

As  shown  in  Figure  3.12,  Models  2,  3 and  4 assume  four  as  the 
average  coordination  number  of  pores  (CN)  in  the  dried  stage.  That 
initial  stage  usually  corresponds  to  a maximum  number  of  branches, 
nodes  and  genus.  As  the  genus  decreases,  during  the  second  stage  of 


Gy  = 1 = By  - Ny  + 1 . 


(3.22) 
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PORE  COORDINATION  NUMBER 


Figure  3.12  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (Bv)  and  genus  (Gy)  as  a function  of  the  pore 
coordination  number  (CN)  when  Pv  = 1. 
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Figure  3.13  Schematic  of  the  pore  coordination  number  (CN) 

varying  from  4 to  0. 
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densification,  the  number  of  nodes  remains  constant  and  the  number  of 
branches  decreases.  Models  2,  3 and  4 assume  that  when  CN  reaches  2, 
further  elimination  of  branches  leads  to  disappearance  of  nodes,  and 
therefore  Bv  and  Nv  decrease,  keeping  CN  constant  at  2 during  the 
third  stage  of  densification.  At  this  stage  Gv  is  equal  to  1 and  is  kept 
constant  during  the  third  stage. 

Models  8,  9 and  10  have  similar  assumptions  to  Models  2,  3 and  4, 
except  that  they  now  consider  the  general  expression 

Gy  = By  - Ny  + Py  , (3.23) 

and  therefore  Pv  can  assume  values  different  from  1 in  the  last  stage  of 
densification.  Therefore,  CN  goes  from  4 to  0,  as  shown  in  Figure 
3.14. 

According  to  Figure  3.14,  Models  8,  9 and  10  assume  4 as  the 
average  coordination  number  of  pores  (CN)  in  the  dried  stage.  As  the 
genus  decreases,  during  the  second  stage  of  densification,  the  number  of 
nodes  remains  constant  and  the  number  of  branches  decreases.  The 
Models  8,  9 and  10  assume  that  when  CN  reaches  2,  the  number  of 
separate  parts  (Pv)  starts  increasing,  the  number  of  nodes  (Nv) 
decreases,  and  the  number  of  branches  (By)  follows  a different  path 
compared  to  its  initial  trajectory. 

In  order  to  evaluate  By,  Nv,  Gy  and  Py  experimentally  during  the 
third  stage  of  densification,  it  is  assumed  that  the  number  of  nodes 
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PORE  COORDINATION  NUMBER 


Figure  3.14  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (Bv),  number  of  separate  parts  (Pv)  and  genus  (Gy)  as  a 
function  of  the  pore  coordination  number  (CN)  when  Pv  is  variable. 
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follows  a linear  relationship  with  CN  as  shown  in  Figure  3.14.  Nv  is 
estimated  by  using  the  expression: 


Nv  = 


V2Sy 

9 D L’ 


(3.24) 


and  the  coordination  number  is  given  by 


CN 


2 Ny 

NyO  ’ 


(3.25) 


where  Nv°  is  the  number  of  nodes  present  in  the  dried  stage,  or  the 
number  of  nodes  in  the  plateau  from  CN  4 to  2. 

The  number  of  branches  per  unit  volume  (By)  can  be  estimated  from 
Eqs.(3.24)  and  (3.25)  as  follows: 


By  = 


CN  Ny 
2 


(3.26) 


Since  Gy  is  equal  zero  in  the  third  stage  of  densification,  Pv  can  be 
calculated  as  follows: 

Gy  = By  - Ny  + Py  = 0 , (3.27) 

thus, 


Py  = Ny  - By  . 


(3.28) 
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3.2.3  Cylindrical  Modeling 


The  cylindrical  geometric  model  is  based  on  the  assumption  that  all 
the  volume  and  all  the  surface  area  is  associated  with  the  cylindrical 
branches  of  the  structure.  In  this  case  a node  represents  the  point 
where  the  branches  meet  and  it  has  no  volume  associated  to  it,  like  the 
node  in  the  tetrahedral  model  does. 

Considering  D as  the  diameter  of  the  cylinder  and  L as  the  average 
size  of  a cylindrical  branch,  the  volume  of  a branch  (or  unit,  Vb)  is 
given  by 


Vb 


K D2  L 
4 


(3.29) 


while  the  surface  area  of  a branch  (Sb)  is  given  by 


Sb  = 7t  D L . 


(3.30) 


Replacing  Vb  by  Vv,  and  Sb  by  Sv,  and  using  Bv  as  the  number  of 
branches  per  unit  volume,  one  can  write  Eq.(3.29)  and  Eq.(3.30)  as 


Vv  = 


7t  By  D2  L 
4 


(3.31) 


and 
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Sy  — 7C  By  D L . 


(3.32) 


One  can  obtain  either  directly  from  the  cylindrical  shape  or  from 
Eq.(3.31)  and  Eq.(3.32)  that  the  value  of  D is  always  a fixed 
relationship  between  Vv  and  Sv: 


D = 


4 Vy 
Sy 


(3.33) 


Therefore,  to  obtain  the  value  of  By  from  Eq.(3.31)  or  Eq.(3.32)  the 
only  unknown  is  the  parameter  L.  Once  L is  estimated  the  number  of 
branches  per  unit  volume  is  obtained  by  either  Eqs.(3.31)  or  (3.32); 
i.e.. 


B 


4 Vy 

v ~ 7t  D2  L ’ 


(3.34) 


or 


By 


4 Sy 
k D L ' 


(3.35) 


Table  3.6  shows  a set  of  assumptions  of  "L"  and  Pv  made  in  order  to 
solve  Eqs.(3.34)  and  (3.35). 
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Table  3.6  Set  of  assumptions  made  to  solve  the  equations 
corresponding  to  the  cylindrical  geometric  model. 


Model  7 

_ 4 (1  - Vy) 
Sy 

Pv  = variable 
Ny  = Nv° 


Model  5 


Model  11 


L = k D 


L = k D 


Pv  = 1 


Pv  = variable 


Model  6 


Model  12 


L 


4 (1  - Vy) 
Sv 


Pv  = 1 


4(1  - Vy) 
Sy 


Pv  = variable 
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3.2.4  Spherical  Modeling 


Even  though  a spherical  model  is  not  adequate  to  describe  the  pore 
structure  of  a silica-gel  monolith,  a spherical  model  is  created  in  order 
to  serve  as  reference  or  limiting  case  for  the  other  models. 

Assuming  the  spherical  pores  as  having  an  average  diameter  D,  the 
number  of  pores  per  unit  volume  (Nv)  can  be  directly  correlated  with 
the  volume  fraction  of  pores  (Vv)  as  follows: 


M 6VV 

Nv  = n D3 


(3.36) 


Ny  can  equally  be  associated  with  the  surface  area  per  unit  volume 
(Sy)  by  the  following  expression: 


Nv  = 


Sy 

71  D2  ' 


(3.37) 


The  average  diameter  of  a spherical  pore  is  directly  obtained  from 
Eqs.(3.36)  and  (3.37): 

6 Vv 

D - "sT-  (3.38) 

Therefore,  by  experimental  determination  of  Vv  and  Sv  it  is  possible 
to  obtain  the  number  of  "spherical  pores,"  Nv. 
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3.3  Application  of  Topological  Modeling  to  Gel-Silica 
3.3.1  Small  Pore  Size  Structures 

The  topological  model  as  seen  in  Figure  3.14  is  applied  to  samples 
type  A (12  A average  pore  radius)  using  the  cylindrical  model  with  a 
variable  Pv.  The  result  is  seen  in  Figure  3.15,  where  Bv,  Nv,  Py  and 
Gy  are  plotted  as  a function  of  the  pore  coordination  number  (CN). 
When  the  topological  parameters  are  plotted  as  a function  of 
temperature  the  results  are  shown  in  Figure  3.16,  for  the  cylindrical 
model  with  a variable  Pv.  Figure  3.16  thus  represents  the  temperature 
dependence  of  the  structural  evolution.  The  initial  situation  in  Figure 
3.16  is  at  a temperature  of  180  °C  (453  K),  corresponding  to  the  dried 
stage.  As  the  densification  temperature  increases  both  the  number  of 
branches  (By)  and  the  genus  (Gy)  decrease,  while  the  number  of  nodes 
(Nv)  remains  constant  at  Nv°  (number  of  nodes  in  the  dried  stage),  and 
Py  remains  constant  at  one.  The  position  where  the  pore  coordination 
number  decreases  to  2 (in  Figure  3.15)  corresponds  in  Figure  3.16  to 
the  point  where  By  reaches  the  value  of  Nv°,  which  occurs  at  a 
temperature  just  above  1000  °C  (1273  K).  At  this  point  the  genus  goes 
to  zero  and  the  number  of  separate  parts  increases  from  one. 

The  application  of  the  cylindrical  model  with  variable  Pv  (Model  12) 
to  other  type  A samples  (12  A average  pore  radius)  is  shown  in  Figures 
3.17,  3.18  and  3.19.  They  all  show  a very  similar  type  of  structural 
evolution.  The  figures  representing  experimental  data  obtained  in  this 
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PORE  COORDINATION  NUMBER 


Figure  3.15  Variation  of  the  number  of  nodes  (Ny),  number  of 
branches  (Bv),  number  of  separate  parts  (Pv)  and  genus  (Gy)  as  a 
function  of  the  pore  coordination  number  (CN)  as  assumed  in  the 
topological  modeling.  The  points  on  the  curves  show  where  the 
experimental  data  obtained  from  samples  type  A are  located  (Model  12). 
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Figure  3.16  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (Bv),  number  of  separate  parts  (Py)  and  genus  (Gy)  as  a 
function  of  densification  temperature  for  samples  type  A1  (Model  12). 
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Figure  3.17  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (Bv),  number  of  separate  parts  (Py)  and  genus  (Gy)  as  a 
function  of  densification  temperature  for  samples  type  A2  (Model  12). 
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Figure  3.18  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (Bv),  number  of  separate  parts  (Pv)  and  genus  (Gv)  as  a 
function  of  densification  temperature  for  samples  type  A3  (Model  12). 
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Figure  3.19  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (Bv),  number  of  separate  parts  (Pv)  and  genus  (Gv)  as  a 
function  of  densification  temperature  for  samples  type  A4  (Model  12). 
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work  are  obtained  either  using  a polynomial  or  simply  connecting  the 
average  points,  in  which  case  local  peaks  or  valleys  in  the  curves  do  not 
necessarily  represent  an  actual  trend  in  the  data.  The  temperature  in 
which  the  genus  goes  to  zero  for  all  the  type  A samples  is  in  the  range 
900  °C  (1173  K)  to  1100  °C  (1373  K). 

3.3.2  Large  Pore  Size  Structures 

Applying  Model  12  (cylindrical  model  with  variable  Pv),  the 
variation  of  number  of  separate  parts  (Pv)  and  genus  (Gv)  for  samples 
type  B (32  A average  pore  radius)  as  a function  of  densification 
temperature  is  shown  in  Figure  3.20.  This  figure  shows  an  increase  in 
Bv  and  Gv  as  the  densification  temperature  goes  from  180  °C  (453  K) 
until  900  °C  (1173  K)  when  both  Bv  and  Gv  begin  to  decrease.  The 
genus  goes  to  zero  at  a temperature  of  1250  °C  (1523  K). 

For  larger  pore  structures  (81  A average  pore  radius),  the  variation 
of  Bv,  Nv,  Pv  and  Gy  is  shown  in  Figure  3.21.  In  this  case  as  in  the 
previous  one  (32  A pore  radius),  both  Bv  and  Gy  increase  as  the 
densification  temperature  increases  from  180  °C  (453  K).  The  genus 
goes  to  zero  at  a lower  temperature  than  for  the  32  A pore  radius 
samples. 
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Figure  3.20  Variation  of  the  number  of  nodes  (Ny),  number  of 
branches  (By),  number  of  separate  parts  (Py)  and  genus  (Gy)  as  a 
function  of  densification  temperature  for  samples  type  B (Model  12). 
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Figure  3.21  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (By),  number  of  separate  parts  (Py)  and  genus  (Gy)  as  a 
function  of  densification  temperature  for  samples  type  C (Model  12). 
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3.4  Topological  Indexing 


3.4.1  Relative  Genus  y 


In  order  to  make  it  easier  to  compare  the  structural  changes  at 
different  stages  of  densification,  the  concept  of  relative  genus  (y)  is 
introduced,  being  defined  as  follows: 


y = 


Gv 
Gv°  ’ 


(3.39) 


where  Gv°  is  the  genus  of  the  initial  as-dried  state,  and  Gy  is  the  genus 
at  any  given  stage  of  densification.  The  usefulness  of  this  topological 
index  is  that  it  permits  comparison  of  the  state  of  interconnection  at 
different  stages  of  densification,  or  even  of  different  samples,  without 
requiring  the  use  of  large  numbers. 

In  the  dried  stage  Gy  = Gy°  and  y = 1 . When  the  genus  goes  to  zero 
the  relative  genus  y = 0,  characterizing  the  third  stage  of  densification. 

Figure  3.22  shows  the  variation  of  relative  genus  (y)  with 
densification  temperature  for  four  sets  of  samples  of  12  A average  pore 
radius.  The  values  of  y were  obtained  by  applying  Model  12 
(cylindrical  model  with  variable  Py).  The  relative  genus  of  sets  of 
samples  of  12  A,  32  A and  81  A average  pore  radii  as  a function  of 
densification  temperature  is  shown  in  Figure  3.23.  The  values  of  y used 
in  this  figure  were  obtained  by  applying  the  Model  12  (cylindrical 
model  with  variable  Py).  The  variation  of  the  genus  of  the  pore 
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Figure  3.22  Variation  of  relative  genus  (y)  as  a function  of 
densification  temperature  for  samples  type  A. 
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Figure  3.23  Variation  of  relative  genus  (y)  as  a function  of 
densification  temperature  for  silica  gel  samples  with 
12  A,  32  A and  81  A average  pore  radii. 
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network  of  silica  gel  monoliths  as  a function  of  densification 
temperature  is  shown  in  Figure  3.24. 

The  variation  of  relative  genus  as  a function  of  volume  fraction  of 
pores  is  shown  in  Figures  3.25  to  3.28  for  samples  type  A (12  A 
average  pore  radius).  Figures  3.25  to  3.28  suggest  that  the  pore 
connectivity  goes  to  zero  when  the  volume  fraction  of  pores  is  around 
0.4.  This  is  a relatively  high  value,  suggesting  that  an  initial  pore 
coordination  number  (CN)  greater  than  4 would  be  more  appropriate. 
This  will  be  discussed  further  in  section  3.7.  Figure  3.29  shows  the 
variation  of  relative  genus  as  a function  of  volume  fraction  of  pores  for 
samples  type  B (32  A average  pore  radius). 

3.4.2  Topological  Index  B 

Because  the  number  of  nodes  is  constant  during  the  second  stage  of 
densification,  and  the  genus  is  constant  during  the  third  stage,  the 
number  of  branches  seems  to  be  a convenient  parameter  to  use  to  follow 
the  evolution  of  the  structure.  The  number  of  branches  (Bv)  varies 
from  an  initial  state  corresponding  to  the  dried  structure  (Bv°)  until  it  is 
zero  when  the  densification  is  completed.  Thus,  the  number  of  nodes 
(Nv)  does  not  reflect  structural  changes  in  the  sample  when  the  plateau 
of  constant  Nv  is  present.  Also,  the  genus  (Gv)  is  insensitive  to 
structural  transformations  during  the  third  stage  of  densification 
because  it  is  zero  during  the  third  stage,  when  changes  in  Bv,  Nv  and  Pv 
still  occur. 
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Figure  3.24  Variation  of  the  genus  (Gy)  as  a function  of  densification 
temperature  for  silica  gel  samples  with  12  A,  32  A and 
81  A average  pore  radii. 
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Figure  3.25  Variation  of  relative  genus  (7)  as  a function  of  volume 
fraction  of  pores  (Vy)  for  samples  type  Al. 
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Figure  3.26  Variation  of  relative  genus  (y)  as  a function  of  volume 
fraction  of  pores  (Vv)  for  samples  type  A2. 
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Figure  3.27  Variation  of  relative  genus  (7)  as  a function  of  volume 
fraction  of  pores  (Vv)  for  samples  type  A3. 
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Figure  3.28  Variation  of  relative  genus  (y)  as  a function  of  volume 
fraction  of  pores  (Vv)  for  samples  type  A4. 
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Figure  3.29  Variation  of  relative  genus  (y)  as  a function  of  volume 
fraction  of  pores  (Vv)  for  samples  type  B. 
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Based  on  the  variation  of  the  number  of  branches  (Bv)  a topological 
index  6 (fraction  of  removed  branches)  is  defined  as  follows: 

Bv 

B = l-B^’  (3-4°) 

where  Bv  is  the  number  of  branches  at  any  stage  of  densification  and 
Bv°  is  the  number  of  branches  of  the  initial  dried  stage,  which  is  used  as 
reference.  Therefore,  in  the  dried  stage  By  = Bv°  and  6 = 0,  meaning 
the  initial  point  of  densification.  When  the  sample  becomes  fully  dense, 
the  number  of  branches  in  the  structure  is  zero  and  then  6=1, 
characterizing  a fully  dense  sample.  Thus,  if  the  dried  stage  corresponds 
to  a maximum  number  of  branches,  the  topological  index  6 will  vary 
from  0 to  1 as  the  densification  proceeds. 

Again,  the  usefulness  of  such  a type  of  topological  index  is  that  it 
makes  it  easier  to  compare  different  stages  of  densification  for  a given 
type  of  sample  or  between  different  samples  by  avoiding  the  use  of 
large  numbers.  Using  this  index  one  can  describe  the  structural 
evolution  of  a material  undergoing  sintering  in  a very  simple  and  direct 
manner. 

The  topological  index  6 is  shown  in  Figure  3.30  as  a function  of 
densification  temperature  for  four  sets  of  samples  of  12  A average  pore 
radius.  Model  12  was  used  to  obtain  the  values  of  number  of  branches 
used  in  this  section.  It  can  be  seen  from  Figure  3.30  that  6 varies 
relatively  slowly  from  180  °C  (453  K)  to  temperatures  in  the  range  800 
°C  (1073  K)  to  900  °C  (1173  K)  above  which  6 increases  towards  1 
very  fast.  The  variation  of  topological  index  6 with  densification 
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Figure  3.30  Variation  of  the  topological  index  6 as  a function  of 
densification  temperature  for  samples  type  A. 
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temperature  for  samples  of  12  A,  32  A and  81  A average  pore  radius  is 
shown  in  Figure  3.31.  If  can  be  seen  from  Figure  3.19  the  different 
behavior  during  densification  of  the  gel-silicas  of  larger  pore  radius  as 
compared  to  the  samples  of  12  A pore  radius.  The  variation  of  the 
number  of  branches  (Bv)  as  a function  of  densification  temperature  for 
silica  gel  monoliths  with  12  A,  32  A and  81  A average  pore  radii  is 
shown  in  Figure  3.32. 

The  variation  of  the  topological  index  6 as  a function  of  the  volume 
fraction  of  pores  for  silica  gel  monoliths  with  12  A average  pore  radius 
is  shown  in  Figures  3.33  to  3.36.  In  these  figures  Model  12  (cylindrical 
with  variable  Pv)  was  used.  Using  Model  10  (tetrahedral  with  variable 
Pv)  one  obtains  Figure  3.37,  for  the  variation  of  8 versus  Vv  of  pores. 

3.5  Comparison  of  Models 

Depending  on  the  shape  assumed  for  the  geometric  model,  it  will 
contain  a certain  volume  and  present  a certain  surface  area  which  will 
vary  from  shape  to  shape.  Since  the  calculation  of  topological 
parameters  depends  on  the  volume  and  surface  area  of  each  geometrical 
unit,  it  is  expected  that  different  shapes  will  give  rise  to  differences  in 
the  values  of  Bv,  Nv,  Gy  and  Py  estimated.  However,  any  difference 
should  be  present  in  such  a way  that  they  do  not  compromise  the 
usefulness  and  description  of  the  structural  evolution.  Creating  several 
different  models  with  different  assumptions  covers  the  expected  range 
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Figure  3.31  Variation  of  the  topological  index  B as  a function  of 
densification  temperature  for  silica  gel  samples  with  12  A,  32  A 
and  81  A average  pore  radii. 
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Figure  3.32  Variation  of  the  number  of  branches  (By)  as  a function  of 
densification  temperature  for  silica  gel  samples  with  12  A,  32  A 
and  81  A average  pore  radii. 
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Figure  3.33  Variation  of  volume  fraction  of  pores  (Vv)  as  a function 
of  the  topological  index  6 for  samples  type  A1  (Model  12). 
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Figure  3.34  Variation  of  volume  fraction  of  pores  (Vy)  as  a function 
of  the  topological  index  6 for  samples  type  A2  (Model  12). 
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Figure  3.35  Variation  of  volume  fraction  of  pores  (Vv)  as  a function 
of  the  topological  index  6 for  samples  type  A3  (Model  12). 
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Figure  3.36  Variation  of  volume  fraction  of  pores  (Vv)  as  a function 
of  the  topological  index  6 for  samples  type  A4  (Model  12). 
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Figure  3.37  Variation  of  volume  fraction  of  pores  (Vv)  as  a function 
of  the  topological  index  8 for  samples  type  A2  (Model  10). 
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of  variations  for  the  topological  parameters  during  densification.  When 
a certain  path  of  evolution  is  recognized  using  different  models  and 
assumptions,  and  when  the  same  ranges  of  numbers  are  obtained  using 
different  conditions,  one  can  be  more  confident  about  the  description  of 
the  structural  transformations  based  on  the  proposed  models. 

The  variation  of  number  of  nodes  (Nv),  number  of  branches  (Bv), 
and  genus  (Gy)  with  densification  temperature  is  shown  in  Figure  3.38, 
for  tetrahedral  (Model  3)  and  cylindrical  (Model  6)  models.  The 
variation  in  the  geometric  shape  (triangular  or  circular  pores)  seems  to 
have  a smaller  effect  on  the  values  of  the  topological  parameters  than 
the  assumption  of  either  L (Model  3)  or  D (Model  4)  as  seen  in  Figure 
3.39.  The  path  of  structural  evolution  in  both  cases  (Figure  3.38  and 
Figure  3.39)  is  about  the  same,  independent  of  the  model  assumed. 

Figure  3.40  shows  the  variation  of  Bv,  Nv,  Pv  and  Gv  for  a set  of  12 
A pore  radius  samples  as  a function  of  temperature.  Model  6 shown  in 
the  figure  assumes  Pv  = 1 and  Model  12  accepts  a variable  Pv.  As 
expected  the  difference  between  those  two  models  (6  and  12)  occurs 
only  in  the  third  stage  of  densification.  That  means  that  the  very  simple 
models  can  be  used  to  describe  structural  changes  during  most  of  the 
process. 

A general  picture  of  the  results  obtained  by  various  different  models 
is  shown  in  Figure  3.41.  The  number  of  pores  (Nv)  for  the  spherical 
model  is  shown  as  open  circles  (legend  1 Nv).  The  number  of  branches 
(Bv)  for  tetrahedral  (3,  4 and  10)  and  cylindrical  (6,  7 and  12)  models 
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Figure  3.38  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (Bv)  and  genus  (Gy)  as  a function  of  densification 
temperature  for  samples  type  A2.  Model  3 is  a tetrahedral 
model  and  Model  6 is  a cylindrical  model. 
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Figure  3.39  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (By)  and  genus  (Gy)  as  a function  of  densification 
temperature  for  samples  type  A2.  Model  3 assumes  L 
and  Model  4 assumes  D. 
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Figure  3.40  Variation  of  the  number  of  nodes  (Nv),  number  of 
branches  (Bv)  and  genus  (Gy)  as  a function  of  densification 
temperature  for  samples  type  A2.  Model  6 assumes  Pv  = 1 and 
Model  12  accepts  a variable  Pv. 
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Figure  3.41  Variation  of  the  number  of  nodes  (Nv)  and  number  of 
branches  (Bv)  as  a function  of  densification  temperature  for  samples 
type  A2.  The  numbers  in  the  legend  correspond  to  the  model  number. 
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are  also  shown  in  Figure  3.41,  with  the  numbers  in  the  legend 
corresponding  to  the  model  number. 

As  pointed  out  before,  the  spherical  model  seems  to  provide  a 
(lower),  limiting  case  for  the  number  of  pores,  even  though  it  is  not  a 
realistic  model.  For  a given  shape,  such  as  cylindrical,  models  with  a 
variable  Pv  (Model  12)  differ  from  models  with  Pv  = 1 (Model  6)  only 
during  the  third  stage  of  densification.  The  cylindrical  Model  7,  which 
assumes  a constant  value  for  Nv  during  the  entire  process,  differs  from 
Model  12  only  in  the  last  stage  of  densification.  For  the  tetrahedral 
model,  assumptions  of  D = 4Vv/Sv,  as  opposed  to  assuming  L = 4(1  - 
Vv)/Sv,  lead  to  relatively  large  differences,  even  though  the  shape  of  the 
curves  are  very  similar. 

3.6  Transmission  Electron  Microscopy  Analysis 

A gel-silica  sample  of  12  A average  pore  radius  in  its  as  dried 
condition  was  analyzed  by  high-resolution  transmission  electron 
microscopy  (TEM)  by  Dr.  Carrole  Perry,  Oxford  University  and 
Brunei  Universities,  London.  The  result  is  shown  in  Figure  3.42.  The 
picture  shows  the  pore  structure  of  the  dried  gel-silica.  As  seen  in  the 
picture  the  scale  of  pore  sizes  is  the  same  scale  as  obtained  by  nitrogen 
gas  adsorption  and  the  scale  of  the  interpore  distance  is  the  same  as  the 
scale  of  interpore  distance  obtained  from  the  topological  models 
proposed.  This  TEM  micrograph  therefore  gives  support  to  the 
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Figure  3.42  High-resolution  transmission  electron  microscopy  photo  of 
a dried  silica  gel  of  12  A average  pore  radius.  In  the  photo  1 cm 
represents  8 nm.  (Photo  courtesy  of  Dr.  Carrole  Perry,  Oxford  and 
Brunei  Universities,  England). 
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proposed  geometric  models,  by  showing  that  the  scale  of  the  structure  as 
obtained  by  using  the  expressions  for  the  topological  models  is  realistic. 

3.7  Sensitivity  Analysis  of  the  Pore  Coordination  Number 

A basic  assumption  of  all  the  topological  models  proposed  is  that  in 
its  initial  (dried)  stage  the  silica  gel  samples  present  an  average  pore 
coordination  number  (CN)  of  four.  The  purpose  of  this  section  is  to 
analyze  the  influence  of  CN  on  the  determination  of  topological 
parameters,  and  to  provide  a tool  so  that  one  can  change  and  adapt  the 
models  if  an  initial  CN  different  than  four  is  expected  to  yield  more 
realistic  results  for  any  given  application. 

A wide  range  of  pore  coordination  numbers  (CN)  is  shown  in  Figure 
3.43.  In  this  figure  the  values  of  CN  were  calculated  using  the  Model  12 
(cylindrical,  with  variable  Pv).  Initial  pore  coordination  numbers 
(CN°)  varying  from  3 to  24  were  chosen.  Figure  3.43  shows  those 
initial  CN's  at  180  °C  (453  K)  and  as  the  densification  temperature  goes 
up,  the  values  of  CN  decrease,  until  they  eventually  become  zero. 

The  variation  of  the  number  of  branches  (Bv)  with  densification 
temperature  is  shown  in  Figure  3.44.  It  can  be  seen  that  the  choice  of 
the  initial  pore  coordination  number  does  not  affect  the  values  of  Bv, 
except  in  the  third  stage  of  densification.  Again,  this  provides  one  more 
argument  for  the  convenience  of  using  the  number  of  branches  (By)  and 
the  topological  index  6 to  describe  the  evolution  of  the  structure  during 
densification.  It  should  be  pointed  out  that  the  curve  representing  the 
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Figure  3.43  Variation  of  the  pore  coordination  number  (CN)  as  a 
function  of  densification  temperature  for  initial  CN's  varying 
from  3 to  24  (Samples  type  A2). 
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Figure  3.44  Variation  of  the  number  of  branches  (Bv)  as  a function  of 
densification  temperature  for  samples  type  A.  The  numbers  in  the 
legend  represent  the  pore  coordination  numbers  in  the 
beginning  of  densification  (Samples  type  A2). 
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variation  of  Bv  with  temperature  is  practically  the  same  for  CN's 
varying  from  4 to  8. 

The  variation  in  the  number  of  branches  Bv  as  a function  of  the 
volume  fraction  of  pores  Vy  is  shown  in  Figure  3.45.  Again,  for  CN°'s 
between  4 and  8 there  is  no  variation  in  the  path  of  structural  change. 

The  variation  of  genus  (Gv)  as  a function  of  densification 
temperature  is  shown  in  Figure  3.46,  when  CN°  varies  from  3 to  24. 
Looking  at  the  shape  of  the  curves,  particularly  to  the  extrapolation  of 
the  last  segment  of  each  curve  and  the  temperature  obtained  by  such 
extrapolation  it  looks  like  CNo's  of  8 to  24  would  require  sharp  changes 
in  their  respective  curves  (of  variation  of  Gy  with  temperature)  in 
order  to  obtain  zero  Gv  at  temperatures  in  the  range  of  1100  °C  (1373 
K).  This  suggests  that  CNo’s  smaller  than  8 would  provide  more 
realistic  results  when  used  in  the  models. 

The  variation  of  genus  (Gv)  with  volume  fraction  of  pores  (Vy)  is 
seen  in  Figure  3.47.  In  this  figure,  if  one  does  not  consider  the  last 
point  of  each  curve  (when  Vv  is  minimal  and  Gv  is  equal  to  zero),  and 
extrapolate  each  curve,  the  range  of  Vv  for  which  Gv  = 0 will  vary 
approximately  between  Vy  = 0.3  to  Vy  = 0.4.  Thus  for  the  wide  range 
of  initial  CN°’s,  pore  closure  should  occur  in  that  range  of  volume 
fraction  of  pores. 
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Figure  3.45  Variation  of  the  number  of  branches  (Bv)  as  a function  of 
the  volume  fraction  of  pores  (Vy)  for  assumed  initial  pore  coordination 
numbers  varying  from  3 to  24  (Samples  type  A2). 
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Figure  3.46  Variation  of  the  genus  (Gy)  as  a function  of  densification 
temperature  when  the  assumed  initial  pore  coordination  number 
varies  from  3 to  24  (Samples  type  A2). 
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Figure  3.47  Variation  of  the  genus  (Gy)  as  a function  of  volume 
fraction  of  pores  (V y)  when  the  assumed  initial  pore  coordination 
number  varies  from  3 to  24  (Samples  type  A2). 
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3.8  Topological  Analysis  of  the  Time  and  Humidity  Effects 

The  variation  of  genus  (Gv)  with  densification  temperature  for  12  A 
average  pore  radius  samples  is  shown  in  Figure  3.48.  A set  of  samples 
was  kept  at  the  (maximum)  densification  temperature  for  a period  of  1 
h (standard  procedure)  while  a second  set  was  kept  at  densification 
temperature  for  10  h.  The  figure  suggests  that  the  longer  the 
densification  time  the  more  the  genus  is  going  to  decrease.  A similar 
type  of  behavior  can  be  seen  in  Figure  3.49  where  the  topological  index 
6 is  plotted  as  a function  of  the  (maximum)  time  at  the  densification 
temperature.  Figure  3.49  suggests  that  the  longer  the  time  used  in 
densification  at  a given  temperature,  a larger  fraction  of  branches  will 
be  removed  from  the  structure. 

In  order  to  evaluate  further  the  effect  of  atmosphere  on  the  sintering 
behavior  of  gel-silica  monoliths,  a set  of  samples  (Al)  was  densified  in  a 
dry-air  (DA)  atmosphere  and  another  set  in  a humid  (OF)  atmosphere. 
The  variation  of  genus  with  densification  temperature  for  those  samples 
is  shown  in  Figure  3.50.  The  figure  suggests  that  samples  densified  in 
humid  atmosphere  close  their  pores  faster  than  samples  processed  in  a 
dry-air  atmosphere.  This  tendency  could  explain  why  samples  sintered 
in  a box-furnace  (open  to  ambient  atmosphere,  with  no  flow  of  dry-air) 
foam  or  bloat  much  more  frequently  and  at  a lower  temperature  than 
samples  densified  in  a closed  furnace  with  a dry-air  flow. 
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Figure  3.48  Variation  of  the  genus  (Gy)  as  a function  of 
densification  temperature  for  samples  type  A1  kept  at 
the  densification  temperature  for  1 h and  10  h. 
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Figure  3.49  Variation  of  the  topological  index  6 as  a function  of  the 
time  the  sample  was  held  at  800  °C  for  samples  type  A1  and  A2. 
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Figure  3.50  Variation  of  the  genus  (Gy)  as  a function  of  densification 
temperature  for  samples  type  A1  densified  in  a dry-air  (DA) 
atmosphere  and  in  a humid  (OF)  atmosphere. 


CHAPTER  4 

MECHANICAL  ANALYSIS 


4.1  The  Strength  of  Brittle  Materials 

Strength  of  materials  is  a science  that  deals  with  the  relation  between 
external  loads,  internal  forces,  and  the  deformation  caused  by  those 
forces  [399].  In  order  for  one  to  study  and  obtain  expressions  relating 
those  loads,  forces,  and  deformations,  the  following  assumptions  are 
usually  made  [399]: 

1)  The  body  is  continuous;  meaning  that  it  contains  no  voids  or 
empty  spaces. 

2)  The  body  is  homogeneous;  meaning  that  it  has  identical  properties 
at  all  points. 

3)  The  body  is  isotropic;  meaning  that  the  property  considered  does 
not  vary  with  direction  or  orientation. 

The  relationships  obtained  using  the  above  assumptions  form  the 
body  of  the  mechanics  of  materials,  or  mechanical  resistance  of 
materials. 

Another  approach  to  the  problem  of  achieving  a treatable 
(mathematical)  aspect  concerning  the  strength  of  materials  is  by 
considering  that  it  is  not  a continuous  body,  but  a discontinuous  body 
formed  by  individual  atoms.  This  is  an  atomistic  approach,  and  it 
relates  to  the  strength  of  a body  as  a whole  to  the  strength  of  the 
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individual  bonding  between  the  atoms  and/or  molecules  that  compose 
the  body.  In  this  case  it  is  also  assumed  that  the  body  is  homogeneous 
and  isotropic,  so  that  calculations  made  in  the  atomic  scale  can  be 
extrapolated  for  the  entire  body. 

Real  materials  have  both  components:  the  atomistic  contribution  and 
the  macroscopic  nature  of  the  body.  Both  types  of  approaches  can  be 
combined  in  order  to  yield  a more  complete  description  of  the 
mechanical  state  of  a material.  The  atomistic  approach  provides  the 
study  of  the  strength  of  brittle  materials  with  an  upper  limit  for  its 
mechanical  strength.  In  the  case  of  covalently  bonded  atoms,  the  force 
associated  with  atomic  interactions  can  be  represented  by  the  Morse 
potential  [400,4011 

U = U0  ( exp(-2B  (x-a))  - 2 exp(  -B  (x-a)))  , (4.1) 

where  U is  the  potential  energy  per  unit  of  surface  area  of  two 
fracture  surfaces; 

U0  and  B are  constants; 

x is  the  distance  between  the  two  fracture  surfaces; 

a is  the  interatomic  distance  of  the  atoms  being  separated. 

The  binding  force  per  unit  area  is  given  by  [399] 
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Assuming  that  Hook's  law  holds,  and  using  Eq.(4.2)  it  can  be  shown 
that  [2] 


where  or  is  the  theoretical  cohesive  strength; 

E is  the  Young's  modulus. 

Therefore,  the  cohesive  strength  of  a material  can  be  calculated  by 
knowing  the  Young's  modulus  and  U0,  which  represents  the  binding 
potential  at  the  equilibrium  atomic  separation  (when  x = a). 

For  a silica  glass,  the  silicon-oxygen  bond  represents  the  strongest 
cohesive  forces  in  the  network;  with  a bond  energy  of  106  kcal/mole 
[2].  Considering  a Young's  modulus  of  72  GPa  for  the  silica  glass,  the 
theoretical  cohesive  strength  is  18  GPa  [2].  Hillig  [402]  and  Proctor  et 
al.  [403]  report  maximum  experimental  values  of  fracture  strength  in 
the  range  13-15  GPa. 

Ordinary  glasses  usually  present  a much  lower  fracture  strength  (30- 
100  MPa)  because  of  the  presence  of  flaws  in  the  surface  of  the  glass 
[2].  The  flaws  cause  concentration  of  stresses,  causing  the  glass  to  fail  at 
an  applied  load  much  lower  than  predicted  theoretically. 

The  effect  of  surface  flaws  can  be  accounted  for  by  using  the 
Griffith  equation  [404,405] 


(4.3) 


(4.4) 
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where  Gf  is  the  stress  required  for  fracture; 
y is  the  surface  energy; 

C is  the  flaw  depth. 

For  an  elliptical  crack  of  length  2c  in  a wide  plate,  the  maximum 
stress  at  the  tip  of  the  crack  Gmax  is  given  the  Inglis  equation  [406] 


The  rate  of  the  crack  propagation  in  a solid  (vc)  can  expressed  as 
[407,408] 


where  vj  is  the  longitudinal  speed  of  sound  in  the  solid; 
c0  is  half  the  critical  length  of  Griffith  flaw; 
c is  half  the  crack  length  at  a given  time; 
k is  a dimension  less  constant. 

Equation  (4.6)  shows  that  the  speed  of  crack  propagation  increases 
when  the  length  of  the  crack  increases.  Reed-Hill  [409]  points  out  that 
this  dependence  of  crack  speed  with  the  speed  of  sound  occurs  because 
the  displacement  of  material  adjacent  to  the  crack  can  only  occur  at  the 
speed  of  sound. 


(4.5) 


where  o is  the  average  tensile  stress; 

Pt  is  the  radius  of  curvature  of  the  crack  tip. 


(4.6) 
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Experimental  ways  of  studying  the  effect  of  cracks  in  the  strength  of 
brittle  materials  and  the  validity  of  theoretical  expressions  correlating 
the  presence  of  flaws  with  strength  usually  involve  deliberately 
introducing  flaws  in  the  material  [407-413].  By  introducing  a flaw,  it  is 
possible  to  evaluate  the  effects  on  the  strength  of  its  size,  shape,  and 
external  media  [414-419].  The  radial  cracking  that  can  occur  around 
Vickers  microhardness  indentations  seem  to  be  particularly  useful  in 
analyzing  crack  growth  characteristics  [420-425]  and  toughness  [426- 
431]. 

Flaws  can  originate  from  a large  number  of  sources  including 
processing,  handling,  cutting,  packaging,  transporting,  etc.,  originating 
materials  with  a wide  variety  of  flaw  lengths  and  a corresponding  wide 
variety  of  strengths  [432,433]. 

As  the  size  of  the  testing  specimen  increases,  the  chances  that  the 
sample  contains  a flaw  also  increases,  meaning  that  larger  samples  are 
likely  to  have  lower  strength  than  smaller  samples.  Besides  the  effect  of 
size,  Rudnick  et  al.  [434]  points  out  that  the  fracture  strength  is 
characterized  by  significant  scatter  of  data,  and  a stress  distribution 
effect;  e.g.,  strength  values  obtained  in  bending,  for  instance,  exceed 
those  obtained  from  uniform  tension  tests. 

In  a porous  material  (like  silica  gels  densified  at  intermediate  ranges 
of  temperature)  the  presence  of  voids  can  influence  the  strength  of  the 
material.  Under  the  application  of  stress  the  voids  present  in  the 
structure  can  enlarge,  and  new  voids  can  be  formed.  Together,  the 
creation  of  new  voids  and  enlargement  of  old  ones  can  cause  the 
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macroscopic  effect  known  as  dilatancy  [435].  However,  the  effects  of 
dilatancy  seems  to  be  more  pronounced  for  relatively  low  volume 
fraction  of  pores  (~  few  percent),  and  also  the  phenomena  is  usually 
related  to  a non-elastic  increase  in  volume  [435]. 

The  effect  of  porosity  in  the  strength  of  silica-gel  monoliths  are 
likely  to  have  at  least  two  aspects:  (1)  the  porosity  decreases  the  actual 
cross-section  of  material  which  supports  the  load;  and  (2)  pores  can  act 
as  stress  risers. 

Because  of  the  characteristics  of  brittle  fracture  (rapid  failure  and  no 
macroscopic  plastic  deformation),  the  statistics  of  the  data  is  a topic  of 
special  importance.  Usually  the  strength  data  are  treated  as  having  a 
Gaussian  (normal)  distribution  [399],  or  a Weibull  distribution 
[399,436-441]: 


where  P(x)  is  the  cumulative  frequency  distribution  of  random 
variable  x; 

Xq  is  the  minimum  allowable  value  of  x; 

0 is  a scale  parameter; 
m is  a shape  parameter  (Weibull  modulus). 

For  design  applications  and  for  dealing  with  probabilities  of  failure 
the  Weibull  distribution  is  usually  preferred,  because  it  does  not  have 
the  unrealistic  long  tails  of  the  Gaussian  distribution  [399]. 


(4.7) 
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4.2  Microhardness 


The  Vickers  hardness  number  (VHN),  also  called  the  diamond- 
pyramid  hardness  test  (DPH),  is  defined  as  the  load  divided  by  the 
surface  area  of  the  indentation  [399].  The  VHN  can  be  determined  by 
the  following  equation  [399]: 


VHN 


1.854  P 


(4.8) 


where  P is  the  applied  load  (kg); 

L is  the  average  length  of  the  diagonals  of  the  indentation  (mm). 

Microhardness  measurements  were  made  using  a semi-automatic 
LECO  DM-400F  Hardness  Tester.  A load  of  0.1  kg  was  used  for  the  12 
A and  32  A pore  radius  samples  and  0.025  kg  was  used  for  the  81  A 
pore  radius  samples.  A loading  time  of  15  s was  selected  for  all  testing. 
The  variation  of  Vickers  hardness  number  (VHN)  with  densification 
temperature  for  sets  of  samples  with  12  A,  32  A and  81  A average  pore 
radius  is  shown  in  Figure  4.1.  It  can  be  seen  from  the  figure  that  for  all 
sets  of  samples  the  VHN  increases  with  increasing  sintering 
temperatures  and  the  rate  of  increase  is  sharper  after  a certain  critical 
temperature  which  increases  for  larger  pore  sizes.  As  pointed  out 
previously,  for  all  the  figures  in  this  work,  the  error  bars  represent  the 
standard  error  of  the  mean;  the  curve  fitting  is  either  a polynomial  or 
just  a curve  that  connect  the  average  points,  not  implying  in  any  trends 
in  the  data  when  a peak  or  valley  occurs. 
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Figure  4.1  Variation  of  the  Vickers  hardness  number  (VHN)  as  a 
function  of  densification  temperature  for  silica  gel  samples  of 
12  A,  32  A and  81  A average  pore  radii. 
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The  initial  VHN  for  the  dried  (180  °C,  453  K)  silica  gel  samples  of 
12  A average  pore  radius  is  higher  than  the  VHN  for  32  A average  pore 
radius  samples,  which  is,  by  its  turn,  higher  than  the  VHN  for  81  A 
pore  radius  samples.  Considering  the  VHN  versus  temperature  curve 
for  12  A pore  radius  samples  it  can  be  seen  that  there  is  a continuous 
increase  in  hardness  as  the  densification  temperature  increase.  When 
the  sintering  temperature  increases  above  900  °C  (1173  K)  the  12  A 
pore  radius  structure  shows  a sharp  increase  in  VHN.  The  hardness  of  a 
12  A pore  radius  sample  densified  at  1150  °C  (1423  K)  is  770  +/-  10 
kg/mm2,  which  is  comparable  to  the  microhardness  of  fused  silica  [442]. 

The  variation  of  VHN  with  densification  temperature  for  32  A pore 
radius  samples  (Figure  4.1)  shows  smaller  values  of  hardness  than  for 
12  A pore  radius  samples  sintered  at  same  temperatures.  The  increase 
in  VHN  with  temperature  is  relatively  slow  up  to  temperatures  in  the 
range  of  1100  °C  (1373  K)  to  1150  °C  (1423  K)  beyond  which  the 
hardness  increases  sharply  to  a value  of  750  +/-  12  kg/mm2  at  125  °C 
(1523  K).  Silica  gels  of  81  A pore  radius  have  a similar  behavior 
compared  to  32  A pore  radius  samples  and  even  lower  values  of  VHN 
for  corresponding  densification  temperatures.  The  observed  hardness 
of  a sample  sintered  at  1250  °C  (1523  K)  is  627  +/-  10  kg/mm2. 

4.3  Flexural  Testing 

Four  points  bending  tests  were  performed  in  a MTS  810/442  testing 
machine.  Example  of  the  type  of  sample  used  in  the  flexural  tests  is 
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shown  in  Figure  4.2.  All  the  tests  were  performed  using  a cross  head 
speed  of  0.016  mm/min  in  ambient  air. 

A typical  plot  of  load  versus  displacement  is  shown  in  Figure  4.3. 
From  the  bending  curves  one  obtains  the  values  of  the  load  of  fracture 
(P)  which  is  used  in  the  following  expression  to  obtain  the  flexural 
strength  (ab)  [443] 


Ob 


3 P (L  - 1) 

2 w t2  ’ 


(4.9) 


where  L is  the  distance  between  the  lower  supporting  points; 

1 is  the  distance  between  the  upper  loading  points; 
w is  the  width  of  the  tested  sample; 
t is  the  thickness  of  the  tested  sample. 

The  variation  of  flexural  strength  with  densification  temperature  is 
shown  in  Figure  4.4,  for  samples  of  12  A average  pore  radius.  This 
figure  shows  a continuous  increase  in  flexural  strength  as  the 
densification  temperature  increases.  The  increase  in  strength  is  more 
pronounced  for  temperatures  above  700  °C  (973  K). 

The  flexural  modulus  (Eb)  of  the  samples  was  calculated  using  the 
following  expression  (see  Appendix  for  the  derivation  of  this 
expression): 

7 P L3 
32  A w t3  ’ 


Eb  = 


(4.10) 
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Figure  4.2  Silica  gel  monoliths  prepared  for  the  tests  indicated. 
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Figure  4.3  Load  (P)  versus  displacement  for  a four  point 
bending  test  of  a dried  silica  gel. 
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Figure  4.4  Variation  of  flexural  strength  as  a function  of  densification 
temperature  for  silica  gel  samples  with  12  A average  pore  radius. 


where  A is  the  maximum  displacement  measured  in  the  center  of  the 
specimen. 

The  variation  of  flexural  modulus  with  densification  temperature  for 
12  A average  pore  radius  samples  is  shown  in  Figure  4.5.  The  flexural 
modulus  increases  continuously  with  increasing  sintering  temperature. 
The  increase  in  flexural  modulus  occurs  more  rapidly  for  temperatures 
above  900  °C  (1173  K). 

4.4  Uniaxial  Compression  Testing 

For  plastic  solids  the  strength  measured  in  uniaxial  compression  is 
the  same  as  the  strength  measured  in  tension.  In  the  case  of  brittle 
solids  (e.g.,  glasses),  the  compressive  strength  is  roughly  15  times 
larger  [444].  According  to  Ashby  and  Jones  [444,  p.167]  the  reason  for 
this  is  that 

cracks  in  uniaxial  compression  propagate  stably  (Figure  4.6)  and 
twist  out  of  their  original  orientation  to  propagate  parallel  to  the 
compression  axis.  In  this  case,  fracture  is  not  caused  by  the  rapid 
unstable  propagation  of  one  crack,  but  the  slow  extension  of  many 
cracks  to  form  a crushed  zone.  It  is  not  the  size  of  the  largest  crack 
that  counts,  but  that  of  the  average  size  of  cracks  (aavg). 

The  uniaxial  compressive  strength  (ac)  can  be  correlated  to  the 

toughness  (Kic)  as  follows  [444]: 


FLEXURAL  MODULUS  (GPa) 
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Figure  4.5  Variation  of  the  flexural  modulus  as  a function  of 
densification  temperature  for  silica  gel  samples  of 
12  A average  pore  radius. 
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Figure  4.6  Propagation  of  cracks  in  uniaxial  compression  [444]. 
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where  ku  is  a constant  (~  15). 

A detailed  analysis  on  the  effect  of  pores  and  small  cracks  on  the 
strength  of  brittle  solids  under  compression  is  given  by  Ashby  and 
Hallam  [445]  and  Sammis  and  Ashby  [446]. 

When  a static  uniaxial  compression  test  is  performed  during  a short 
period  of  time  and  at  relatively  low  loading  rates  (e.g.,  0.07  MPa/s  to 
0.7  MPa/s)  the  test  is  referred  to  as  time-independent,  as  contrasted  to 
time-dependent  (creep)  tests  [447]. 

Two  types  of  compression  samples  were  prepared,  as  shown  in 
Figure  4.2.  Those  samples  were  densified  in  the  range  of  180  °C  (453 
K)  to  1000  °C  (1273  K)  in  dry-air  atmosphere,  and  at  1150  °C  (1423  K) 
in  a chlorine  atmosphere.  The  specimens  were  tested  under  uniaxial 
compression  in  a MTS  810/458  at  a cross  head  speed  of  0.02  mm/min. 
The  displacement  was  directly  measured  using  a clip  on  gage 
extensometer. 

A typical  result  of  a uniaxial  compression  test  is  shown  in  Figure 
4.7.  The  values  of  load  (P)  are  used  to  obtain  the  compression  strength 
(ac)  as  follows: 

<*c  = X'  (4.12) 


where  A is  the  cross  sectional  area  of  the  tested  specimen. 

The  variation  of  compression  strength  with  densification  temperature 
for  a set  of  samples  of  12  A average  pore  radius  is  shown  in  Figure  4.8. 
The  figure  shows  an  increase  in  compression  strength  as  the 
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Figure  4.7  Variation  of  the  uniaxial  compression  load  (P)  versus 
displacement  for  a dried  silica  gel.  The  straight  line  represents 
the  region  in  the  plot  which  was  used  to  calculate  the  modulus. 
The  initial  curved  region  is  associated  with  the 
deformation  of  the  Teflon®  lubricant. 
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Figure  4.8  Variation  of  the  compression  strength  as  a function  of 
densification  temperature  for  samples  type  A1  and  C. 
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densification  temperature  increases.  As  shown  in  Figure  4.8,  the 
compression  strength  for  a set  of  dried  silica  gels  with  81  A average 
pore  radius  (samples  type  C)  is  lower  than  the  strength  of  silica  gels 
with  12  A average  pore  radius  (samples  type  Al). 

The  compression  modulus  (Ec)  was  calculated  using  the  following 
expression: 


where  A ac  is  the  variation  in  stress  and  Ae  the  variation  in  strain 
during  the  test,  measured  in  the  elastic  region  of  the  curve  stress  versus 
strain  (not  shown).  Figure  4.9  shows  the  variation  in  compression 
modulus  with  densification  temperature  for  two  sets  of  samples  of  12  A 
average  pore  radius,  as  well  as  for  a set  of  dried  samples  of  81  A pore 
radius,  and  a set  of  fused  quartz  samples. 

The  effect  of  densification  atmosphere  in  the  compression  strength  is 
shown  in  Figure  4.10.  Samples  sintered  in  a dry-air  (DA)  atmosphere 
seems  to  have  a higher  compression  strength  than  samples  densified  in  a 
humid  atmosphere  (OF).  The  effect  of  atmosphere  on  the  compression 
modulus,  however,  seems  to  be  smaller,  as  shown  in  Figure  4.11.  In 
this  figure  it  is  seen  that  both  the  compression  modulus  of  samples 
processed  in  a dry-air  and  in  a humid  (ambient)  atmosphere  are  very 
close  for  every  tested  temperature. 


COMPRESSION  MODULUS 
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Figure  4.9  Variation  of  the  compression  modulus  as  a function  of 
densification  temperature  for  samples  type  Al,  A3  and  C. 

FQ  represents  fused  quartz  samples. 
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Figure  4.10  Variation  of  compression  strength  as  a function  of 
densification  temperature  for  samples  type  A1  densified  in  a 
dry-air  (DA)  atmosphere  and  in  a humid  (OF)  atmosphere. 
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Figure  4.11  Variation  of  compression  modulus  as  a function  of 
densification  temperature  for  samples  type  A1  densified  in  a dry-air 
(DA)  atmosphere  and  in  a humid  (OF)  atmosphere. 
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4.5  Brazilian  Testing 

The  Brazilian  test  (also  known  as  the  diametral-compression  test) 
was  developed  in  the  early  1950s  by  Cameiro  and  Barcellos  [448]  to 
evaluate  the  strength  of  concrete.  Because  of  both  the  simplicity  of  the 
test  and  ease  of  sample  preparation,  and  the  fact  that  this  test  yields  an 
indirect  measurement  of  tensile  strength  (otherwise  of  difficult 
experimental  determination  for  brittle  materials)  the  Brazilian  test  has 
been  used  ever  since  for  a variety  of  materials  [449-455]. 

The  stress  distribution  across  a loaded  diameter  of  a cylindrical 
sample  is  shown  in  Figure  4.12  [452].  The  maximum  tensile  stress  (at) 
which  acts  across  the  diameter  is  given  by  [452] 

at  = JiDH’  <4-14) 

where  P is  the  applied  load; 

D is  the  diameter  of  the  specimen; 

H is  the  length  of  the  specimen. 

Figure  4.12  shows  the  load  being  applied  over  a region  of  1/10  of 
the  diameter  of  the  cylinder.  If  the  load  arc  is  increased  to  20°  the 
reduction  in  maximum  tensile  stress  will  be  only  4%  compared  to  the 
ideal  situation  when  the  load  is  applied  on  a line  [456].  Therefore, 
Eq.(4.14)  can  be  used  to  determine  the  tensile  strength  of  cylindrical 
samples  even  in  the  presence  of  a small  contact  area  between  the  sample 
and  the  loading  fixtures. 
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Figure  4.12  Stress  distribution  across  a loaded  diameter  of 
a cylindrical  sample  in  the  Brazilian  test  [452]. 
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Two  sets  of  silica  gel  samples  of  12  A average  pore  radius  samples 
were  tested  for  the  Brazilian  test  in  a MTS  810/458  at  a cross  head 
speed  of  0.02  mm/min.  Teflon®  tape  was  used  as  lubricant  during  the 
tests,  in  order  to  prevent  failure  of  the  samples  at  the  contact  area  with 
the  loading  fixtures.  A typical  result  obtained  from  the  Brazilian  test  is 
shown  in  Figure  4.13. 

The  variation  of  tensile  strength  with  densification  temperature  for 
gel-silica  samples  densified  in  the  range  of  180  °C  (453  K)  to  1150  °C 
(1423  K)  is  shown  in  Figure  4.14.  There  is  a continuous  increase  in 
tensile  strength  as  the  sintering  temperature  is  more  pronounced  for 
temperatures  above  900  °C  (1173  K).  Figure  4.14  shows  tensile 
strength  values  for  a set  of  fused  quartz  samples,  which  are  around  the 
values  obtained  for  gel-silica  samples  sintered  at  1150  °C  (1423  K). 
The  figure  also  presents  tensile  strength  values  for  dried  samples  of  81 
A average  pore  radius,  which  are  lower  than  the  corresponding  values 
for  samples  of  12  A average  pore  radius. 

4.6  Fracture  Analysis 
4.6.1  Aspects  of  a Fracture  Surface 

Usually  an  analysis  of  the  surface  of  fracture  is  carried  out  with  the 
intention  of  determining  the  cause  of  the  failure.  This  is  important  for 
aiding  in  materials  selection,  for  design  applications,  or  for  improving 
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Figure  4.13  Load  (P)  versus  displacement  for  a dried  silica  gel 
with  12  A average  pore  radius  in  the  Brazilian  test. 
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Figure  4.14  Variation  of  tensile  strength  (Brazilian  test)  as  a function 
of  densification  temperature  for  samples  type  Al,  A2  and  C. 

FQ  represents  fused  quartz  samples. 
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processing  or  operation  conditions.  In  the  case  of  this  research,  the 
purposes  of  the  fracture  analysis  are  the  following: 

1)  Check  the  validity  of  the  tests.  When  the  specimen  breaks  at 
contact  points  with  the  loading  fixtures,  for  instance,  the  results 
are  not  considered.  In  most  cases  a visual  inspection  is  enough 
for  that  purpose. 

2)  Characterize  the  type  of  surface  of  fracture  for  the  silica  gel 
samples  densified  in  a wide  range  of  temperatures. 

3)  Obtain  indications  on  where  the  cracks  initiate  and  how  they 
propagate  during  the  loading. 

Some  features  commonly  observed  in  the  fracture  surface  of  glasses 
are  shown  schematically  in  Figures  4.15  and  4.16  [442,457].  These 
figures  show  schematically  the  presence  in  the  fracture  surface  of  the 
fracture  mirror,  mist  and  hackle  (Figure  4.15),  and  Wallner  lines  and 
striations  (Figure  4.16). 

In  the  following  sections  fracture  surfaces  of  samples  tested  in 
flexural  test  and  Brazilian  test  will  be  shown.  A comprehensive  analysis 
of  fracture  surfaces  for  the  uniaxial  compression  test  will  not  be  carried 
out  because  for  most  samples  a general  crushing  of  the  specimens 
occurred  during  the  test,  leaving  no  convenient  fracture  surface  to  be 
analyzed.  Figure  4.17  shows  the  fracture  surface  of  a dried  silica  gel 
sample  tested  under  uniaxial  compression.  In  this  case,  the  crushing  of 
the  sample  was  not  complete  so  that  a piece  of  the  ruptured  sample 
could  be  analyzed.  Figure  4.17  is  oriented  in  such  a way  that  the  border 
shown  is  perpendicular  to  the  loading  fixtures.  The  figure  suggests  a 
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Figure  4.15  Schematic  of  a fracture  surface  showing  the  region  of 
fracture  origin,  the  fracture  mirror,  the  mist  and  the  hackle 
(Based  on  [442]  and  [457]). 


243 


(a) 


(b) 


Figure  4.16  Schematic  of  fracture  surfaces  showing  (a)  Wallner 
lines  and  (b)  striations  (Based  on  [442]  and  [457]). 
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Figure  4.17  SEM  fracture  surface  of  a dried  silica  gel  sample  with 
12  A average  pore  radius  submitted  to  uniaxial  compression. 
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figure  suggests  a possible  region  of  crack  initiation,  situated  at  the  upper 
right  part  of  the  photomicrograph.  The  results  observed  for  the 
uniaxial  compression  tests,  in  terms  of  mechanics  of  fracture  and 
fracture  propagation  can  be  well  represented  by  Figure  4.6. 

4.6.2  Flexural  Testing 

Analysis  of  fractured  surfaces  of  mechanically  tested  (flexural  and 
compression)  samples  was  carried  out  in  a scanning  electron  microscope 
(SEM)  JEOL-JSM-35CF.  All  samples  were  coated  with  a Au-Pd  layer. 
A typical  fracture  surface  of  a dried  gel-silica  sample  tested  in  four 
points  bending  is  shown  Figure  4.18.  The  figure  shows  the  presence  of 
a sharp  edge  formed  during  the  fracture  process,  suggesting 
delamination  of  the  sample  caused  by  the  bending  test,  besides  a possible 
effect  of  the  state  of  stress  during  the  test.  The  same  type  of  effect  is 
present  for  every  temperature  tested,  as  exemplified  by  Figure  4.19 
(sample  sintered  at  900  °C,  1173  K). 

The  delamination  phenomena  suggests  a structural  gradient  inside  the 
sample.  Indeed,  microhardness  measurements  along  the  transverse 
direction  in  the  sample  show  a profile  in  which  the  hardness  is  higher  in 
the  center  when  compared  to  the  hardness  of  the  external  surface  of  the 
sample.  The  variation  of  VHN  with  the  transverse  position  in  a dried 
sample  is  shown  in  Figure  4.20.  At  higher  temperature  of  sintering  the 
hardness  gradient  increases,  as  seen  in  Figure  4.21  for  a sample 
densified  at  900  °C  (1173  K). 
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Figure  4.18  SEM  fracture  surfaces  of  a dried  silica  gel  sample  with 
12  A average  pore  radius  submitted  to  four  point  bending.  The 
photos  represent  the  opposing  sides  of  the  fracture  surface. 
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Figure  4.19  SEM  fracture  surfaces  of  a silica  gel  sample  with 
12  A average  pore  radius  densified  at  900  °C  and  submitted 
to  four  point  bending.  The  photos  represent  the 
opposing  sides  of  the  fracture  surface. 


NHA 


248 


POSITION  (mm) 


Figure  4.20  Variation  of  the  Vickers  hardness  number  (VHN)  with 
transverse  position  in  a dried  silica  gel  sample  of  12  A average 
pore  radius.  Positions  0 and  4.5  represent  the  external 
limiting  surfaces  of  the  specimen. 
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Figure  4.21  Variation  of  the  Vickers  hardness  number  (VHN)  with 
transverse  position  in  a silica  gel  sample  of  12  A average  pore 
radius  densified  at  900  °C.  Positions  0 and  4.5  represent  the 
external  limiting  surfaces  of  the  specimen. 
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A possible  cause  for  the  gradient  in  hardness  is  associated  with  the 
propagation  of  a diffusion  front  formed  by  the  gas-liquid  interface 
during  the  drying  of  the  sample.  As  the  drying  front  propagates,  from 
the  external  surface  towards  the  center  of  the  sample,  it  creates  a 
structural  gradient  which  is  amplified  during  the  densification  process. 

The  fracture  behavior  of  silica  gel  monolith  samples  tested  in  four 
point  bending,  as  observed  in  these  experiments,  can  be  summarized  by 
Figure  4.22.  Figure  4.22(a)  represents  the  regions  under  compression 
and  tension  during  the  application  of  load.  The  crack  initiates  at  the 
tensile  region.  Figure  4.22(b)  represents  a typical  shape  of  the  fracture 
surface  observed.  This  shape  is  probably  influenced  by  a structural 
gradient  inside  the  sample.  Figures  4.22(c)  and  4.22(d)  represent  the 
regions  where  cracks  were  observed  to  initiate. 

4.6.3  Brazilian  Testing 

Scanning  electron  microscopy  analysis  of  the  fracture  surface  of  gel- 
silica  samples  densified  in  the  range  of  180  °C  (453  K)  to  1150  °C  (1423 
K)  tested  in  the  Brazilian  test  show  characteristics  of  brittle  fracture, 
i.e,  the  presence  of  no  macroscopic  deformation.  All  tested  samples 
show  characteristic  striations  (or  river  patterns),  as  shown  in  Figures 
4.23  to  4.29. 

Fracture  surface  of  dried  gel-silica  sample  is  shown  in  Figure  4.23. 
Figure  4.23  was  taken  in  one  of  the  borders  of  the  sample.  It  indicates 
the  presence  of  river  patterns  in  a direction  perpendicular  to  the  loading 
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Figure  4.22  Schematic  of  the  fracture  in  four  point  bending,  (a) 
Compression  and  tension  regions,  (b)  Fracture  surfaces,  (c)  Crack 
initiation  at  the  middle  of  the  surface  under  tension. 

(d)  Crack  initiation  at  one  comer  of  the  specimen. 
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Figure  4.23  SEM  fracture  surface  of  a dried  silica  gel  sample  with 
12  A average  pore  radius  submitted  to  Brazilian  test. 
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Figure  4.24  SEM  fracture  surface  of  a silica  gel  sample  with  12  A 
average  pore  radius  densified  at  600  °C  and  submitted  to  Brazilian  test. 
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Figure  4.25  SEM  fracture  surface  of  a silica  gel  sample  with  12  A 
average  pore  radius  densified  at  600  °C  and  submitted  to  Brazilian  test. 
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Figure  4.26  SEM  fracture  surface  of  a silica  gel  sample  with  12  A 
average  pore  radius  densified  at  700  °C  and  submitted  to  Brazilian  test. 
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Figure  4.27  SEM  fracture  surface  of  a silica  gel  sample  with  12  A 
average  pore  radius  densified  at  900  °C  and  submitted  to  Brazilian  test. 
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Figure  4.28  SEM  fracture  surface  of  a silica  gel  sample  with  12  A 
average  pore  radius  densified  at  1150  °C  and  submitted  to  Brazilian  test. 
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Figure  4.29  SEM  fracture  surface  of  a silica  gel  sample  with  12  A 
average  pore  radius  densified  at  1150  °C  and  submitted  to  Brazilian  test. 
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fixtures.  Figure  4.24  shows  the  radial  propagation  of  river  patterns  out 
of  a circular  type  of  mirror  present  in  the  fracture  surface  of  a silica 
gel  sample  densified  at  600  °C  (873  K).  Figure  4.25  shows  an  aspect  of 
the  striations  present  in  the  fracture  surface  of  a silica  gel  sample 
densified  at  600  °C  (873  K). 

The  fracture  surface  of  a silica  gel  sample  densified  at  700  °C  (973 
K)  is  shown  in  Figure  4.26.  This  figure  shows  the  striations  originating 
in  the  center  (region  under  tension)  of  the  specimen,  and  propagating 
towards  the  border,  where  the  sample  touches  the  loading  fixture.  The 
fracture  surface  of  a silica  gel  sample  densified  at  900  °C  (1173  K)  is 
shown  in  Figure  4.27.  As  in  Figures  4.23  and  4.26,  Figure  4.27 
indicates  that  the  fracture  in  the  samples  tested  in  the  Brazilian  test 
originates  in  the  center  of  the  sample,  where  the  tension  stress  is 
maximum.  Figure  4.28  shows  the  fracture  surface  of  a sample  densified 
at  1150  °C  (1423  K),  with  the  presence  of  striations.  Figure  4.29  shows 
an  aspect  of  the  fracture  surface  of  a silica  gel  sample  densified  at  1150 
°C  (1423  K),  with  the  presence  of  "smooth"  and  "sharp  edge"  striations. 
Figure  4.30  shows  the  fracture  surface  of  a fused  quartz  sample  which 
also  shows  river  patterns. 

The  striations  are  present  in  the  fracture  surfaces  under  different 
aspects.  Figure  4.31(a)  shows  a river  pattern  present  in  the  fracture 
surface  of  a dried  silica  gel.  Figure  4.31(b)  shows  a magnified  detail  of 
Figure  4.31(a). 

A schematic  of  the  observed  fracture  surfaces  of  the  silica  gel 
samples  tested  in  the  Brazilian  test  is  shown  in  Figure  4.32. 
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Figure  4.30  SEM  fracture  surface  of  a fused  quartz  sample  submitted 

to  Brazilian  test. 
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Figure  4.31  SEM  fracture  surface  of  a dried  silica  gel  sample  with  12 
A average  pore  radius  submitted  to  Brazilian  test,  (a)  General  view;  (b) 
detail  of  the  secondary  striations  shown  in  (a). 
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Figure  4.32  Schematic  of  the  observed  fracture  surfaces  of  silica  gel 
samples  in  the  Brazilian  test,  (a)  Loading,  (b)  Fracture  surface  (the  size 
of  the  striations  is  shown  exaggerated  for  clarity). 
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In  general,  it  can  be  seen  that  from  the  dried  stage  (180  °C,  453  K) 
throughout  densification,  all  samples  exhibit  brittle  fracture 
characteristics,  namely  the  absence  of  macroscopic  deformation  and  the 
presence  of  river  patterns  in  cleavage-like  (glasses  do  not  have  actual 
cleavage  planes)  surfaces  of  fracture.  Even  though  the  densification 
process  increases  the  strength  of  silica-gel  monoliths,  it  apparently  does 
not  change  the  mechanics  of  fracture  of  the  material,  at  least  to  a extent 
in  which  such  a change  could  be  directly  observed  on  the  fracture 
surface. 


4.7  Dynamic  Mechanical  Analysis 

Dynamic  mechanical  analysis  (DMA)  is  a technique  of  mechanical 
spectroscopy.  It  measures  the  energy  stored  and  dissipated  during 
deformation. 

A damped  oscillation  in  the  material  will  be  created  when  it  is 
deformed  and  then  released  of  the  load.  The  energy  stored  during 
deformation  will  be  released  in  such  a way  that  is  a characteristic  of 
each  material.  The  frequency  of  the  damped  oscillation  is  a function  of 
the  elastic  modulus  of  the  material  [458]. 

In  an  ideal  elastic  material  all  the  energy  stored  is  elastic.  In  a real 
material,  a viscoelastic  behavior  is  observed  in  which  part  of  the  stored 
energy  will  be  dissipated  in  other  forms,  such  as  heat  [458].  The 
greater  the  tendency  for  energy  dissipation,  the  larger  the  damping  of 
the  induced  oscillation. 
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The  resonant  frequency  (fQ)  of  the  mechanical  system  formed  by  the 
DMA  and  the  sample  is  related  to  the  Young's  modulus  (E)  of  the 
sample  as  follows  [458]: 


where  J0  is  the  moment  of  inertia  of  the  arm  that  clamps  the  sample; 

k is  the  spring  constant  of  the  flexure  pivots  of  the  clamp  arms; 

Lv  is  the  clamping  distance; 

H is  the  sample  width; 

D is  the  sample  length; 

B is  the  sample  thickness. 

More  details  on  the  DMA  and  on  the  equations  used  in  dynamic 
mechanical  analysis  are  found  in  references  458  to  463. 

Gel-silica  samples  with  densification  temperatures  varying  from  180 
°C  (453  K)  to  1000  °C  (1273  K)  were  tested  in  a dynamic  mechanical 
analyzer  machine,  type  DMA-982  (Dupont  Instruments).  The  tests  were 
performed  at  temperature  varying  from  25  °C  (298  K)  to  400  °C  (673 
K). 

Figure  4.33  shows  the  variation  of  tensile  modulus  at  room 
temperature  for  a set  of  gel-silica  samples  with  12  A average  pore 
radius  as  a function  of  densification  temperature.  The  tensile  modulus 
increases  continuously  as  the  densification  temperature  increases,  and 
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Figure  4.33  Variation  of  tensile  modulus  (DMA)  at  25  °C  as  a 
function  of  densification  temperature  for  silica  gel  samples 
of  12  A average  pore  radius. 
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the  increase  is  more  pronounced  for  sintering  temperatures  above  800 
°C  (1073  K). 

The  variation  of  tensile  modulus  measured  in  the  range  of  25  °C 
(298  K)  to  400  °C  (673  K)  for  samples  with  12  A average  pore  radius 
samples  densified  in  the  range  180  °C  (453  K)  to  1000  °C  (1273  K)  is 
shown  in  Figure  4.34.  As  the  testing  temperature  increases  the  tensile 
modulus  decreases  for  all  samples.  Figure  4.35  presents  similar  results 
for  the  shear  modulus. 

It  is  worthwhile  to  note  that  the  values  of  tensile  modulus  obtained 
from  the  DMA  tests  (Figure  4.33)  are  very  close  to  the  compression 
modulii  obtained  with  the  uniaxial  compression  tests  (Figures  4.9  and 
4.11).  The  values  of  tensile  modulus  obtained  from  the  Brazilian  tests 
were  very  similar  to  the  results  obtained  from  the  DMA  and  uniaxial 
compression  tests. 
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Figure  4.34  Variation  of  tensile  modulus  (DMA)  as  a function  of 
densification  temperature  for  silica  gel  samples  of  12  A average  pore 
radius.  The  temperatures  in  the  legend  indicate  the  test  temperature. 


SHEAR  MODULUS  (GPa) 


268 


TEMPERATURE  (°C) 


Figure  4.35  Variation  of  shear  modulus  (DMA)  as  a function  of 
densification  temperature  for  silica  gel  samples  of  12  A average  pore 
radius.  The  temperatures  in  the  legend  indicate  the  test  temperature. 


CHAPTER  5 
THERMAL  ANALYSIS 


5.1  Thermal  Expansion  of  Materials 


There  are  several  variables  that  can  influence  the  thermal  expansion 
behavior  of  a material.  The  level  of  such  a influence  varies  from 
material  to  material.  Among  those  variables,  the  following  usually  have 
some  effect  on  the  thermal  expansion: 

1)  State  of  the  phases  (crystalline  or  amorphous)  [464-467]. 

2)  Residual  stress  [468-470]. 

3)  Chemistry  of  the  surface  [471-473]. 

4)  Fictive  temperature  of  glassy  phase  [474-477]. 

5)  Occurrence  of  cracks  [478]. 

6)  Density  [476,479]. 

In  some  cases  the  level  of  porosity  does  not  seem  to  affect 
significantly  the  thermal  expansion  [480-482]. 


The  coefficient  of  volumetric  thermal  expansion  (av)  is  defined  as 
[483] 


(5.1) 


where  V is  the  volume; 
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T is  the  absolute  temperature; 

P is  the  (constant)  pressure. 

Replacing  in  Eq.(5.1)  the  volume  (V)  by  the  linear  dimension  (L)  the 
coefficient  of  linear  thermal  expansion  (ccl)  can  be  written  as  follows: 


For  isotropic  materials  (like  glasses,  for  instance)  av  = 3ccl. 

The  coefficient  of  linear  thermal  expansion  (av)  is  referred  to 
simply  as  coefficient  of  thermal  expansion,  or  CTE,  in  most  of  this 
work.  The  average  coefficient  of  thermal  expansion  is  called  CTEm  in 
this  work.  The  domain  of  averaging  for  CTEm  is  from  30  °C  (303  K) 
to  the  specified  testing  temperature. 

Usually  material  expands  when  the  temperature  is  increased.  This  is 
a direct  consequence  of  the  asymmetry  of  the  potential  energy  curve  for 
two  adjacent  atoms  [484].  When  the  temperature  increases,  the 
amplitude  of  the  vibrating  atoms  also  increase.  Since  the  amplitude  of 
the  vibrating  atom  is  symmetric  about  a median  position,  and  at  the 
same  time  the  potential  energy  curve  is  asymmetric,  the  vibrating  atom 
moves  further  apart  from  its  neighbor,  causing  the  thermal  expansion. 

There  are  many  different  methods  available  to  measure  the  thermal 
expansion  of  solids,  as  follows: 


(5.2) 


1)  Push-rod  methods  [485-487]. 

2)  Optical  methods  [488-496]. 

3)  Interferometer  methods  [497-503]. 
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4)  X-ray  methods  [504-506]. 

5)  Parallel  plate  capacitor  methods  [507-511]. 

The  coefficient  of  thermal  expansion  (CTE)  of  materials  usually 
varies  depending  on  the  temperature  of  the  measurement.  Hummel 
[512]  reports  an  average  CTE  for  fused  silica  of  0.5x1 0-6  K_1  (in  the 
range  20-1000  °C,  293-1273  K);  Zwikker  [513]  reports  a CTE  of 
0.3xl0-6  K-1  (at  0 °C,  273  K).  Many  investigators  [514-519]  use  a 
value  around  0.4x1 0-6  K_1  for  the  CTE  of  fused  silica  at  25  °C  (298  K). 
These  values  are  very  low  when  compared  to  other  glasses  (e.g.,  for 
borosilicate  glass  the  CTE  at  0 °C  (273  K)  is  3xl0-6  K_1  [513]). 
According  to  Biscoe  et  al.  [520], 

the  small  coefficient  of  thermal  expansion  of  vitreous  silica  is  due  to 
the  fact  that  in  the  continuous  random  network  a change  in  volume 
can  come  only  from  a change  in  the  silicon-oxygen  distance.  The 
silicon-oxygen  bond  is  such  a strong  one  that  any  such  change  will  be 
small. 

Small  differences  in  the  CTE  of  vitreous  silica  occur  depending  on 
the  commercial  grade  of  the  glass  [521].  Different  commercial  grades 
will  present  different  levels  of  impurities,  bubbles,  and  hydroxyl 
groups,  all  of  which  affect  the  thermal  expansion  of  the  glass  [521]. 

The  CTE  of  fully  dense  silica  gel  is  reported  to  be  smaller  than  the 
CTE  for  fused  silica  glass  [29].  A possible  reason  for  that  apparent 
behavior  was  assumed  to  be  a lower  concentration  of  cation  impurities 
and  a larger  intermolecular  space  between  the  structural  units  in  the 
fully  dense  silica  gel  [29]. 
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According  to  Hench  and  West  [28]  the  average  bond  length  between 
neighboring  silicon  atoms  increases  when  water  is  absorbed  in  the 
structure.  Because  of  this,  the  dimensional  change  of  a silica  gel  can  be 
divided  into  two  components:  (a)  an  intrinsic  thermal  expansion,  due  to 
the  thermal  behavior  of  the  structural  network  of  the  silica  gel;  and  (b) 
an  extrinsic  component,  due  to  the  effects  of  the  adsorbed  water  [28]. 

Since  factors  other  than  the  intrinsic  thermal  expansion  can  influence 
the  behavior  of  a material  when  it  is  heated,  the  more  general 
designation  "dimensional  change"  (DC)  will  be  used  in  this  work  instead 
of  thermal  expansion.  The  expression  "coefficient  of  dimensional 
change"  (CDC)  will  be  used  in  this  work,  meaning  the  instantaneous 
slope  of  the  curve  dimensional  change  versus  testing  temperature.  In  a 
similar  manner,  the  expression  "average  coefficient  of  dimensional 
change"  (CDCm)  will  represent  the  slope  of  the  curve  dimensional 
change  versus  temperature  when  averaged  in  the  range  of  30  °C  (303  K) 
to  the  specified  testing  temperature.  While  CDC  is  a more  general 
expression  for  CTE,  CDCm  is  a more  general  expression  for  CTEm. 

The  purpose  of  the  dimensional  change  analysis  carried  out  in  this 
research  is  to 

1)  Characterize  the  dimensional  change  behavior  of  silica  gel 
samples  densified  in  a wide  range  of  temperatures  (180-1150  °C, 
453-1423  K). 

2)  Use  the  dimensional  change  behavior  as  an  indicator  of  the  degree 


of  densification. 
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3)  Obtain  indirect  information  on  the  structural  transformations 
during  sintering  of  silica  gel  monoliths. 

5.2  Dilatometrv 

5.2.1  Dimensional  Change  Behavior  of  Silica  Gel  Monoliths 

The  dimensional  change  experiments  reported  in  this  work  were 
carried  out  using  an  automatic  recording  differential  push-rod 
dilatometer  Dilatronic  II-R  from  Theta  Industries. 

A typical  curve  of  dimensional  change  as  a function  of  testing 
temperature  for  a gel-silica  sample  densified  at  600  °C  (873  K)  is  shown 
in  Figure  5.1.  The  (linear)  dimensional  changes  involved  in  most  of  the 
experiments  reported  in  this  research  superimpose  the  effects  of 
shrinkage  (due  to  densification)  and  intrinsic  thermal  expansion.  As 
seen  in  Figure  5.1,  the  initial  condition  is  at  room  temperature,  in  which 
the  sample  is  at  zero  dimensional  change.  When  the  sample  is  heated  at 
a rate  of  3 °C/min,  it  begin  to  expand  up  to  a maximum  expansion 
which  occurs  at  a certain  temperature,  Tmax.  Beyond  Tmax  the  gel- 
silica  monolith  start  shrinking  and  the  expansion  crosses  the  line  of  zero 
expansion  at  a certain  temperature,  Tzero*  As  the  testing  temperature 
increases  to  700  °C  (973  K),  the  sample  continues  to  shrink.  When  the 
temperature  reaches  700  °C  (973  K)  the  furnace  (and  the  sample)  begin 
the  cooling  step,  at  a cooling  rate  of  3 °C/min.  At  high  temperatures. 
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Figure  5.1  Variation  of  the  dimensional  change  with  testing 
temperature  for  samples  type  A4  previously  densified  at  600  °C. 
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the  sample  still  shrinks  significantly  during  cooling,  but  below  a certain 
temperature  the  size  of  the  sample  remains  fairly  unchanged  upon 
cooling.  When  the  cooling  temperature  reaches  about  150  °C  (423  K),  a 
small  expansion  occurs,  because  of  water  adsorption. 

The  variation  of  dimensional  change  with  testing  temperature  for  12 
A average  pore  radius  samples  densified  in  the  range  180  °C  (453  K)  to 
1150  °C  (1423  K)  is  shown  in  Figure  5.2.  It  is  seen  in  this  figure  that 
the  dried  sample  shrinks  the  most  upon  heating  and  the  amount  of 
shrinkage  decreases  as  the  sintering  temperature  increases.  The  area 
between  the  curves  for  heating  and  cooling  (area  of  the  thermal 
hysteresis  loop)  can  be  used  to  compare  thermal  stability  among  sets  of 
samples.  Another  way  comparing  thermal  stability  is  describing  the 
dimensional  change  at  a given  test  temperature,  e.g.  at  690  °C  (963  K). 

Similarly,  the  coefficient  of  dimensional  change  (CDC)  is  useful  not 
only  to  characterize  the  dimensional  change  behavior  of  a material  but 
as  a tool  to  compare  the  evolution  of  thermal  properties  as  the  sintering 
proceeds.  The  variation  of  the  average  coefficient  of  dimensional 
change  (CDCm)  with  testing  temperature  for  12  A average  pore  radius 
samples  previously  densified  in  the  range  of  180  °C  (453  K)  to  1150  °C 
(1423  K)  is  shown  in  Figure  5.3.  The  coefficient  of  dimensional  change 
is  continuously  less  negative  as  the  densification  temperature  increases. 

Analyzing  the  curve  for  the  dried  sample,  one  sees  that  at  the 
beginning  of  the  heating  program,  the  coefficient  of  dimensional  change 
is  positive,  corresponding  in  the  Figure  5.2  to  the  initial  stage  where 
there  is  a positive  dimensional  change  in  the  sample.  The  coefficient  of 
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Figure  5.2  Variation  of  the  dimensional  change  with  testing 
temperature  for  samples  type  A4  previously  densified  in 
the  range  of  180  °C  to  1150  °C. 
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Figure  5.3  Variation  of  the  average  coefficient  of  dimensional  change 
(CDCm)  as  a function  of  testing  temperature  for  samples  type  A4 
previously  densified  in  the  range  of  180  °C  to  1 150  °C. 
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dimensional  change  decreases  continuously,  and  after  a certain 
temperature  it  becomes  negative,  indicating  shrinkage  of  the  sample  in 
reference  to  its  initial  size  before  the  test  begin.  As  the  temperature  of 
the  test  goes  up,  the  coefficient  of  dimensional  change  becomes 
increasingly  negative,  indicating  increasing  shrinkage  of  the  sample. 

5.2.2  Large  Pore  Size  Structures 

Similarly  to  the  small  radius  samples,  the  81  A pore  radius  samples 
present  greater  thermal  stability  as  the  densification  temperature 
increases  (Figure  5.4).  Using  the  dimensional  change  at  690  °C  (963  K) 
as  an  indication  of  the  thermal  stability  of  the  samples,  the  81  A pore 
radius  gel-silicas  are  thermally  more  stable  than  the  12  A pore  radius 
samples,  as  shown  in  Figure  5.5. 

Figure  5.6  shows  the  variation  of  the  average  coefficient  of 
dimensional  change  (CDCm)  with  testing  temperature  for  81  A pore 
radius  samples  densified  in  the  range  180  °C  (453  K)  to  1150  °C  (1423 
K).  Considering  the  curve  of  variation  of  CDCm  with  temperature  for 
the  dried  sample  (180  °C),  one  sees  that  at  low  temperatures  the  CDCm 
is  positive,  indicating  that  the  specimen  is  expanded  compared  to  its 
initial  length.  As  the  testing  temperature  increases,  the  coefficient  of 
dimensional  change  becomes  increasingly  negative,  indicating  shrinkage. 
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Figure  5.4  Variation  of  the  dimensional  change  with  testing 
temperature  for  silica  gel  samples  with  81  A average  pore 
radius  previously  densified  in  the  range  of  180  °C  to  1150  °C. 
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Figure  5.5  Variation  of  the  dimensional  change  at  690  °C  as  a function 
of  densification  temperature  for  samples  type  A4  and  C. 
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Figure  5.6  Variation  of  the  average  coefficient  of  dimensional  change 
(CDCm)  as  a function  of  testing  temperature  for  silica  gel  samples  with 
81  A average  pore  radius  densified  in  the  range  of  180  °C  to  1150  °C. 
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5.2.3  Isothermal  Experiments 

In  order  to  evaluate  the  influence  of  time  on  the  sintering  behavior 
of  gel-silica  monoliths,  samples  were  held  in  the  dilatometer  at  700  °C 
(973  K)  for  10  hours.  The  heating  rate  (3  °C/min)  was  the  same  as  the 
previous  experiments. 

The  variation  of  dimensional  change  with  time  at  700  °C  (973  K)  for 
12  A pore  radius  samples  densified  at  900  °C  (1173  K)  and  1150  °C 
(1423  K)  is  shown  in  Figure  5.7.  For  the  sample  densified  at  1150  °C 
(1423  K)  as  the  testing  temperature  increases  (time  0 to  228  min)  there 
is  a positive  dimensional  change.  After  the  temperature  reaches  700  °C 
(973  K)  the  sample  stops  expanding  and  its  dimensions  remain 
approximately  unchanged  with  time.  After  10  hours  at  700  °C  (973  K) 
(at  the  time  828  min),  the  cooling  step  begins  and  the  sample  shrinks 
back  to  its  original  dimensions.  Considering  now  the  sample  densified 
at  900  °C  (1173  K),  it  can  be  seen  from  Figure  5.7  that  it  shrinks  upon 
heating  from  room  temperature  to  700  °C  (973  K),  corresponding  to  a 
variation  of  time  from  zero  to  228  min.  When  the  heating  temperature 
reaches  700  °C  (973  K)  and  is  held  at  that  level,  the  rate  of  shrinkage 
changes  dramatically.  From  time  228  min  to  the  end  of  the  isothermal 
treatment  at  time  828  min,  the  sample  densified  at  900  °C  (1173  K) 
shrinks  slowly.  The  silica  gel  sample  does  expand  slightly  upon  cooling 
from  700  °C  (973  K)  to  room  temperature  due  to  water  adsorption. 

Figure  5.8  shows  the  variation  of  dimensional  change  with  time,  at  a 
temperature  of  700  °C  (973  K),  for  a sample  densified  at  1150  °C 
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Figure  5.7  Variation  of  the  dimensional  change  with  time  for  samples 
type  A4  densified  at  900  °C  and  1150  °C.  The  samples  were  heated  to 
700  °C  in  228  min  and  cooled  to  room  temperature  after  10  h at  700  °C. 
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Figure  5.8  Variation  of  the  dimensional  change  with  time  for  a 
sample  type  A4  densified  at  1 150  °C  and  a NBS  silica  standard. 
The  samples  were  heated  to  700  °C  in  228  min  and  cooled 
to  room  temperature  after  10  h at  700  °C. 


285 


(1423  K)  compared  with  a NBS  silica  standard.  In  the  case  of  these 
samples,  the  dimensional  change  corresponds  to  intrinsic  thermal 
expansion.  Those  samples  present  very  similar  behavior  during  the 
isothermal  treatment  and  the  total  expansion  of  both  samples  are  very 
close.  Again,  the  expansion  between  time  zero  and  228  min  is  due  to 
the  heating  step  from  room  temperature  to  700  °C  (973  K).  The 
shrinkage  after  time  828  min  is  due  to  the  cooling  step. 

As  shown  previously  in  Figure  5.5,  the  gel-silica  samples  of  81  A 
pore  radius  are  thermally  more  stable  than  the  12  A pore  radius 
samples.  That  same  observation  can  be  drawn  by  using  isothermal 
experiments.  As  seen  in  Figure  5.9,  the  12  A pore  radius  dried  sample 
shrinks  much  more  than  the  81  A pore  radius  dried  sample  under  the 
same  temperature  profile.  Again  the  influence  of  time  on  the  shrinkage 
of  gel-silica  can  be  seen  in  Figure  5.9.  As  the  time  at  700  °C  (973  K) 
increases  (from  time  228  min  to  828  min)  both  (81  A and  12  A) 
samples  shrink  but  at  a much  lower  rate  than  during  the  heating  period. 

5.2.4  Aspects  of  the  Dimensional  Change  Curve 

As  described  before  (section  5.2.1)  for  the  dried  gel-silica,  as  the 
testing  temperature  increases,  the  specimen  first  expands  and  after  a 
certain  temperature  (Tmax)  it  begins  to  shrink.  The  sample  reaches  its 
original  length  at  a temperature  Tzero,  and  upon  further  heating  the 
sample  continues  to  shrink  until  the  maximum  test  temperature  of  700 
°C  (973  K).  Figure  5.10  shows  the  variation  of  the  temperature  for 
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Figure  5.9  Variation  of  the  dimensional  change  with  time  for  dried 
silica  gel  samples  types  A4  and  C.  The  samples  were  heated  to  700  °C  in 
228  min  and  cooled  to  room  temperature  after  10  h at  700  °C. 
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Figure  5.10  Variation  of  the  temperature  for  zero  expansion 
during  heating  (Tzero)  and  temperature  for  maximum  expansion 
during  heating  (Tmax)  as  a function  of  densification 
temperature  for  samples  type  A4. 
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maximum  expansion  (Tmax)  and  the  temperature  for  zero  expansion 
(Tzero)  as  a function  of  densification  temperature  for  12  A average  pore 
radius  samples.  It  can  be  seen  from  Figure  5.10  that  Tmax  and  Tzero 
decreases  as  the  densification  temperature  increases.  For  samples 
densified  at  temperatures  of  900  °C  (1173  K)  and  1000  °C  (1273  K)  the 
initial  expansion  cannot  be  seen  using  the  initial  testing  temperature  of 
25  °C  (298  K).  The  sample  densified  at  1150  °C  (1423  K)  does  not 
present  Tmax  and  Tzero  because  it  expands  continuously  from  room 
temperature  until  700  °C  (973  K). 

5.3  Thermogravimetric  Analysis 

Thermogravimetric  analysis  (TGA)  was  performed  in  a 
Thermogravimetric  Analyzer  Model  951  (Dupont  Instruments).  Silica 
gel  monolithic  samples  densified  in  the  range  of  180  °C  (453  K)  to  1150 
°C  (1423  K)  were  crushed  to  fine  powders  and  submitted  to  a TGA 
analysis  in  the  testing  temperature  range  of  25  °C  (298  K)  to  1150  °C 
(1423  K). 

Before  testing  the  samples,  the  powders  were  put  in  an  oven  at  140 
°C  for  a period  of  one  day,  in  order  to  obtain  a homogeneous  set  of 
samples. 

A typical  TGA  plot  for  a dried  gel-silica  is  shown  in  Figure  5.11. 
That  figure  shows  the  increasing  weight  loss  that  occurs  when  the  dried 
gel-silica  is  heated.  There  is  a sharp  decrease  in  the  weight  of  the 
sample  when  it  is  heated  from  room  temperature  to  about  150  °C 
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Figure  5.1 1 Variation  of  the  TGA  weight  loss  as  a function  of  testing 
temperature  for  a dried  silica  gel  sample  type  A4. 
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(423  K).  Beyond  150  °C,  the  weight  continues  to  decrease,  but  at  a 
much  slower  rate. 

The  pronounced  weight  loss  that  occurs  before  150  °C  (423  K)  is 
associated  with  elimination  of  water  from  the  sample.  The  weight  loss 
that  occurs  beyond  150  °C  (423  K)  is  usually  associated  with  loss  of 
water  and  other  volatile  materials  in  the  sample,  such  as  organic 
residues. 

The  variation  of  TGA  weight  loss  as  a function  of  densification 
temperature  for  a set  of  silica  gel  samples  of  12  A pore  radius  (A2)  is 
shown  in  Figure  5.12.  The  figure  shows  both  the  weight  loss  that 
occurs  from  room  temperature  to  1100  °C  (1373  K),  and  the  weight 
loss  from  150  °C  (423  K)  to  1100  °C  (1373  K).  For  simplification,  the 
difference  between  the  two  types  of  weight  loss  shown  in  Figure  5.12  is 
called  in  this  work  the  "water  content"  of  the  silica  gels.  The  expression 
"water  content"  is  used  here  in  regard  to  the  large  contribution  in  the 
weight  loss  due  to  the  physical  water  content  when  the  silica  gels  are 
heated  from  room  temperature  to  about  150  °C  (423  K).  When  the 
silica  gels  are  heated  above  150  °C  (423  K),  water  still  contributes  to  the 
weight  loss,  but  mainly  in  the  form  of  elimination  of  OH  groups. 
Figure  5.13  shows  that  as  the  densification  temperature  increases  the 
water  content  in  the  silica  gels  decreases.  This  fact  can  be  directly 
associated  with  the  dimensional  change  behavior  shown  in  Figure  5.10. 
The  results  shown  in  Figures  5.10  and  5.13  suggest  that  the  initial 
expansion  of  the  samples  densified  at  lower  temperatures  is  associated 
with  the  presence  of  water  in  the  samples.  As  the  temperature  of  the 
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Figure  5.12  Variation  of  the  TGA  weight  loss  measured  from 
25  °C  and  from  150  °C  as  a function  of  testing  temperature  for  a 
dried  silica  gel  sample  type  A4. 
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Figure  5.13  Variation  of  the  water  content  as  a function  of 
densification  temperature  for  samples  type  A4. 
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specimen  increases  the  water  present  in  the  sample  expands, 
contributing  to  the  observed  expansion.  Figure  5.14  shows  the  variation 
of  the  temperature  of  maximum  expansion  (Tmax)  as  a function  of  the 
water  content  in  silica  gel  monoliths  of  12  A average  pore  radius. 
Figure  5.15  shows  the  variation  of  the  temperature  of  zero  expansion 
(Tzero)  as  a function  of  the  water  content  for  the  same  samples  referred 
to  in  Figure  5.14.  The  variation  of  the  dimensional  change  measured  at 
Tmax  as  a function  of  the  water  content  in  silica  gel  monoliths  with  12  A 
average  pore  radius  is  shown  in  Figure  5.16. 

Figure  5.17  shows  the  variation  of  TGA  weight  loss  as  a function  of 
densification  temperature  for  a different  set  (A  1)  of  12  A pore  radius 
samples.  The  same  tendency  of  decreasing  weight  loss  as  the 
densification  temperature  increases  is  observed  for  that  set  of  samples. 

According  to  the  results  obtained  in  this  chapter,  the  dimensional 
changes  that  occur  during  the  heating  of  partially  dense  silica  gel 
monoliths  have  at  least  four  components: 

1)  Intrinsic  thermal  expansion. 

2)  Effects  of  the  presence  of  water. 

3)  Shrinkage  due  to  densification. 

4)  Effect  of  structural  defects  (bubbles,  cracks,  etc.). 

In  a fully  dense  sample,  with  no  adsorbed  water  and  free  of  defects, 
the  dimensional  changes  that  occur  during  heating  coincide  with  the 
intrinsic  thermal  expansion  effects.  This  is  shown  by  Figure  5.8, 
assuming  that  the  NBS  standard  is  defect  free  and  has  no  adsorbed 


water. 
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Figure  5.14  Variation  of  Tmax  as  a function  of  the  water  content 
in  silica  gel  samples  type  A4. 
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Figure  5.15  Variation  of  Tzero  as  a function  of  the  water  content 
in  silica  gel  samples  type  A4. 
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Figure  5.16  Variation  of  the  dimensional  change  measured  at  Tmax  as  a 
function  of  the  water  content  in  samples  type  A4. 


297 


0 200  400  600  800  1000  1200 

TEMPERATURE  (°C) 


Figure  5.17  Variation  of  the  TGA  weight  loss  measured  from 
25  °C  and  from  150  °C  as  a function  of  testing  temperature 

for  samples  type  Al. 
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For  a fully  dense  silica  gel  the  values  of  coefficient  of  dimensional 
change  (CDC)  coincide  with  the  values  of  coefficient  of  thermal 
expansion  (CTE).  In  order  to  calculate  CTE’s  for  partially  densified 
silica  gel  samples,  the  effects  of  the  presence  of  water,  shrinkage  and 
structural  defects  have  to  be  taken  into  account.  Therefore,  the 
observed  variation  in  length  (Al)  in  silica  gel  monoliths  can  be 
expressed  as  follows: 

Al  = Alt  + Als  + Alw  + Aid,  (5.3) 

where  Alt  is  the  thermal  expansion  component; 

Als  is  the  shrinkage  component; 

Alw  is  the  component  due  to  the  presence  of  water; 

Aid  is  the  component  due  to  the  presence  of  defects. 


CHAPTER  6 
STRUCTOLOGY 


6.1  Concept  of  Structologv 

In  order  for  a structure  to  be  completely  described,  both  metric  and 
topological  properties  have  to  be  taken  into  account.  The  science  of 
correlating  structural  parameters  with  the  physical  properties  and, 
ultimately,  with  the  behavior  of  materials  was  named  by  Rhines  [365]  as 
"microstructology."  However,  because  of  the  many  different  scales  in 
which  structures  exist  (e.g.,  molecules,  angstrom-structures;  silica  gels, 
nanometer-structures;  polycrystalline  ceramics  and  metals,  micrometer- 
structures;  galaxies,  megaparsec-structures  [522]),  the  designation 
"structology"  seems  to  be  more  appropriate.  In  this  sense  structure 
means  the  set  of  topological  and  metric  properties  necessary  to 
characterize  the  geometry  of  a given  n-dimensional  manifold. 

Structology  deals  with  the  definition  and  experimental  determination 
of  mathematical  relationships  between  the  structural  parameters  and  the 
physical  properties  of  a material,  as  well  as  the  description  of  variations 
of  structural  properties  as  a function  of  other  structural  properties  (the 
"path"  of  a structural  transformation),  and  variations  of  structural 
properties  with  time  (the  "kinetics"  of  a structural  transformation). 
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Many  examples  of  application  of  the  concepts  of  structology  are 
given  by  Rhines  [365].  The  following  are  some  few  examples  of 
phenomena  that  can  be  associated  with  structural  properties: 

1)  Grain  growth  [523]. 

2)  Precipitation  of  a second  phase  [524]. 

3)  Sintering  [525,526]. 

4)  High  temperature  oxidation  [527,528]. 

5)  Grain  boundary  hardening  [529-531]. 

6)  Tensile  strength  [532]. 

The  application  of  structology  involves  the  following  steps: 

Step  1 - Measurement  of  structural  properties. 

Step  2 - Measurement  of  physical  properties. 

Step  3 - Description  of  relationships  between  the  structural  and 
physical  properties,  with  the  purpose  of  discovering 
regularities  of  behavior. 

Step  4 - Explanation  of  the  mechanism  (or  mechanisms)  that  lead  to 
the  structure-property  relationships  observed. 

The  objective  of  this  research  is  to  accomplish  Steps  1,  2 and  3. 
However,  possible  explanations  for  the  observed  structure-property 
relationships  will  be  given  whenever  possible. 

6.2  Structure  and  Mechanical  Properties 

The  variation  of  Vickers  hardness  number  (VHN)  as  a function  of 
volume  fraction  of  pores  is  shown  in  Figure  6.1,  for  a set  of  silica  gel 
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Figure  6.1  Variation  of  Vickers  hardness  number  (VHN)  as  a 
function  of  volume  fraction  of  pores  (Vv). 
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samples  of  12  A average  pore  size.  The  figure  suggests  that  smaller  a 
Vv  of  pores  are  associated  with  increasing  hardness  of  the  gel-silica 
samples.  Figure  6.2  shows  the  variation  of  VHN  as  a function  of  the 
topological  index  6.  It  can  be  noted  that  the  relationship  between  VHN 
and  6 is  not  as  linear  as  the  relationship  between  VHN  and  Vv  of  pores. 
That  fact  suggests  that  the  hardness  of  the  gel-silica  samples  depends 
more  on  the  level  of  porosity  of  the  specimen  than  on  the  number  of 
pores. 

The  connectivity  of  the  pores  influences  the  hardness  in  the  manner 
shown  in  Figure  6.3.  The  hardness  is  nearly  independent  of  pore 
connectivity  up  to  a value  of  relative  genus  of  0.2.  When  the  relative 
genus  (y)  becomes  smaller  than  0.2  the  microhardness  increases 
considerably.  During  the  third  stage  of  densification,  the  relative  genus 
is  constant  at  zero  and  the  hardness  continues  to  increase  from  300  to 
770  VHN.  This  can  be  seen  more  clearly  in  Figure  6.4,  which  shows 
the  variation  of  VHN  as  a function  of  the  pore  coordination  number 
(CN).  The  figure  shows  that  during  the  second  stage  of  densification 
(CN  varying  from  4 to  2)  VHN  does  not  change  as  dramatically  as  when 
it  does  during  the  third  stage  (CN  varying  from  2 to  0).  This 
observation  suggests  that  pore  closing  promotes  stiffening  of  the  gel- 
silica  structure  to  a very  large  degree. 

The  variation  of  flexural  strength  of  silica  gel  monoliths  as  a 
function  of  volume  fraction  of  pores  is  shown  in  Figure  6.5.  This 
figure  suggests  a fairly  good  correlation  between  flexural  strength  and 
Vv.  However,  Figures  6.6  and  6.7  suggest  slightly  better  correlations 
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Figure  6.2  Variation  of  Vickers  hardness  number  (VHN)  as  a 
function  of  topological  index  6. 
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Figure  6.3  Variation  of  Vickers  hardness  number  (VHN)  as  a 
function  of  relative  genus  (y). 
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Figure  6.4  Variation  of  the  Vickers  hardness  number  (VHN)  as  a 
function  of  the  pore  coordination  number  (CN). 
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Figure  6.5  Variation  of  flexural  strength  as  a function  of 
volume  fraction  of  pores  (Vy). 
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Figure  6.6  Variation  of  the  flexural  strength  as  a function  of 

topological  index  6. 


308 


RELATIVE  GENUS 


Figure  6.7  Variation  of  flexural  strength  as  a function  of 

relative  genus  (y). 
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of  flexural  strength  with  the  topological  parameters  (topological  index  6 
and  relative  genus  y).  Although  it  is  not  conclusive  that  6 and  y 
influences  flexural  strength  more  than  the  porosity,  Figures  6.5,  6.6  and 
6.7  show  that  the  topological  parameters  derived  are  suitable  for  use  in 
the  study  of  structure-properties  relationships. 

The  flexural  modulus  varies  with  volume  fraction  of  pores  as  shown 
in  Figure  6.8.  Again  the  topological  parameters  6 and  y seem  to  yield 
better  relationships  with  flexural  modulus,  as  shown  in  Figures  6.9  (6) 
and  6.10  (y). 

Figure  6.11  shows  the  variation  of  compression  strength  of  the  silica 
gel-monoliths  (for  samples  type  Al)  as  a function  of  volume  fraction  of 
pores.  As  shown  in  the  figure,  a linear  relationship  between  Vv  and 
compression  strength  does  not  exist  in  the  experimental  range  studied, 
and  also  Vv  varies  within  a very  small  range.  A linear  expression  is 
derived  just  as  a way  to  permit  comparison  between  the  applicability  of 
Vv,  6 and  y.  Figure  6.12  shows  the  variation  of  compression  strength 
as  a function  of  topological  index  6.  This  figure  suggests  that  the 
compression  strength  correlates  significantly  better  with  variations  in 
topological  index  6 than  with  variations  in  Vv.  The  correlation  of 
compression  strength  with  the  pore  connectivity  seems  to  be  the  best  for 
the  experimental  range  studied,  as  shown  by  the  variation  of 
compression  strength  with  relative  genus  in  Figure  6.13. 

The  variation  of  compression  modulus  (for  samples  type  Al)  with 
volume  fraction  of  pores  is  shown  in  Figure  6.14. 
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Figure  6.8  Variation  of  flexural  modulus  as  a function  of 
volume  fraction  of  pores  (Vv)- 
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Figure  6.9  Variation  of  flexural  modulus  as  a function  of 

topological  index  6. 
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Figure  6.10  Variation  of  flexural  modulus  as  a function  of 

relative  genus  (y ). 
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Figure  6.11  Variation  of  compression  strength  as  a function  of 
volume  fraction  of  pores  (Vv)  (Samples  type  Al). 
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Figure  6.12  Variation  of  compression  strength  as  a function  of 
topological  index  6 (Samples  type  Al). 
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Figure  6.13  Variation  of  compression  strength  as  a function  of 
relative  genus  (y)  (Samples  type  Al). 
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Figure  6.14  Variation  of  compression  modulus  as  a function  of  volume 
fraction  of  pores  (Vv)  (Samples  type  Al). 
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As  the  densification  takes  place,  Vv  changes  and  the  compression 
modulus  changes  relatively  fast.  Figure  6.14  suggests  that  for  the  later 
stages  of  densification  (smaller  values  of  Vv),  variations  in  Vv  produces 
decreasing  effects  in  the  compression  modulus.  The  variation  of 
compression  modulus  with  topological  index  B is  shown  in  Figure  6.15. 
This  figure  suggests  a better  correlation  between  compression  modulus 
and  B when  compared  to  the  relationship  between  compression  modulus 
and  Vy.  The  correlation  between  relative  genus  and  compression 
modulus  seems  to  be  the  best  for  the  experimental  range  studied,  as  seen 
in  Figure  6.16. 

The  variation  of  compression  modulus  for  a second  set  of  samples 
(A2)  as  a function  of  Vy,  B and  y is  shown  in  Figures  6.17,  6.18  and 
6.19  respectively.  For  this  set  of  experiments  again  the  topological 
parameters  seem  to  describe  better  than  Vv  the  physical  evolution  of  the 
gel-silica  samples.  However,  for  this  case  the  topological  index  B, 
rather  than  the  relative  genus,  yields  a better  relationship. 

Figure  6.20  shows  the  variation  of  tensile  strength  (as  obtained  with 
the  Brazilian  test,  for  samples  type  Al)  as  a function  of  volume  fraction 
of  pores.  The  linear  correspondence  between  tensile  strength  and  Vv  is 
poor,  as  shown  by  the  low  coefficient  of  correlation  (R2  = 0.749)  and 
the  presence  of  a relatively  isolated  data  point  (at  Vv  of  around  0.40) 
which  influences  dramatically  the  construction  of  the  least  squares 
straight  line.  In  this  case  the  whole  range  of  variation  of  Vy  is  very 
small.  The  variation  of  tensile  strength  with  topological  index  S is 
shown  in  Figure  6.21.  The  structure-property  relationship  correlates 
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Figure  6.15  Variation  of  compression  modulus  as  a function  of 
topological  index  8 (Samples  type  Al). 
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Figure  6.16  Variation  of  compression  modulus  as  a function  of 
relative  genus  (y)  (Samples  type  Al). 


320 


Figure  6.17  Variation  of  compression  modulus  as  a function  of 
volume  fraction  of  pores  (Vv)  (Samples  type  A2). 
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Figure  6.18  Variation  of  compression  modulus  as  a function  of 
topological  index  6 (Samples  type  A2). 
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Figure  6.19  Variation  of  compression  modulus  as  a function  of 
relative  genus  (y)  (Samples  type  A2). 
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Figure  6.20  Variation  of  tensile  strength  (Brazilian  test)  as  a function 
of  volume  fraction  of  pores  (Vv)  (Samples  type  Al). 
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Figure  6.21  Variation  of  tensile  strength  (Brazilian  test)  as  a function 
of  topological  index  6 (Samples  type  Al). 
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better  with  6 than  Vv.  The  relationship  between  tensile  strength  and 
relative  genus  (y)  is  shown  in  Figure  6.22.  For  this  set  of  samples  (Al) 
both  6 and  y yield  better  relationships  with  tensile  strength  than  Vv 
does. 

The  relationships  between  tensile  strength  with  Vv,  B and  y for  a 
second  set  of  samples  (A2)  are  shown  in  Figures  6.23,  6.24  and  6.25 
respectively.  In  the  case  of  this  set  of  samples,  the  tensile  strength  has  a 
better  linear  correspondence  with  Vv  than  with  the  topological 
parameters  B and  y.  It  should  be  the  noted  however,  that  the  linear 
equation  representing  the  variation  of  tensile  strength  with  Vv  is 
markedly  influenced  by  the  presence  of  an  isolated  data  point  (at  65 
MPa,  as  shown  in  Figure  6.23).  The  variation  of  tensile  strength  with 
relative  genus  shown  in  Figure  6.25  suggests  again  a sensitive  effect  of 
pore  closing  on  the  physical  strength  of  gel-silica  monoliths.  The  effect 
of  pore  closing  on  tensile  strength  is  similar  to  the  effect  of  pore  closing 
on  the  microhardness,  as  seen  earlier  (Figure  6.3).  In  the  case  of 
Figures  6.23,  6.24  and  6.25  a third  order  polynomial  was  chosen  to 
provide  a better  fit  for  the  experimental  points. 

The  variation  of  tensile  modulus  (DMA)  as  a function  of  volume 
fraction  of  pores  is  shown  in  Figure  6.26.  The  figure  suggests  a fairly 
good  correlation  between  tensile  modulus  and  Vv.  The  variation  of 
tensile  modulus  with  surface  area  of  pores  per  unit  volume  (Sv)  is 
shown  in  Figure  6.27.  The  relationship  of  tensile  modulus  with  Sy 
seems  to  be  better  than  with  Vv,  as  indicated  by  the  linear  coefficient  of 
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Figure  6.22  Variation  of  tensile  strength  (Brazilian  test)  as  a function 
of  relative  genus  (y)  (Samples  type  Al). 
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Figure  6.23  Variation  of  tensile  strength  (Brazilian  test)  as  a function 
of  volume  fraction  of  pores  (Vv)  (Samples  type  A2). 
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Figure  6.24  Variation  of  tensile  strength  (Brazilian  test)  as  a function 
of  topological  index  6 (Samples  type  A2). 
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Figure  6.25  Variation  of  tensile  strength  (Brazilian  test)  as  a function 
of  relative  genus  (y)  (Samples  type  A2). 
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Vv  pores 


Figure  6.26  Variation  of  tensile  modulus  (DMA)  as  a function  of 
volume  fraction  of  pores  (Vv). 
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Sv  pores  (10A6  cm-1) 


Figure  6.27  Variation  of  tensile  modulus  (DMA)  as  a function  of 
surface  area  per  unit  volume  (Sv). 
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correlation  (R2).  The  variation  of  tensile  modulus  with  the  topological 
parameters  6 and  g are  shown  in  Figures  6.28  and  6.29  respectively. 
Again,  the  topological  parameters  seem  to  correlate  better  to  the 
physical  property  than  does  the  volume  fraction  of  pores  (Vv). 

Figures  6.30  and  6.31  present  the  variation  of  tensile  modulus 
(DMA)  measured  in  the  range  of  25  °C  (298  K)  to  250  °C  (523  K)  as  a 
function  of  Vv  and  B respectively.  The  linear  relationship  between 
tensile  modulus  and  the  structural  parameters  is  very  poor  for  higher 
testing  temperatures.  However,  the  decrease  in  the  gel-silica  tensile 
modulus  as  the  testing  temperature  increases  can  be  clearly  seen  in  both 
Figures  6.30  and  6.31. 

6.3  Structure  and  Thermal  Properties 

The  variation  of  dimensional  change  at  690  °C  (963  K)  as  a function 
of  volume  fraction  of  pores  is  shown  in  Figure  6.32.  At  high  volume 
fraction  of  pores  small  variations  in  Vv  lead  to  large  variations  in  the 
dimensional  change  of  gel-silica.  For  smaller  values  of  Vv  the  gel-silica 
specimens  are  more  thermally  stable,  as  shown  in  Figure  6.32  by  the 
relatively  small  variations  in  dimensional  change.  Thus,  dimensional 
change  at  690  °C  (963  K)  dependence  on  Vv  is  best  represented  as  the 
intersection  of  two  straight  lines. 

The  variation  of  dimensional  change  at  690  °C  (963  K)  with 
topological  index  B is  similar  to  that  presented  by  Vv?  i.e.,  an 
intersection  of  two  linear  relations  (Figure  6.33).  For  the  beginning  of 
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Figure  6.28  Variation  of  tensile  modulus  (DMA)  as  a function  of 

topological  index  6. 
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Figure  6.29  Variation  of  tensile  modulus  (DMA)  as  a function  of 

relative  genus  (y). 
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Vv  pores 


Figure  6.30  Variation  of  tensile  modulus  (DMA)  as  a function  of 
volume  fraction  of  pores  (Vv)  for  testing  temperatures  varying 

from  25  °C  to  250  °C. 
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Figure  6.31  Variation  of  tensile  modulus  (DMA)  as  a function  of 
topological  index  6 for  testing  temperatures  varying  from 

25  °C  to  250  °C. 
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Figure  6.32  Variation  of  dimensional  change  at  690  °C  as  a 
function  of  volume  fraction  of  pores  (Vv). 
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Figure  6.33  Variation  of  dimensional  change  at  690  °C  as  a 
function  of  topological  index  6. 
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densification  (6  varying  from  0 to  0.2)  the  silica  gel  dimension  changes 
significantly.  As  6 increases  from  0 to  0.2,  the  dimensional  change  at 
690  °C  decreases  from  -60x10-3  cm/cm  to  less  than  -1x10-3  cm/cm. 
That  change  illustrates  that  a sample  with  higher  6 will  shrink  less  upon 
heating  than  a sample  with  a lower  6.  For  values  of  6 between  0.2  and 
1.0  the  dimensional  change  variations  are  relatively  small  as  compared 
to  samples  in  earlier  stages  of  densification. 

A similar  effect  is  observed  when  one  examines  the  dependence  of 
dimensional  change  at  690  °C  (963  K)  on  the  degree  of  connectivity  of 
the  gel-silica  structure.  Figure  6.34  shows  the  variation  of  dimensional 
change  at  690  °C  (963  K)  as  a function  of  relative  genus.  For  the  dried 
sample  (relative  genus  of  1)  there  is  a large  shrinkage  when  the 
specimen  is  heated  to  690  °C  (963  K).  As  the  gel-silica  densifies 
(relative  genus  going  from  1 to  0)  the  sample  becomes  increasingly 
more  thermally  stable. 

The  observation  that  the  silica  gel  monoliths  increase  in  thermal 
stability  as  the  densification  takes  place  is  valid,  but  also  is  arbitrary 
because  the  dimensional  change  at  an  arbitrary  temperature  (690  °C, 
963  K)  was  chosen  as  reference.  As  shown  in  Figure  6.35,  the  average 
coefficient  of  dimensional  change  (CDCm)  assumes  different  values  for 
different  reference  temperatures.  The  variation  of  CDCm  as  a function 
of  densification  temperature,  assuming  reference  temperatures  in  the 
range  30  °C  (303  K)  to  690  °C  (963  K),  is  plotted  in  Figure  6.35.  As 
shown  by  Figure  6.35,  the  curve  of  variation  of  CDCm  with 
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Figure  6.34  Variation  of  dimensional  change  at  690  °C  as  a 
function  of  relative  genus  (y). 
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30  °C 
100  °C 
200  °C 
300  °C 
400  °C 
500  °C 
600  °C 
690  °C 


Figure  6.35  Variation  of  the  average  coefficient  of  dimensional  change 
(CDCm)  as  a function  of  densification  temperature,  for  reference 
temperatures  varying  from  30  °C  to  690  °C. 
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densification  temperature  depends  on  the  assumed  reference 
temperature. 

Figure  6.36  shows  the  variation  of  dimensional  change  as  a function 
of  testing  temperature,  for  samples  densified  in  the  range  of  500  °C 
(773  K)  to  1150  °C  (1423  K).  The  range  of  testing  temperature 
indicated  in  the  figure  is  from  200  °C  (473  K)  to  460  °C  (733  K).  The 
lower  limit  of  200  °C  (473  K)  comes  from  the  choosing  of  a 
temperature  higher  than  Tzero  (for  the  sample  densified  at  500  °C,  773 
K).  The  upper  limit  of  460  °C  was  chosen  to  be  smaller  than  the 
temperature  of  the  sample  densified  at  the  lower  temperature  (again  500 
°C,  773  K).  Therefore,  for  the  set  of  samples  densified  in  the  range  500 
°C  (773  K)  to  1150  °C  (1423  K)  the  testing  temperature  range  of  200  °C 
(473  K)  to  460  °C  733  K)  should  provide  an  unbiased  region  of 
dimensional  change  for  comparisons  between  the  thermal  behavior  of 
that  set  of  samples.  As  seen  in  Figure  6.36,  the  curves  of  variation  of 
dimensional  change  with  testing  temperature  for  all  silica  gel  samples 
present  good  linearity,  as  shown  by  the  linear  coefficients  of 
correlation. 

Plotting  the  slopes  of  the  dimensional  change  versus  temperature 
curves  (shown  in  Figure  6.36),  as  a function  of  densification 
temperature  results  in  Figure  6.37.  According  to  Figure  6.37  there 
seems  to  be  a good  linearity  between  the  dimensional  change  slope  (or 
slope  between  the  dimensional  change  and  testing  temperature  curve 
measured  between  200  °C  (473  K)  and  460  °C  (733  K))  and  the 
densification  temperature.  The  dimensional  change  slope  (DCS)  can 
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Figure  6.36  Variation  of  dimensional  change  as  a function  of  testing 
temperature  for  samples  densified  in  the  range  of  500  °C  to  1150  °C. 
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Figure  6.37  Variation  of  dimensional  change  slope  (DCS)  as  a 
function  of  densification  temperature. 
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thus  be  used  to  characterize  the  dimensional  change  properties  of  gel- 
silica  samples  densified  in  the  range  of  500  °C  (773  K)  to  1150  °C  (1423 
K).  It  is  important  to  point  out  that  the  notation  "dimensional  change 
slope"  was  used  instead  of  the  coefficient  of  dimensional  change 
designation.  The  reason  for  that  is  because  the  dimensional  change 
slope  (DCs),  as  defined  and  used  in  this  work,  differs  from  both  the 
average  coefficient  of  dimensional  change  (which  uses  the  original 
length  of  the  sample  as  reference)  and  the  instantaneous  coefficient  of 
dimensional  change  (evaluated  as  a limit  for  each  point  of  the  curve  of 
dimensional  change  versus  temperature).  It  should  be  noted  however 
that  since  all  curves  shown  in  Figure  6.36  are  quite  linear,  the 
instantaneous  coefficient  of  dimensional  change  measured  at  a 
convenient  temperature  inside  the  range  200  °C  (473  K)  to  460  °C  (733 
K)  (say  300  °C,  573  K)  should  provide  a good  approximation  to  the 
"dimensional  change  slope"  used  in  this  work. 

The  variation  in  the  average  coefficient  of  dimensional  change 
(CDCm)  as  a function  of  temperature  is  shown  in  Figure  6.38.  The 
range  of  testing  temperature  (between  340  °C  (613  K)  and  460  °C  (733 
K))  was  chosen  so  that  linear  variations  of  the  coefficient  of  dimensional 
change  with  temperature  would  occur  inside  the  limiting  temperatures. 
In  the  case  of  Figure  6.38  it  was  introduced  data  for  dried  gel-silica. 
From  the  curves  of  CDCm  versus  temperature  (Figure  6.38)  the  slopes 
were  calculated,  as  shown  in  Figure  6.38.  The  variation  of  the 
coefficient  of  dimensional  change  slope  (CDCS)  as  a function  of 
densification  temperature  is  plotted  in  Figure  6.39.  This  figure  shows 
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Figure  6.38  Variation  of  the  average  coefficient  of  dimensional  change 
(CDCm)  as  a function  of  testing  temperature  for  silica  gel  samples 
densified  in  the  range  180  °C  to  1150  °C. 
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Figure  6.39  Variation  of  coefficient  of  dimensional  change  slope 
(CDCS)  as  a function  of  densification  temperature. 
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that  larger  variations  in  the  coefficient  of  dimensional  change  with 
densification  temperature  are  present  for  samples  densified  at  lower 
temperatures.  Densification  at  650  °C  (923  K)  and  above  yields  a stable 
CDC  behavior,  within  the  range  of  200-460  °C  (473-733  K). 

The  variation  of  the  dimensional  change  slope  (DCS)  as  a function  of 
volume  fraction  of  pores  is  shown  in  Figure  6.40.  It  is  seem  from  this 
figure  that  in  the  beginning  of  densification  small  variations  in  Vv  cause 
relatively  large  dimensional  variations  in  silica  gel  monoliths.  For  a 
large  variation  in  Vv  (from  0.4  to  0.0)  there  are  relatively  small 
variations  in  dimensional  change.  Therefore,  the  dimensional  change 
behavior  of  gel-silica  does  not  have  a direct  relationship  with  Vv,  but 
appears  to  be  the  sum  of  two  relationships.  The  variation  of 
dimensional  change  slope  (DCS)  with  6 is  shown  in  Figure  6.41.  In  this 
case,  6 provides  a better  linear  correlation  (not  shown  in  Figure  6.41) 
with  dimensional  change  slope  than  Vv  does.  However  the  relationship 
between  pore  connectivity  (as  measured  by  y)  and  dimensional  change 
slope  is  the  best,  as  seen  in  Figure  6.42.  As  densification  proceeds 
(from  y=l  to  y^0)  there  is  a direct  correspondence  between  decrease  in 
pore  connectivity  and  dimensional  change.  For  the  third  stage  of 
densification  (when  y is  constant  at  zero)  the  variation  in  the  number  of 
pores  present  in  the  structure  (as  expressed  by  B)  can  be  associated  with 
variations  in  the  dimensional  change  behavior  of  gel-silica. 
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Figure  6.40  Variation  of  dimensional  change  slope  (DCS)  as  a 
function  of  volume  fraction  of  pores  (Vv). 
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Figure  6.41  Variation  of  dimensional  change  slope  (DCS)  as  a 
function  of  topological  index  6. 
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Figure  6.42  Variation  of  dimensional  change  slope  (DCS)  as  a 
function  of  relative  genus  (7). 
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6.4  Applications 

As  showed  previously  both  metric  and  topological  parameters  can  be 
used  to  describe  the  evolution  of  the  gel-silica  structure  during 
sintering.  The  relationships  between  structural  parameters  and  physical 
properties  obtained  suggest  that  it  may  be  possible  to  predict  physical 
properties  based  on  knowledge  of  the  structural  state  for  sets  of  samples 
processed  under  similar  conditions. 

The  applicability  of  either  metric  or  topological  parameters  is 
function  of  specific  conditions  of  the  problem.  Physical  properties  can 
be  correlated  with  both  metric  and  topological  parameters.  The  amount 
of  void  spaces  in  a structure  relates  with  mechanical  strength  since  pores 
decrease  the  actual  cross  sectional  area  that  supports  the  physical  force 
and  pores  also  act  as  stress  raisers.  The  connectivity  on  the  other  hand 
can  be  used  to  differentiate  mechanical  states  associated  with  pore 
closing. 

For  catalytic  applications  of  gel-silica,  knowledge  of  the  surface  area 
of  pores  per  unit  volume  is  of  paramount  importance.  However,  large 
surface  areas  in  a structure  of  closed  pores  would  not  result  in  large 
catalytic  activity.  The  connectivity  in  this  case  is  associated  with  the 
number  of  paths  the  substances  can  take  in  order  for  the  catalytic 
reactions  occur. 

In  doping  and  impregnation  applications  of  porous  gel-silica  the 
volume  fraction  of  pores  is  a structural  parameter  of  obvious 
importance.  However,  in  a closed  pore  structure  no  doping  or 


353 


impregnation  can  be  taken  place  (except  through  slower  solid  diffusion 
processes).  Therefore,  the  connectivity  of  the  structure  should  be 
known  as  requisite  for  doping  and  impregnation  applications.  Figure 
6.43  shows  the  variation  of  helium  flow  as  a function  of  gas  pressure 
for  samples  densified  in  the  range  of  180  °C  (453  K)  to  1100  °C  (1373 
K).  This  experiment  was  carried  out  using  samples  of  1.9  mm  thick 
attached  to  a sample  holder  directly  connected  to  a helium  gas  cylinder. 
The  gas  pressure  was  measured  directly  by  a regulator  attached  to  the 
gas  cylinder.  The  gas  leaving  the  gas  cylinder  (at  room  temperature)  is 
forced  (by  its  own  pressure)  through  the  sample.  The  gas  flow  is 
measured  by  a bubble  flowmeter  connected  to  the  sample  holder.  The 
experiment  for  each  sample  begins  with  a low  pressure  and  for  each 
pressure  the  gas  flow  through  the  sample  is  measured.  Figure  6.43 
indicates  a larger  flow  of  helium  through  samples  densified  at  lower 
temperatures.  For  the  sample  densified  at  1100  °C  no  gas  flow  was 
observed.  A sample  densified  at  1150  °C  (1423  K)  also  presented  no 
gas  flow  in  the  experiment.  The  slope  of  the  curves  in  Figure  6.43  are 
plotted  as  a function  of  densification  temperature  in  Figure  6.44.  The 
variation  of  the  slope  of  helium  flow  rate  as  a function  of  volume 
fraction  of  pores  is  presented  in  Figure  6.45.  The  variation  of  the  slope 
of  helium  flow  rate  with  relative  genus  is  shown  in  Figure  6.46. 
Comparing  Figures  6.45  and  6.46  and  the  shape  of  Figure  6.44  with 
Figures  2.27  (Vv  versus  temperature)  and  3.21  (relative  genus  versus 
temperature)  it  can  be  seen  that  the  gas  flow  through  gel-silica 
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Figure  6.43  Variation  of  helium  flow  as  a function  of  gas  pressure  for 
samples  densified  in  the  range  180  °C  to  1100  °C. 
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Figure  6.44  Variation  of  the  slope  of  the  curves  helium  flow  versus 
pressure  (Figure  6.43)  as  a function  of  densification  temperature. 
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Figure  6.45  Variation  of  the  helium  flow  slope  as  a function  of 
volume  fraction  of  pores  (Vv). 


HELIUM  FLOW  SLOPE 


357 


Figure  6.46  Variation  of  the  helium  flow  slope  as  a function  of 

relative  genus  (7). 
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monoliths  correlates  better  with  the  pore  connectivity  than  with 
porosity. 


6.5  Structological  Analysis 

Structures  are  usually  complex  and  involve  many  different  geometric 
aspects.  Therefore,  it  should  not  be  expected  that  simplistic  approaches 
to  the  problem  of  correlating  structures  with  physical  properties  are 
likely  to  yield  meaningful  relationships.  The  more  complex  a structure 
is  the  larger  the  number  of  structural  parameters  required  to  describe 
it.  The  introduction  of  topological  parameters  in  describing  gel-silica 
structures  adds  to  the  understanding  of  how  a particular  structure 
correlates  to  its  physical  properties.  Table  6.1  summarizes  the  results 
obtained  for  the  structure-properties  relationships  derived  in  this  work. 
The  use  of  topological  parameters  should  be  intended  to  complement 
information  given  by  metric  parameters  and  not  to  replace  then. 

Considering  for  instance  the  problem  of  impregnating  the  porous 
gel-silica  with  a certain  type  of  molecules,  many  structural  features  have 
to  be  brought  up.  First,  one  has  to  consider  whether  the  volume 
fraction  of  pores  is  compatible  or  not  with  the  amount  of  impregnating 
material  intended  to  be  introduced  in  the  matrix.  If  the  porosity  is 
adequate  that  is  not  a guarantee  that  the  impregnation  procedure  is 
going  to  be  successful.  If  the  porosity  is  adequate  but  the  pores  are 
closed,  there  will  be  no  impregnation.  Therefore,  the  knowledge  of  the 
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Table  6.1  Coefficients  of  correlation  (R2)  between  physical 
properties  and  metric  (Vv)  and  topological  parameters 
(6  and  y).  A (3)  after  the  number  means  that  R2  was  taken 
out  of  a third  order  polynomial.  No  number  indicates  a 
first  order  polynomial  (straight  line). 


TEST 

R2 

Vv 

R2 

B 

R2 

Y 

VHN 

0.963 

0.991  (3) 

0.937  (3) 

Flexural  Strength 

0.954 

0.994 

0.994 

Flexural  Modulus 

0.809 

0.907 

0.907 

Compression  Strength  (Al) 

0.555 

0.860 

0.899 

Compression  Modulus  (Al) 

0.337 

0.731 

0.796 

Compression  Modulus  (A2) 

0.873 

0.910 

0.895 

Tensile  Strength  (Al) 

0.749 

0.881 

0.880 

Tensile  Strength  (A2) 

0.997  (3) 

0.995  (3) 

0.982  (3) 

Tensile  Modulus 

0.910 

0.967 

0.967 

Dimensional  Change 

0.695 

0.863 

0.957 

Helium  Flow 

0.395 

0.451 

0.711 
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pore  interconnectivity  is  a must  for  controlled  impregnation  and 
doping  procedures.  Besides  the  porosity  and  connectivity,  the  size  of 
the  pores  also  plays  a decisive  role  during  doping  and  impregnation.  If 
the  porosity  and  connectivity  of  the  sample  are  adequate  but  the  size  of 
the  pores  are  too  small  for  the  impregnating  molecules  to  go  through, 
again  the  impregnation  will  not  occur.  Therefore,  a combination  of 
structural  parameters  (metric  and  topological)  are  necessary  to  create  a 
more  complete  picture  of  a given  practical  situation. 

The  VHN  data  obtained  for  silica  gel  monoliths  correlate  better  with 
the  volume  fraction  of  pores  than  with  the  connectivity  of  pores.  This 
result  is  possibly  related  to  the  fact  that  a hardness  test  measures  directly 
the  resistance  of  the  material  against  the  penetration  of  the  diamond 
indenter.  That  means  that  the  resistance  to  the  penetration  is  directly 
proportional  to  the  cross  sectional  area  that  is  opposing  the  diamond 
indenter;  i.e.,  the  actual  area  fraction  of  material  (Figure  6.47).  Since 
the  areal  fraction  is  equal  to  the  volume  fraction,  this  parameter  (Vv)  is 
appropriate  to  describe  the  expected  hardness  of  a material. 

The  data  on  strength  and  modulus  of  silica  gel  monoliths  correlate 
better  with  the  topological  parameters  (6  and  y)  than  with  the  volume 
fraction  of  pores.  It  is  possible  that  these  results  come  from  a higher 
concentration  of  stresses  in  a pore  with  several  branches  than  in  a pore 
with  one  or  no  branch  (Figure  6.48).  It  was  shown  that  the  mechanical 
properties  variations  with  densification,  can  be  divided  into  two  distinct 
regions:  the  first  one  is  at  lower  temperatures,  when  the  strength  of  the 
silica  gel  monoliths  increases  relatively  slowly  with  densification 
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indentation 


pores 


Figure  6.47  Schematic  of  the  actual  area  in  a specimen  that  resists  the 
penetration  of  a hardness  indenter. 
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Figure  6.48  Schematic  of  the  reduction  in  concentration  of  stresses 
when  the  pore  connectivity  decreases,  (a)  CN  = 4;  (b)  CN  = 2; 

(c)  CN  = 0. 


363 


temperature;  the  second  region  is  at  higher  temperatures,  when  there  is 
a more  pronounced  increase  in  strength  as  densification  proceeds.  In 
the  first  region  the  pores  have  many  branches,  which  act  as  stress  risers, 
decreasing  the  strength  of  the  material.  When  the  pores  start  closing  in 
the  second  region  (during  the  third  stage  of  densification),  the  strength 
increases  dramatically,  because  of  the  decrease  in  concentration  of  stress 
as  the  branches  and  nodes  disappear. 

As  observed  in  this  research,  the  dimensional  change  correlates 
better  with  the  pore  connectivity  than  with  volume  fraction  of  pores.  It 
is  possible  that  this  effect  be  related  with  the  reasoning  that  follows. 
When  Vv  of  pores  vary,  but  the  solid  matrix  is  continuous  (Figure  6.49) 
it  is  expected  that  the  dimensional  change  of  the  solid  matrix  will  not  be 
affected.  By  other  hand,  when  the  pore  connectivity  changes  (Figure 
6.50),  even  though  the  matrix  is  still  connected  in  the  three-dimensional 
space,  the  branched  pores  change  the  thermal  behavior  of  the  material. 
These  changes  could  come  from  interferences  in  the  heat  flow  through 
the  sample  caused  by  movement  of  vapor  inside  the  connected  pores, 
and  from  interference  in  the  overall  dimensional  change  caused  by 
expansion  (or  contraction)  of  the  connected  pores. 

The  flow  of  gas  through  the  silica  gel  monoliths  correlates  better 
with  the  pore  connectivity  than  with  Vv.  This  result  is  expected,  since  a 
gas  can  flow  through  a connected  pore,  and  the  flow  stops  when  the 
pores  are  closed. 

The  introduction  of  topological  characterization  to  the  study  of  gel- 
silica  monoliths  opens  up  opportunities  for  better  understanding  of  the 
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Figure  6.49  The  schematic  structure  shown  in  (a)  has  a larger  Vv  of 
pores  than  the  structure  shown  in  (b),  however  both  structures  have 

a continuous  matrix. 
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Figure  6.50  The  schematic  structure  shown  in  (a)  has  a larger  pore 
connectivity  than  the  structure  shown  in  (b). 
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problem  of  structure-properties  correlations.  A more  comprehensive 
description  of  the  structure  is  more  likely  to  lead  to  more  accurate 
relationships  with  physical  properties.  The  ultimate  goal  of  such  a type 
of  structural  description  is  to  make  feasible  the  designing  and  obtaining 
of  sets  of  desired  properties  in  a material  through  controlled  processing. 


CHAPTER  7 

CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  WORK 


The  objectives  of  this  research  were  to  study  the  evolution  of 
structural  and  physical  properties  during  densification  of  silica  gel 
monoliths  and  to  provide  new  knowledge  on  the  relationships  between 
structure  and  properties.  The  main  conclusions  of  this  research  can  be 
summarized  as  follows: 

1)  Nitrogen  gas  adsorption  experiments  show  that  as  the 
densification  of  silica  gels  proceeds,  the  volume  fraction  of  pores  (Vv) 
and  the  surface  area  of  pores  per  unit  volume  (Sv)  decrease  towards 
zero  porosity  for  the  fully  dense  material.  Meanwhile,  the  average  size 
of  pores  remains  approximately  constant  during  densification  of  silica 
gel  monoliths  with  pores  of  12  A initial  average  pore  size.  For  silica 
gels  of  larger  initial  pore  sizes,  the  average  pore  size  slightly  decreases 
during  the  beginning  of  densification. 

2)  As  the  densification  proceeds,  there  is  a significant  variation  in 
the  structural  density  of  the  silica  gels.  The  structural  density  initially 
increases  from  the  dried  samples  and  assumes  maximum  values  in  the 
temperature  range  of  800-1000  °C  (1073-1273  K),  depending  on  the 
average  pore  size  of  the  gels.  Increasing  the  densification  temperature 
beyond  that  maximum  the  structural  density  falls  off  again. 

3)  The  variation  in  structural  metric  parameters  for  silica  gels  shows 
that  the  structural  evolution  during  densification  can  be  roughly  divided 
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into  two  major  regions,  which  depend  on  the  pore  size  of  the  silica  gel. 
For  silica  gels  of  12  A average  pore  size  the  first  region  goes  from  180 
°C  (453  K)  (dried  sample)  until  a temperature  of  about  800  °C  (1073 
K).  In  this  region  the  structural  transformations  with  densification 
temperature  are  less  pronounced  than  for  the  second  region,  which  goes 
from  about  800  °C  (1073  K)  until  the  sample  is  fully  dense  (at  about 
1150  °C,  1423  K). 

4)  It  is  possible  to  use  model  geometric  shapes  together  with 
experimentally  obtained  values  of  Vv  and  Sv  to  obtain  topological 
information  on  a structure.  Three  types  of  models  (spherical, 
tetrahedral  and  cylindrical)  were  used  to  describe  the  evolution  of 
topological  properties  of  silica  gels  during  densification.  As  the 
densification  proceeds,  the  number  of  nodes  and  branches  associated 
with  the  gel  structure  decreases.  The  genus  decreases  during 
densification  and  goes  to  zero  at  a temperature  which  depends  on  the 
pore  size.  The  region  where  the  genus  becomes  zero  characterizes  the 
beginning  of  the  third  stage  of  densification.  The  pore  closing 
associated  with  the  third  stage  of  densification  is  responsible  for  water 
trapping  and  foaming  when  water  is  not  removed  before  the  genus 
becomes  zero. 

5)  In  order  to  make  it  easier  to  compare  the  structural  state  and 
characterize  the  structural  evolution  during  sintering  two  topological 
indices  are  used:  (1)  the  relative  genus  y and  (2)  the  fraction  of  removed 
branches  6.  The  relative  genus  is  used  to  identify  the  different  stages  of 
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densification.  The  topological  index  6 is  associated  to  pore  network 
transformations  during  the  entire  process. 

6)  In  general  the  mechanical  property  variations  during  densification 
for  silica  gels  samples  show  the  presence  of  two  distinct  regions:  at  low 
temperatures  the  variations  in  the  mechanical  properties  are  less 
pronounced  than  in  a second  region  corresponding  to  higher  sintering 
temperatures. 

7)  Comparisons  between  the  thermal  properties  behavior  of  silica 
gels  are  carried  out  by  using  the  coefficient  of  dimensional  change 
measured  in  a temperature  range  that  lies  at  the  same  time  above  the 
temperature  range  in  which  the  presence  of  physical  water  in  the 
samples  cause  an  expansion,  and  below  the  temperature  in  which  the 
sample  was  previously  heat  treated. 

8)  Application  of  "structology,"  or  the  obtaining  of  relationships 
between  structural  properties  and  physical  properties,  reveal  a general 
good  agreement  between  the  changes  in  physical  properties  and  the 
topological  transformations  that  occur  in  the  silica  gel  samples.  It  is 
suggested  that  the  first  region  of  variation  of  physical  properties  is 
related  with  the  region  in  which  the  pores  are  connected.  The  second 
region  is  associated  with  pore  closure.  For  instance,  as  the  connectivity 
decreases,  the  pore  coordination  number  decreases  and  the  stress 
concentration  decreases,  rising  the  strength  of  the  material. 

This  research  showed  the  importance  of  a complete  understanding  of 
the  structure,  involving  both  metric  and  topological  properties,  in  order 
for  the  behavior  of  materials  be  understood.  However,  many  additional 
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questions  and  problems  are  brought  up  by  this  work.  Important 
questions  and  problems  that  deserve  further  study  are  the  following: 

1)  The  full  densification  of  silica  gel  monoliths,  although  a feasible 
process  (as  shown  in  this  research),  still  is  a time-consuming  and 
elaborate  process.  In  order  for  sol-gel  derived  materials  of  optical 
quality  become  more  attractive  commercially,  more  research  is 
necessary  in  order  to  improve  and  decrease  the  costs  of  processing.  The 
initial  structure  of  silica  gels  seems  to  have  a profound  effect  on  the 
behavior  of  the  material  during  densification.  Therefore,  it  would  be 
useful  a study  on  the  effect  of  the  transformations  that  occur  during 
gelation  and  aging  on  the  structural  and  physical  evolution  of  silica  gel 
monoliths  during  the  densification  process. 

2)  The  type  of  atmosphere  inside  the  furnace  affects  the  densification 
of  silica  gel  monoliths.  It  is  of  practical  interest  to  explore  the  effects 
of  various  gases,  including  combinations  of  gases,  on  the  densification 
of  silica  gel  monoliths.  The  pressure  inside  the  furnace  and  the  partial 
pressure  of  the  gases  used  are  important  variables  that  need  to  be  better 
understood. 

3)  Silica  gel  monoliths  with  relatively  large  pore  size  structures  (>50 
A of  pore  radius)  have  potential  use  in  impregnation  and  doping 
applications.  Therefore,  a comprehensive  study  on  structural  and 
physical  (mechanical  and  thermal)  properties  of  these  materials  is 
necessary. 

4)  This  research  indicates  the  importance  of  topological  properties  in 
determining  physical  properties  of  silica  gel  monoliths.  The 
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introduction  of  topological  parameters  in  the  description  of  the 
structure  opens  up  opportunities  for  the  control  of  impregnation, 
doping  and  diffusion  applications.  It  is  suggested  that  impregnation, 
doping  and  diffusion  experiments  be  carried  out  in  silica  gel  monoliths 
of  known  topological  properties,  in  order  to  check  for  the  applicability 
of  the  proposed  topological  models,  and  to  derive  working  relationships 
between  the  structural  properties  and  processing  conditions. 

5)  Further  studies  are  necessary  to  explain  the  physical  mechanisms 
responsible  for  the  structure-property  relationships  observed  in  this 
research.  For  example,  a mechanical  analysis  to  describe  the  stress 
concentration  effects  associated  with  pore  connectivity  is  necessary,  as 
well  as  to  explain  the  correlation  between  the  local  levels  of  stress 
concentration  and  the  strength  of  the  material.  The  relative  and 
quantitative  importances  of  pore  coordination  and  curvature  of  pores  in 
the  local  level  of  stress  have  to  be  evaluated  in  order  to  clarify  the 
mechanisms  responsible  for  the  mechanical  strength  of  silica  gel 
monoliths. 


APPENDIX 


For  the  four  point  bending  the  maximum  deflection  (A)  at  the  center 
of  the  beam  can  be  given  by  the  tangential  deviation  (t)  which  is  given 
by  [533] 


A = tyu 


U 

fM  Xj 
Tr-rdx 


(A.l) 


where  tyu  is  the  tangential  deviation  between  two  generic  points  V 
and  U; 

M is  the  bending  moment; 
x is  the  position; 

E is  the  elastic  modulus; 

I is  the  moment  of  inertia. 

Considering  the  rectangular  beam  shown  in  Figure  A.l,  the  reactions 
Rl  and  R2  are  equal,  because  of  symmetric  loading.  Consequently  both 
Rl  and  R2  are  P/2,  where  P is  the  total  testing  load. 

The  variation  of  bending  moment  with  position  between  points  A and 
B (Mab)  is  given  by 


Mab  = j x* 


(A.2) 
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Figure  A.l  (a)  Schematic  of  a four  point  bending  beam  with  the 
cross  section  shown  in  (b). 
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The  variation  of  bending  moment  between  points  B and  C (midspan  of 
the  bending  beam)  is  given  by 


(A.3) 


where  L is  the  length  of  span. 

Introducing  Eqs.(A.2)  and  (A.3)  into  Eq.(A.l)  one  obtains 
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(A.4) 


The  solution  of  Eq.(A.4)  is 

_ 7 p L3 
A “ 384  E I* 


(A.5) 


The  moment  of  inertia  (I)  of  a rectangular  beam  is  given  by  [534] 
w t^ 

I = > (A.6) 


where  w is  the  width  of  the  beam; 

t is  the  thickness  of  the  beam. 
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Substituting  Eq.(A.6)  into  Eq.(A.5)  one  obtains: 


7 P L3 

~ 32  E w t3  • 


(A.7) 


The  elastic  modulus  for  four  points  flexural  testing  (Eb)  is  directly 
obtained  from  Eq.(A.7)  as 


_ 7 P L3 

" 32  A w t3  • 


(A.8) 
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